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§5. The correction to the current associated with 


the mass type term 

We shall now discuss the integral (37b) 
ape 2) =i% al (de )( de 19" (#)Pa BD (a= 2”)BBD oa? — 2") Dac 2H) 
+ $!(") BBD (2! —2!) PyDae( a" 2") BDx(2! —2) Bu) 
+L am{ (ar O)R BD lee ADAH O)BD A= ASE) (BT) 


which involves the effect of mass renormalization on the property of the mesonic 
_ charge-current. We begin with the evaluation of 


4 z K(2)=B,BD (2) B,Dor(*)- 57) 


_ Substituting the Fourier integral representations (22) and (35) for D,(«) and 
D,,(x) respectively, (67) becomes — 


L:@/) +L WB) —m— PB, 


p 
Sei ce i Ao taxa 2(k? + m2) a 


fe This. can also be written as 


BLi(9B—&B) ++ ( p8—AB)—m—— (P—- #8, 


4 
: BOP) Sto a (dk) (dp)e”" od 
mi: | (2 3s hL(p—ky +m) dis 


on introducing 
piewhatter (70) 


_ The denominator in (69) can be simplified by writing it as 


ve 
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(dul (pu)? + mtat + ( f+nb)u(l—n)}?. (71) 
0 
The transformation 
hy hy +Upy (72) 
brings (71) into the form 
{dull + mtu + (7° 4+ m’)u(l—x) |°. (73) 
0 


Now this transformation (72) as well as the integration over 4 must be carried 
out carefully, because the integral (69) is divergent. Thus, after the transforma- 
tion (72), a becomes 


— Gayl CO Canela uta + yt Pose ae) 


_ Polit 1) 28-48 + £1 (—a)pp—apy—m—= | 1) p—2 FI, 


(74) 
(Piet (+m) u(1—x) F 


Here the first and the second derivatives in the integrand represent the contribu- 
tions from surface terms that must be added to take into account the effect of 
the displacement (72) at large values of &, where the integrand does not tend 
to zero rapidly enough. Derivatives higher than the second vanish identically, 
since the integral is quadratically divergent with respect to the 4-integration. 
The integration over & in (74) can now be performed in a spherically symmetric | 
region of the energy-momentum space. 

Performing the differentiation with respect to 4, the integrand in (74) becomes 


(E +4) *{8—w) (28) += (1—u)*( 98)? (1-82) — mB 
+ BAB} + ud B+ #)15( 98) += Aw) (98)°—-82)} 
sre reicagen a? 
+2 -(A)*(E+ A) 2( 28)" + 978} LB) 


+ 2p") (B+ A) 9-34 (+ A) Hi i(L—w) (98( += (Ln)? (28)*(1—B,?) —mB,}] 
(75) 
_ where 
A= i! + (2° +m2)u(l—u). (76) 
In deriving (75), the following substitutjons 
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k, —>0 
duis: 
Ruky —-»> a0 om 
ky kk, —> 0 
kha —> 5 ()* (Byrd po DpDs0+ Duodoy) (77) 


have been made, in virtue of the dependence of the denominator of (74) on # 
alone. 


It is convenient to write (4°+/?)~ as 
(2 + meu?) + (2 +A)? — (2+ ne?) ~ 


= (+70)? —2( + m’*)u(l—z) [ase niates (P+ w?)u(1—u)|*. = (78) 


~ Then one obtains 


K(4)=- al (dp) >| dul | (dé) C28 wad 0 ie hy rE : (et epee = “ii 


x {i(1—n) (98) ++ ny" (98)*—Be) — mB + B18} 


+ liu 18) +501 (6—5x) (28) Be —z—0*9'8A—Ps)} | (79) 
using the formulas (57) and 
tee ent 
J) ea aC me 


Now, the integrand in the square bracket can be rearranged in the form of 
ascending powers with respect to 7#8+m as follows 


[J crey Ce mtaey-*{ —m (11) (=u) — mB 242 —u) + BAP} 
+E {— mu——-mu(6—Du) (1—£,’) +m hil i —#2)} | 
+ pp+m)| | (dé) (42+ m2) “| (L—n) (3 —2u) — 28,701 —u)*} 
+ FE (4 2u(6—5u) (1-P.)} | 

- (f° +m')| —2u(1—n) | (ae) (A+ mi) {—m(1—n) (2—u) —mB,2u(2—x) 


6m 


+H 2a —0)| — 2% wt920—B8) |+0(Gip8-+m)” (s1) 
2 7 
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where if8+m may be laid on the right-hand side of the square brackets. The 
integration with respect to 2 is performed in a trivial way for the lowest power 
term of (p°+’)-expansion. Therefore (79) becomes 


KG)= Sone (x) {eal \ (dk) (42+ mit) {— m(1—u) (2—u) —mB,2u(2—n) 


+ a B2(1 — B2)} + FE {nL me 6-51) (1—B,2) + pmt3,2(1—B2)} | 


4 : 1 : ave _ ae 
— Gai P+) 3) (ictal | (ae) (2+ mu’ {(1—u) (8—2u) —2p(1—x)*} 
+7 lu 2 u(6— 5) a-p,)} | 
4 2 9 1 . : r ee 
ent +m) (2) aul — (1-8, (1-88) | (@) (E+ mu 


+i {(-w)*@—w) + Aiu(1—n) 2—x) +4u(1-20)a (1 —A5)} | (Sm) 


The first term in (82) is one which should be cancelled by the counter self- 
energy term in (87b). In order to show that this is really the case, we shall 
begin with the consideration of self-energy. The electromagnetic self-energy 
of a free meson can be separated easily from the S-matrix, which is in the 
é’-approximation 


H"(2) = 19" (2)9(2) +8'@)9@)] (83) 
with 
d(x) = —inu{ (dt")B, BDp(x—2') 8, Dor(x—2')(2'). (84) 


This is treated almost in the same way as above with the only difference that 
we can now identify 78+ in (82) with zero on account of the fact that (x) 
satisfies the Duffin-Kemmer equation for the free meson. We can thus find that 


H™"(x) = ans) {a {(ae) (42+ mud)? —m(1—n) (2—u) —mB 2u(2—n) 


Be ¢ n° 
+ rer, (y— 48’) } + ey — mu — smi (6—5x) (1—8,°) + som (n—48,°) | Wc) 
a (85) 
ere 


- 6x1 for vector mesons 
=(-B) Bi| 
4x1, for scalar mesons Ned, 
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and 1, is a unit matrix in a #-dimensional space. Remembering the remark 


concerning the formula (26), we see that (85) gives rise to the same effect in 
the S-matrix as 


Ome" (x) (x) (87) 
with 
_ dtu (* 2 2, 2\ —2 : 2) ee 
m= Gayla JHA +s) {—m(1—n) (2-0) — 2 mu@—n)+ = @—10)} 
+ FE mm {ou— 2? + 4-10} | (88) 


which is obtained by replacing f2 in (85) with 5/2. What we want to show 
is that (87) can be used as the counter self-energy term of our theory inhte 
é-approximation. We shall write, for convenience, the self-energy (85) as 


H*(x)=Ap' (x) p(x) + A'P"(4) Boh (*), (89) 
then 
5 
In AA (90) 


The first term in (82), ie., 
4 (44 A'B2) (2) (91) 
Tu 


together with the. counter self-energy term then gives to the current d7:(7) the 
following contributions 


iF ae { (de!) (te?)| i) gD ale 2) L(A AAA HHCY 


+ ge) (A+ AB Va 2") BD —2) 8 | 


+ Ldml (de) PO) PeBD oe — 2) 02) +8" BDA 2) BHA) 


=£a'{(as| ea, BDde— 25 —Be e+ (er 282) BDde'—298.H0 |. 
(92) 

Now, since it is proved that 
pie! 5p \omB Dae —2)—id e'—2) = GOH Be AeA) BDole 2) 


(93) 
employing (1’) and 


(80+ m) BD, (x) =210 (2), (94) 
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(92) becomes 


x a'| (dx') [y" (x)B,0(4 —2’) (> a 82) (2) + $'(2’) (2 —F 7 )ae— *)Puf (#) ] 


=f. Aly (x) By (> Bp (z)+ w(a\(2- Be) Bus (=) |, (95) 
which is equal to zero on account of 


B.Be+ B22, =P. (96) 


This proof can easily be extended to include general cases. To see that it 
is really so, we shall observe a single meson line which consists of several 
internal meson lines and zero or two external meson lines and which has no 
vertex without a photon line. We now insert a graph of ¢-self-energy as well as 
that of the counter-self-energy term into this diagram in all possible ways. Then 
we obtain a series of modified graphs whose representative parts are described 
in Fig. 5. When the self-energy term is expanded into the ascending powers 
of Bd+m, the first term is easily found to involve the same factor 4+ A’8,* as 
K(x) (see (82) and (91)). 


(original graph) (modified graphs) 
Fig. 5. 


Therefore, we have only to show that we can generally determine the value 


of Om as A+ 2d! 


By analogous discussion to the one made in deriving (92), the modified 
graphs in Fig. 5 are found to give tise to terms 


POP emer meee eens 


cee. BDx(2'—x,)(2— B2)BDs(x,—2")B,.BD,(x"—2'") eh 
sree BD, (4! — 2") 8, BD,( 2" — x.) (>- 8,*)B Dy (%.—3!"") «+++ (97) 
piss cheese note tasixe | 


in the integrand of the matrix element in consideration, where only the terms 


involving the vertex x” are given explicitly and dm is equated to A+ ey, 


2 
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The integral of (97) is equal to that of 


eee eet osecceces 


ome To. 
+ {£3 = B)\BD,( x! —x"")B,BDx(2"”—2"") 
5 2 
++ BD,(x'—2")( = _ BP \BABD 2" 20") | 
i 5 
1 | 2p, (2! —x'’) Bul > — ,2)BDy(x"—2"") 


+4 BD! —2")8,BD ox" 2") — 8) | (98) 


eee 


using the identity 
benitt: on bt age: 
BD;(z'!—x) (2 Mi B2)BDr (@=2")-48@"— *)(> x B2)BD,(2—2") 
Ath i = bee 2 A 
~BD,(x +)(5 B, Jace 2!) 
1 9 9° 
pepe a ate P,P. —Pwp }BDe(4—2")] (99) 


(98) is easily found to vanish identically on employing (96) when our meson 
line involves no external line. When it is not the case, we must apply (93) in 
place of (99) to the vertex which is connected to the external meson line. The 
above conclusion is, of course, unchanged. This is the result which we have 
intended to prove. 

To make our statements complete, we shall give here the explicit form for 
the electromagnetic mass of the free mesons. By a partial integration of (88) 


with respect to ~, one obtains 


_ 2ta 1 <4 SrtA Toy, Se a0 42)-2 
m= od a 10) (de) (B+ m2) 147% (40) | (dk) (+m) 


1 
— 2intm| dul(2 +ut+ hart — 11’) 
0 


3 12 
+ nl du Qu — 0? + 414-10) | (100) 
2 Jo 2 3 
which is equal to 
af 1 4 K+K, 
— 2) - (y—10)( KK,—m | ° 
ales @ ( ea OB m 
m K+K, 
M (40) (1 Ae <1) 
tin” )(log — 
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m m (17 ) 
ay Cy tok pad pelle 
oe ee ae 7 


= 2m + 0-7) A484 £647) Loire he 1 @r-23) | (101) 


Ir me m 


using the formulas (55) and (56). When expressed in respective representations, 
the electromagnetic masses are therefore 


m= | KK 4 log Ath 3] (102) 
Ar uw m 12 
for scalar mesons and 
6m=— ml Ko +joghto_TT (103) 
Qi Ww 6 m 2 


for vector mesons. 


Since the term (92) has thus disappeared, only the remaining parts of (82) 
give contributions to the current 


Hi2 =) = BEI de") (de 18 (BBD ee) 9 +) (a — 2") PAVE") 


+" (2!) P(B) (Ba" + m) 32" —2') BD, (x! —x) Bp (*) } 


. 1Ue 


rea (dr')(de")1$" a) 5, BDAa—x') (—L"+ 9?) 0(2’ 2") Oe") 


+ $'(2") OCB) (—L)'" + m2) 6 (2" — 2’) BD; (2'— 2x) 8,9 (2) } 
ap cece ee eee cones (104) 


where we have employed the following abbreviations 
1 
PB) =P+ Piss] di {(1—n) @—2u) — 282} f (de) (B+ mbae) 
rd {ut = u(6—5u) (1—p,2)} 5 
2, 3 . é p } ’ (10! ) 


QO(2)=2,+ O.82=( auf — 1 u(l—x) (y7—48,° )| (dk) (2+ m2) ~2 


+= {( 1—x)’(2—u) + 8,°u(1—u) (2—x) +p0(1 —3)(y—42,)} | 


(106) 
After the partial integration with respect to 2’ the integration of (104) can now 
easily be carried out with the result 


Que 


(an yf (*) 1B. P(B) + PB) By —ByB(#) O(8) — Q(8)B(—x) 8, }6(#) (107) 
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where the formulas (24) and (94) have been employed. Bt(— x) means that the 
differentiation involved in it operates not on ¢(%) but on ¢'(x). (107) can be 
simplified as 


Gayl! OPP (2) QP. + 57,4 40,4 100) (108) 


on making use of the Duffin-Kemmer equations (1) and (1’) satisfied by the free 
meson field. In (108) only the first term of (107) survives. The divergent 
constants P’s and Q's are evaluated in a familiar way, the results being 
K+K, 83 

m 36/7" 


rnin Mog KKM), 


yesf t 
P,=i'(— log 


m 9 


GATE fT hogy Kot Ke 21 WA gt tee 
Dies pn MO en 3°72 


nt (2 
Q,=% (2 togh they 7 (109) 
m3 m g 

where 4&,,, is an invariant minimum light quantum wave number and therefore 
the term log(m/2kmi,) indicates the divergence associated with zero frequency 
light quanta, i.e., the “infra-red catastrophe”. Inserting these values into (108), 
one obtains the following result 

“ K+K, 

m 


(2 
Roe sr fea log 
ats 9) log 


+8 logo” + +. (713) bi,(2). (110) 
Bagels 

It is to be noted, however, that (110) is not the only result that can be expected 

from (104). In fact, there are many other ways of computing it. For example, 

we may integrate (104) after the partial integration with regard to z''’ (not to 

x’) obtaining a vanishing value for it. 

The correct procedure to treat (104) has been prescribed by Dyson™ and 
explained in detail by Karplus and Kroll’? in the case of electrons which can of 
course be brought unchanged into our case. The above way of computing (104) 
is the same with evaluating the modified function BD;(x) corresponding to the 
propagation function of a virtual meson BD,(*). Among the correction terms 
of the latter is a (divergent) term which is proportional to the original BD;(*) 
and therefore BD;(%) is expressed as 


BDL x) =Z,BD (x) $e (111) 


with 
Zo=1+dR , (112) 


where R is in a close relation to the coefficient of the current in (110). Now in 
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. . * . ‘ a 

order to presetve the unitarity of the matrix S, it is necessary to divide by Z, 

the matrix element of S between states coataining ” mesons, as Was shown by 
: : 1, 1/2 

the above authots. This is equivalent to putting equal to Lg = Leg) dae 


renormalization factors in eqs. (113) ; 


f' (4) =Z,'"9 (4) 
and 
f" (4) =Z)"°9"(4), (113) 


which represent the modification of the meson wave functions brought about by 
virtual interaction with the electromagnetic field. (113) is exact since g(x) and 
¢'(z) obey the Duffin-Kemmer equations for free mesons. (A straightforward 
calculation of ¢/() yie'ds of course a result with the coefficient Z,.) It is now 
evident that the correct expression for 072(%) is obtained by multiplying (110) 


with ‘1/2 since the second term of Z,'/?=1+ FaR+ involves a factor 1/2. 


Therefore, one finally obtains 
ye 2 K+ K, m 1 : 
d(x) = 12 7 log =~" +81 = —13)} 
jh (2) = SST —7) log =F Moy Slog s™— + (y—13)} iz). (14) 


§ 6. The correction of Lamb shift type 
We shall finally be concerned with the third cortection (37c) 


Oj; (x) == ue| (dx’) (2") f (2!) 8 BD, (a’—+) 8, BD (4— 2") Bib (2) Dax 2! — 2") 


(37c) 
or 
YE (2) = 17 ve (@x’) (dP) K, (2! —2, x— 2") (2) (115) 
with 


K, (2! —4, x—x") =B,BDj(x!—x)9,BD(4—2")B.Dop(x’—2") (116) 


which can be evaluated inserting the Fourier integral representations (22) and 
(35) for the D-functions involved. On performing this substitution, (116) becomes 


Ki, (4! —2x, x—2x'’) ek Ak) (dK) (dhl! ett AD 22) glk RINE = 21) 
[2° + 2*) (2+ mm) e) 
a(R! x LIAN? __ 1 2 a 4 
x YiCEB) + Wy —m— Lanta liwra) + wray— m—1vr}s, 
(117) 


It is convenient to replace Ay and &! by 


Pi=thy Pe=h tk (118) 


- 
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which is effectively subject to the restriction 
pitn=p'"+ m=, (119) 


since K,(#’—x, 2—2'’) is to be multiplied by ¢'(2’) and ¢(#”) and integrated 
with respect to x’ and x’. As a result of this transformation, (117) becomes 


K,(4'—2, x—+!') — au (dk) (dp’) (dp") giPl (AI 2) ppl (x—2!1) 


x (P+ mH Be me Ve 
x B, {i( 2/848) ++ (2'8—£8)* —m—— ('—2)} By 
Mm Mm 


x {i(p"—48) ++ (9 B—#8)'—m——(p'—A)'}8,, (120) 
mM Mm 


which can be rewritten as 


i 


Se “| (de) (dp) (ap! emer erro dead hf p +p! +u(p' —p") lu}-3 


‘G Ja 


xB {i(g'8—£8) +> (p48) —m—— (0-0) 


x {i(#"8—48) + © (p"8— apy—m—=(p" 8) 18, (121) 
since 
[PA 0H (pO +m }eP = [alu Ap +p" +009") Wl. 
J, (122) 


K,(#'—2, x—x") diverges logarithmically on integration with respect to #, and 
therefore it must be treated carefully. Thus the transformation 


Ry ky + 5 {be +put+v(pu—Pu) } (123) 
brings (121) into 
(— 27) d. ” fei a ip! (xt —x) giv! a—2!!) 
SN | (ae) (dp) Cap") te 
x[1— 24 pe Be! +0 pe — B12 | ft mires (0p ca—o)} "| 
p 
(124) 
with 
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@,=Ai{(1—% +2) (A) — (12) (9°8)} HB) —™ 
+1 (1-£049)(9M—£0-0) (0A) 
—1{(1-*a+0))9- £0-2"f | 
x Au[i{—% 142) )(2'8) +(1- 4-2) )(0"8) | #8) —m 
+4{-£(1+0)( 8) +(I- 20-9) )o")} 


—={-50+ v)~ +(1-% 1-2) "a, (125) 


where the derivative with respect to 4 in the integrand wepresents the surface 
integral contribution which is necessary to identify (124) with (121). Higher 
_ derivatives are unnecessary, since K,(2/—+z, x—2’’) diverges only logarithmically. 
To evaluate this integral we must first rearrange the. integrand into the form in 
which 2~’8+m is placed on the’left-hand side while if’’"84 on the right. Now, 
a left-hand factor 7f’8+m is equivalent to 0’+m operating on A,(2’—2z, 24-2"), 
which annihilates ¢'(2’) on integration by parts. Similarly, a right-hand factor 
ip''B8+m annihilates ¢(7’’). In peiforming this rearrangement, the following 
formulas 


Ou P (2) =B, (89) $2"), Ou" (2) =9,'$"(2") 88, 
(Oy + syydV +m3,)$(2") =0, Ons" (2’) +3,'$"(2’) yp —mg"(2’) 8, =0 


| (126) 
or their Fourier transforms 
Pu=Pv" BB, Pur=p'BBu, 
Put Surpy"” —imB, =0, Pu—Suvp,' —im3, =0 (127) 
are useful which can easily be obtained from eqs. (1) and (1’). Suv is given by 
Suv=Puby—Pr By. (128) 
After a‘ straightforward but tedious calculation, (124) becomes 
2)" aaneapyap'y| dof nd exp |i 242" 2) (y" he: 2. 
(Qn)# }~ diete so exp |i" p)(2— 7) 


x [ea maefi +257 +p")*(L--2%) i 
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s | #4, {10-8 —2n(1—n) +1 us, | 
m 4 


fp 


‘ (p,'—p.") ¢ 2 2 Q 
+7 Say a {2+ (9-27) + A2(—4+ 30) 


+ By{mQ—n)?+ (pp Ow) (Ian g “4 LALO 23) La) 
: oS (= 297210) 1—-)| 
timsy, (Df — 2) { 2 @-n)2+ Beae@—n) + (=0 (7 -S 6) 
+ {2—P (62-0) 04) 
— {14 LAP et) fmt Bae 29-8 +01") 
+ imsy(p) — 9.) 2 +38} | (129) 


In writing (129) we have exploited the symmetry of the denominator in (124) 
in connection with the & integration, and discarded terms involving #,’s in odd 


powers, while £,. is replaced by ete Terms of fourth power in & are found 


to cancel algebraically each other. Further, the terms linear in v are omitted in 
(129), since they will vanish on integration. 
It is now useful to introduce the new variables 


PEL Ath), =n"! (130) 


since the integrand of (129) involves px and p," merely as a combined form 
p,'—p,!'. Then the P-integration can be carried out immediately, yielding d(4/—+”). 
The divergent part of K(2/—2, x—+’’) can be isolated by a partial integration 
of the #-terms in the square bracket of (129) with respect to v, «as indicated by 


41 

9 -3 
fete 22 + mtu + F 1—x*)} | 
yu Am 


+1 

9 4 

=2( #2 + mie) 3 pute 4 mie + 0-” I (131) 
—1 


then 


K,(4'—2, #— es ae 3 (2! — 2") {(@) (dp) een adult (4+ mu?) 
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x i 10 87—2y(1—n) +2 (— 1 +~u)| 


ig, pforto-m +82-4490]] 


ee 6(a/— x") { (a) (dp) enon def ud + nie 1+ £0 —v) \y 


x | By{m"2— utp Q—a(1-u+ S)4% Late —3x)(1—v") + (FY 1—ml—} 


: 2 (5 s 
~ ims, pf 4 =u)? +> nu (2—u) +2°(1—x) a—0)(2—#e) 


*—x) G—2)} 


we {ry (1-7 | me "Biu(2n—8 +t z)- ints pa 2 — 4+ 28 a 


ae B(x! —x""){ (dk) caper duel nt + f0-“) Ir 


x "| B,{10—8u—2y(1—n) +2£(-1 + +1) 
—= sup, 2+ (0-27) + 82(—4+ 30)} | (132) 


The integration with oa to & anc w may now be performed with the result 


(a'—2"")| (dp)eme” 


K, (a —4,2-—2") =— fe y: 


M(Go-m +h Mates 
10-4 SFM Daft 


ae | ay) d ip(x—2!) | 4 loo m Ls 
tome - 4 aK Pe {a Bet = oe 2 +54 et 


9 


~ 


& | ve(14 L-*) 
Teac iG i bissinas 48 Hae) (428) 
An? 


; fel 1 a 
So ey ra > aa es? 
1% 4 6 1+ f4(1—e4) =i 12 7 —B,)) 


On the Interaction of Mesons with the Electromagnetic Field. Il 7163 


— = {8( 20-8 +2 Z)-2 —Syy p(e—n itn ya 


9 


pees 
d(x! -a")\(ap) Zt ly ee; 4mm’ 


re: ee) 
x[A( S244 Z)—ton a(G—2-Ae)| (133) 


On integration with respect to ~ and v (if possible) followed by an appropriate 


rearrangement, (133) becomes 


Ki’ 22, 4-2") = —Baw'—2A,l {> (Fito & KK a 08 a 
+4 -13)} ae") +2 ( Se +(4-2 1) @—2") 
em CE {(-!es oes Re) Ges nt “” (x—x)} | 
A (2-2) -2(2— ain (x—2')+ Sra A 


where F(x) is defined by (60) and 


— aay 


(4 oe (134) 


log(1+ 21-*)) 


Ha)=G ail el Po aes (135) 
Finally, inserting (134) into (115), we obtain 
g(a) =— {4 7-2) log k*% 44 log 3+ 37-18) 4) 
HN (—tog K+ Ms 4 2) Ti (a) — | (ee inl) eee 
+O(- log at ag) Pole *’) + F (2-2) -—3 (94> TA (e—#)}| 
ee a2, ee 2 Mo eee sit as” Jra(2) + {log X42 aah m2) | 


-# al Eb 2 {-3 (= iF (4-2) aaa pie dr") 


= 
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+12 AG-#')— Cp, (x2!) bn (x) | (136) 
12 Ant 
in whicn 
1 
My, (x) = _e f(r) —su (4) 
2mm 1 
and | 
ini (2) = —£-9'(2) 583 (2) (137) 
Qi 1 


are both antisymmetric tensors representing electric and magnetic dipole moment 
densities. The first term of (136) is proportional to the current operator and 
therefore can be removed by the procedure of charge renormalization. Actually, 
this is just equal to the current 67%(7) with opposite sign, as can easily be seen 
when one compares (136) with (114), and thus it is not necessary at all to 
renormalize the charge when these terms are added together, quite similarly to 
the situation encountered in the case of electrons.'?" 
(136) represent deviations from the current operator 7, (+), induced by the coupling 
with the vacuum electromagnetic field. They must therefore be observable if they 
have finite fixed coefficients. In fact, there appear logarithmically divergent factors 
multiplying m,,(*) and m,,(#) besides Jed Ju(#) previously discussed. The 
former also cannot be removed by the usual technique of renormilization and thus 
give another example of defects of the present meson theory. This is however a 
result expected from the beginning, as has been discussed at the end of §2. 
(136) has common form both for scalar and vector mesons which are only 
distinguished by the special representations of $-matrices (10-10 matrix for vector 
and 5-5 matrix for scalar) and further by 9 (which is equal to 4 or 6 for scalar 
or vector mesons respectively). To interpret (136) literally is, however, fairly 
dangerous. For example, (136) seems to show that the scalar meson has an 
anomalous magnetic moment, although its spin angular momentum vanishes. In 
fact this is not curious at all since the eq. (1) should be regarded as a wave 


\ 
equation of paitictes with spin 1 lying in the special state described by the 5-rowed © 
representation of the #-matrix. 


The remaining parts of 


To discuss the scalar and vector mesons separately, we must go from the 
Duffin-Kemmer formalism to the usual tensor forms (described by the Klein- 
Gordon or the Proca equations) by the well-known relations between the components 
of wave functions. We shall first consider the Klein-Goidon case. 


hi 
Ju(%) is equal to In this case, 


7e (Ou8* (x) -$(2) —9* (x) 0,9(2)), (138) 


where ¢(x) and ¢*(x) are wave functions satisfying the Klein-Gordon equation, 
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The tensor correspondents of dy,(v) and d,m,,(%) are both easily found equal 
to that of 


Er ; 
ape ny ANG: r). (139) 


When (114) and (136) are added together, they can be rewritten simply as 
Opn (#) + Ou (x) 
= £\ dey (w)| (442 D)ae—21) 
4a mr 


+log 5” —- i n@—ay+A@—2)}] ado) 
Bemin 700° 

Thus there appears no magnetic moment in the Kleia-Gordon formalism, Moreover, 

the ultra-violet divergences in (136) are those of mere appearance in this case 

since no such divergence remains in (140). 

Our result, (140), is marred only by the fact that there remains yet a 
logarithmic divergence associated with zero frequency light quanta, the so-called 
infra-red catastrophe. This is however not discouraging since one knows that the 
logarithmic dependence on 4,4, will cancel when all contributions to a certain 
3 scattering process are added together as was indicated by the work of Bloch and 
Nordsieck.” 

In the Proca case, the current operator 7,(4) is equal to 


ie{ (AuGs* (x) —O,bu* (x) )9,(4) —9,* (4) (Anh, (#) — On (2) Fs (141) 


where ¢,(7) and ¢,*(%) are four-potentials of the Proca field. Od, (2) and 
0,m,,(%) are replaced by 


1 
Tee 7. (x) +(1-3 us a My (4) +9,8eQuv0 (#) (142) 
and 
1 aps Xx 
5 LT ult )+(2 a, M,,(x) +29,0cQy,6(*) (143) 


where 7, (4) is given by (141) and 
Myx “ie eae (x) $,(2)) (144) 


Q ,* (4) box) —$o* (7) Aub, (4) ) (145) 


are operators representing the dipole and quadrupole moments of the spin 1 field 
respectively. Inserting (142) and (143) into (136), one obtains the foimula of 
correction to the current in terms of 7,(7), Myv(*) and Q,,.(%). Since no 
simplification is attained through this substitution, we shall not give here its form 
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explicitly. We only remark that jul), Mu(#), and Quve(¥) are all multiplied 


with divergent coefficients. 


§7. Discussion and summary 


In the preceding three sections we have analysed in detail the reactive besa 
arising from the interaction between meson and electromagnetic fields in ee 
second order approximation. The aim of this rather lengthy and labore’ 
investigation is, as was mentioned in the Introduction, more of an academic 
nature than practical. The obtained results even seem to accentuate this situation, 
since the theoretically observable correction to the original mesonic charge- 
current has turned out to be generally divergent and consequently failed to furnish 
us with quantitative conclusions about the electromagnetic properties of meson 
field, at least in the case of vector type. It is somewhat discouraging to find 
out that the Tomonaga-Schwinger theory, which brought about so great a success 
in the quantum electrodynamics of electrons cannot claim its range of efficacy to 
be very wide. But this is not unexpected. We know that the electron field 
with its Fermi statistics and first order character of the equation of motion is 
most favorable since the singularities of interaction are minimized. We should 
not expect the same thing from meson fields which obey- Bose statistics and second 
order wave equations. There may lie in this fact a déeper reason which is more 
than mathematical, though we cannot express it in a convincing way. 

But apart from quantitative considerations, we have obtained some qualitative 
and formal conclusions which are rather important. The program of quantum 
electrodynamics developed by Tomonaga, Schwinger, Feynman and Dyson can 
also be extended to the present case without essential modifications, at least from 
a formalistic point of view. However, a little more precaution and labor as well 
as minor deviations are necessary. First of all, the integrability requirements of 
the interaction representation necessitates us to supplement the original interaction 
Hamiltonian with a suitabie additional term. This causes a peculiar trouble when 
we inseit a Counter-self-energy term in the interaction, since the “mass” of the 
meson field is not an integrable quantity by itself and does not commute with 
the charge-current either, in the sense of Duffin-Kemmer formalism. We therefore 
can not bring about the interaction representation with counter-term from the 
original Schroedinger representation by two steps of transformation, namely first 


eliminating tlte free, field and next attaching the counter term, as was the case 


for electron, but have to do it in one step, thereby using for the free field an 
experimental mass instead of a mechanical mass. 
from the fact that there are two scal 


matrices so that we have to adopt bot 
effects. 


A further complication ‘arises 
ar quantities (1 and 8,°) built up of spin 


; h of them in order to eliminate self-energy 
This was avoided by making use of a relation by w 


me hich they are con- 
nected together wijth a divergence expression which can be 


counted upon as 


- —- 
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vanishing for free field. Translated in tensor formalism, the above fact is reflected 
in the appearance of terms quadratic in ¢ and 0,6. The latter, however, can be 
identified with the former by a partial integration, or alternatively, by a renor- 
malization of the Lagrangian. 

The cumbersome procedure of using correct integrable form of the interaction 
Hamiltonian with mass-term which is not yet given explicitly, is then replaced 
by a simpler one. The latter method consists in an appropriate definition of 
Dyson’s ordering operator P when two points under consideration lie in a space- 
like relation with an infinitesimal distance. In actual expression, we have only to 
adopt the interaction Lagrangian density itself as the interaction operator and use 
the corrected Green function, the first term of (23), which represents the vacuum 
expectation value of ordered field operators. Tlis has been suggested by 
several authors,” and is very interesting as well as instructive since it offers us 
possibilities for a consistent extension of the present localizable theory to a non- 
localizable one where the old Hamiltonian description would lose its validity. 

‘Making use of these manipulations, the program of mass and charge renor- 
malizations has been effected without difficulty, and the observable deviation from 
the original current has been determined. The charge renormalization is required 
only for that part which describes the polarization of the meson field due to an 
external electromagnetic field, as is the case for electron, at least in the ¢- 
approximation. When we go to higher order corrections, this may not hold in 


general. 
The actual result in the present case has turned out to be divergent. But 


for scalar meson this divergence is only superficial, since we obtain a convergent 
expression when we go back to the tensor formalism. Further, it has no spin 
magnetic moment in spite of its seeming appearance in the Duffin-Kemmer 
representation. For vector meson, on the other hand, we have essential diver- 
gences which cannot be removed. Moreover, there appears a true anomalous 
magnetic moment as well as a term which may be interpreted as a quadrupole 
moment. The latter did not happen in the electron because of the simpler 
structure of the Dirac’ equation. All these terms, however, are multiplied by 
divergent factors, which makes it difficult to come to a conclusion about ex- 
perimental verification of such effects. To relieve ourselves of this point, various 
devices may be suggested. For example, we may start with an expression for 
zeroth order current in which all these anomalous terms are contained with un- 
known factors yet to be determined so that all divergence may cancel out in the 
final result. This, however, fails since such an inclusion of anomalous terms give 
tise to still higher divergences in the next step, as was revealed by Feldman,” 
so that it must be continued ad infinitum until at last we reach a meaningless, 
ot at best an integral rather than an differential, expression for the current. 
Another more effective means is furnishd by the Pauli regulator, or its more 
realistic version of Umezawa and Feldman, With suitable conditions, we can 
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certainly scratch all undesirable terms from the result. But again sn shall have 
no excuse for its highly arbitrary character, so long as we can give no further 
theoretical or experimental foundation about it. Therefore we steal? be hrestes 
content only with pointing out these difficulties as they are, hoping that this 
may help towards their correct and thorough solution. 

The prospect over any application of the present result to practical problems 
is thus rather poor. As is well known by Christy and Kusaka’s work on burst 
production the interaction between vector meson and radiation field becomes 
increasingly large at high energies, yielding an essential difference from scalar 
meson.’ .But at these energies, the higher order corrections are in general also 
large, and we may expect that fhe interaction does not lead to a divergent 
catastrophe at extreme relativistic region. In reality such corrections are by 
themselves divergent, and no definite refinement of Christy-Kusaka’s results seems 
to be possible for vector mesons. 

We have, however, more conservative ways of analysis, as will be seen below. 
Terms are divergent, but we know their physical properties as well as the effects 
they would cause if the infinite factors multiplied on them were actually finite. 
The order of magnitude of these “finite”’ coefficients may be determined by 
analogy to other really convergent cases. In some cases we could even predict. 
their sign from the divergent coefficients. With such convention, we get for the 
charge-current of a vector meson a picture which consists of an elliptically ex- 
tended: intrinsic charge-current and a spinning current. The former includes a 
quadrupole moment of the vector meson. It would be interestiag if such effects 
could really be observed. Their magnitude, however, would be too small, and 
the vector meson itself too short-lived, if exists at all. Hence the only possibility 
may be found in the deuteron, in which vector mesons might play a role and 
contribute a magnetic and a qurdrupole moment. But it may be premature to 
enter into such minute effects in the present stage where more basic facts about 
the deuteron remain to be explained. These anomalous terms of the meson 
would also cause an anomalous behavior of a single nucleon with which it interacts. 
But it would again be a higher order effect, even if it turned out convergent. 
In this connection it may be noted that a single nucleon cannot have an 
anomalous quadrupole moment so long as it obeys the Dirac equation. 
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§ 1. introduction 


Previously, we? synthetically analysed the problem of vacuum polarization 
by considering simultaneously various charged particles. By this: analysis, it 
becomes clear that the problem of gauge-invariancy of vacuum polarization closely 
connects with the definition of the vacuum in quantum field theory.” It isthe 
purpose of this paper to analyse this problem of the definition of the vacuum under 
more comprehensive version and to investigate the problem of gauge-invariancy 
and various divergencies in the present quantum field theory. In order to make 
clear the physical meaning, correspondence theoretical argument is used in this 
paper. This argument provides some interesting suggestions on the nature of 
vatious difficulties of present quantum field theory. Thus we will propose a 
correspondence theoretical argument relating to the vacuum and expect in future 
its complete formulation in quantum field theory. 

Recently two important progresses were made in the region of quant'm 
electrodynamics, one of which is the completely relativistic formulation of the 
theory,2® and the other the classification of types of divergencies on the basis 
of this formalism.?® It was clarified by means of the above mentioned method 
(or by the improvement of usual perturbation calculation) that the self-energy of 
an electron has the transformation propetty of the mass and that no diverging 
term except the mass-renormalization term and the charge renormalization term 
appears in the system of electrons and electromagnetic fields. On the basis of 
this fact, the renormalization method (Self-Consistent Subtraction Method) was 
proposed in order to evade the divergence difficulties and to obtain the finite 
terms to be compared with the experiments.” The essential point of this 
method is that one may put divergencies of the mass type and of the charge 
type into the observed mass and charge of electron respectively, as the alterations 
of the mass and the charge are, if they are finite, to be already included in the 
observed mass and charge. This method led, in fact, to the Lamb-Retherford 
shift of the energy level of hydrogen atom which agrees almost with the 
experiment. However, there still remain many problems to be investigated. 
It must be, at first, investigated whether this renormalization method is valid 
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not only for the system consisted of electrons and electromagnetic fields, but also 
for more general system. More important problem is to analyse the character . 
of divergence difficulties evaded in this method. It is demanded to make efforts 
to disclose the internal conflict in the present theory of elementary particles in 
the course of its adjustment of various difficulties and to grasp the correct laws” 
of motion in future. Especially, various divergence-difficulties appearing in the 
theory of elementary particles seem to be related with the structure of elementary 
particles. 

As the first step of the analysis” of the structure of elementary particles, 
hitherto, we have had Prof. Sakata’s view point, who suggested the synthetical 
studies on the quantum field theory by taking account of various fields simultane- 
ously. The theories”?® developed in this direction have revealed various features 
in the theory of interactions between elementary particles. From these “concrete 2 
theories, it may be possible to construct ‘‘ abstract” theories by abstracting 
each of these features. (Of Course such a abstract theory is frequently of transient 
character but, by means of this theory, we may be able to investigate the 
difficulties from the more extensive standpoint). For instance, the one aspect 
of. contents of ‘C-meson theory for solving the divergence-difficulties of the self- 
energy of an electron which has the transformation property of mass term was 
abstracted into the self-consistent subtraction method.” Such a theory, that was 
established by abstracting only one feature of a concrete theory, however, cannot 
‘explain all difficulties. The renormalization method fails, in fact, in the problem 
of the self-stress of an electron.” Through the way to construct many abstract 
theories by extracting various contents of the concrete theory, we can make clear 
up the physical character of these difficulties. Recently, according to such a 
method of analysis, we mentioned that one characteristic of mixture of charged 
particle field which we? and also Feldman” and Rayski®™ introduced in the 
region of the vacuum polarization, leads to Pauli’s regularization method and the 
other aspect is related to the problem of the definition of the vacuum2 “In this 
paper we will discuss in details this second point and show that it is possible 
to petlone some new important analysis on the problem of gauge-invariancy and 
divergencies from wider field of view based on. this argument relating to the 
vacuum. Then, we will propose this argument on the vacuum in quantum field 
theory, independently of validity of the mixture-theory of various charged particles” 
from which we stafted. In this paper, this argument is developed correspondence 
theoretically and we will in future expect its complete formulation in quantum 
field theory. 
Si boat i uaa teeing ae Rea 
vanishes in future correct theory ax iit nes Bere en aie persis: 
x(x*=0), if its self-energy have eee? wid wee: he ae ad 

nergy have a term which does not vanish when x0, this 

term may vanish as in the case of the mass term of a photon in future correct 
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theory.* Thus, in the case of various fields, non-zero mass terms which must 
vanish in future correct theory may remain in the calculation of present field 
theory. We may call such teims as ‘“ generalized non-gauge-invariant terms ” 
((g.n.g.) terms). The (g.n.g.) term may appear in all problems which relate to 
the vacuum, e.g. the problem of self-energies of elementary particles, the problem 
of meson decay through the medium of vacuum particles, etc. 

In this paper, the correspondence theoretical method is proposed which leads 
to vanishing (g.n.g.) terms. In §2, We will investigate the physical origin of 
the self-energy of a photon and give the method in which it vanishes. In 
following chapter, by means of this method, problems of self-energies of elemen- 
tary particles and meson decay** through the medium of vacuum particles are 
investigated. 


§2. The self-energy of a photon 


Denoting the self-energies of a photon with the momentum due to various 
charged fields by W, the mass term of photon is appeared in the Hamiltonian 
as AS so that, if Wis not zero, it destroys the condition of gauge-invaiiance. 
Previously, we had calculated the self-energies of a photon due to the various 


charged fields (mass #) and obtained the following results” : 


+1 for scalar charged field, 
m= 4—2 for spinor charged field, (2) 
+3 for vector charged field. 
From these results we” had proposed the following fields-mixture for eliminating 


the self-energy of a photon: m scalar, x spinor and / vector fields with masses 


pf, # and p" respectively : 
—2n+m+3l=0, 


(4) +3 (3) 
— 23) (pi)* + 3 4a)" +32) (8) 2=0. 


* In the case of neutral scalar field, too, by means of introducing mass coordinate in the five di- 
mensional representation, the condition of gauge-invariance can be introduced. 
In this case we must care about the essential difference between the electromagnetic field and 
the neutral vector field with non-vanishing mass. cf. F. J. Belinfante, Phys. Rev. 75 (1949), 
337. 

** In this paper, system of unit #=c=1 is used. The spatial component of four dimensional 
vector £(4,, 40, kz, 24=iko) is represented by bold type letter A. (20), £@)) and (KM, K)) 
means scalar products in four dimensional space and three dimensional space respectively. 
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These conditions for mixture fields have two interesting characteristics. One is 
that the second formula of (3) is the relation between masses. This feature 
gives the realistic meaning for Pauli’s regularization method. The other important 
aspect of this mixture theory, is that the self-energies of a photon due to various 
charged fields are in the same proportion and of the same signs as the zero-point 
energies of these charged particles : 


ieee SEE 
Exo=n| VP+e (xy , (4) 


where 2 has the same meaning as in Eq (2). It gives some suggestions for 
the problem of gauge-invariancy to abstract this aspect of affairs from the «mixture 
theory and to analyse it. 

In order to carry out such a analysis, we return to Dirac-Heisenberg- 
Weisskopf’s positron theory. Dirac” and Heisenberg? used the density matrix 
R, and the Hartree approximation in their calculations. When ¢(2) means the 
wave function of the electron, R, is: 


(RI = TA a* @) fale”) — Bh ba 2" Pa 2) ) 


>} or 5) means the summation over the states which are occupied or unoccupied 
occ, unocc, 


by the election, respectively. They assumed, that the nonsingular quantity, which 
is obtained by subtracting the singular quantity S from the &, with singularity, is 
the actually observed quantity : where S is deteimined such as: R,—S=0O in vacuum. 
Therefore, S is the charge-density in vacuum: consequently, the determination 


of S is equivalent to defining the vacuum. Dirac and Heisenberg obtained the 
following result for S: 


wx v x 
(2'|S|2"’) adap SE ita) + w log eal (2=2'’ —x''), (6) 


wu, v and w are the quantities without singularities. Especially, 


z 


nome (tf ts) ® 


in which the integral means the line integral along the straight line joining two 
points + and g’. 


Thus, Dirac and Heisenberg constructed the gauge-invariant positron theory. 
The physical meaning of this subtraction theory was given by Weisskopf.” 
According to Weisskopf, 


to subtract (6) from &, is equivalent to following 
subtractions : 


(1) Subtraction of the energy of the vacuum electrons in free space. 


(2) Subtraction of the charge and current density of the vacuum electrons 
in free space. 
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(3) Subtraction of the chronologically and spatially constant polatizability 
of vacuum. 

Subtraction (3) cotresponds to the charge renormalization, and (1) and (2) to 
the definition of the vacuum. Then according to. Weisskopf, however, it is not 
adequate to use the ordinary zero-point energy as the energy of vacuum electrons. 
The factor w appearing in (6) suggests that one has to use the following quantities 
as the energy of the vacuum: 
| 2 
(zy: 


The factor xz, as its form shows, is important for obtaining the gauge-invariant 
result. Accordingly, the definition (8) of the energy of the vacuum is essential: 
for constructing the gauge-invariant theory. From the definition (8) of the energy 
of the vacuum we get following postulate : 

(1) As the self-energy of the vacuum we must take one in which J, is 
replaced by £,—eA, instead of the ordinary zero-point energy; more generally 
we must defiae the vacuum with p.—eA, and not with ~,. But to define the 
vacuum by (8) is only the correspondence theoretical argument, without considering 
the second quantization of electromagnetic field. Accordingly we will investigate 
the physical meaning of (1) by using the classical theory. 

In classical dynamics, for the charged particle. with momentum P and velocity 
v, we have following relation: | 


TG, ee 


|v Pa eay tp eA}. (8) 


VP-ed) +f —ehy= 7h, 2 


i.e., postulate (1) specifies the state of charged particle by its velocity v. More 
interesting fact is that since the ight side of (9) has the transformation propeity 
of the fouith component of the energy-momentum four vector, W%?+4 72 has not 
such a transformation property ; accordingly, it is natural that in the usual theory 
defining the energy of the vacucm by ordinary zero-point energy, only the spatial 
Component of the electromagnetic potential contributes to the mass term of a 
photon and gives therefore the relativistically non-invariant result, and also that 
the usual definition on the vacuum (4) provides non-gauge-invariant result, 
because in Dirac-Heisenberg’s gauge-invariant subtraction theory, (8) is used as 
the energy of the vacuum. 

It is also possible to make clearer the relation between the self-energy 
of a photon and the definition of the vacuum. By means of postulate (I), the 
energies of the vacuum of various charged fields are defined as follows: 


Lex crch {¥(P—¢cA)? +2 —ed,}aP. (10) 


Expanding (10) into power series for the charge ¢, as in the peiturbation 
calculation, and taking the terms up to ¢”, we obtain?; 
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Uj= (Exero— Ap) a U;: ¥ 


; 11) 
ne-{ 1 ve ( 
U.= | Sen ie 

w=(P, P!) =2P°?+=E£p— (P, -P). ql’) 


U, is just the mass term (1) of a photon. When we use, therefore, (10) as the 
energy of the vacuum to be subtracted in the calculation, the self-energy of a 
photon vanishes and gauge-invariant vacuum polarization can be obtained.* 


* Tt is well known that the value of the self-energy of a photon is very ambiguous and depends on 
the method of calculation. The result (1) is obtained when integration over energy is firstly 
performed in the calculation in four-dimensional energy-momentum space. This definition of integra- 
tion calculation gives (4) as zero-point energies. Therefore, it is adequate to compare (11) with (1). 

Previously we pointed out®!6)17) that above-mentioned relation between zero-point energies of an 
electron (that is W/cox) is established independent of definition of integration calculation. 

Pauli!®) and Wentzel5) calculated the self-energies of a photon by using the integral represen- 
tation of 4-function. 


Fe(ey ee fa (k,2+p2)+tke, (i) 
cae (aay [a gota (ky + 
Then, we gbisin as the self-energies HW” of a photon due to various charged particles by using (i): 
Wim j= oe ) * AtAy |v Je] o- 
Wea Be ag | ae OC Carden) + a8Bay fete ES (ecm =l), Gn 
where 
Cu= fete" pate (dk). (iii) 


Cyuy is ambiguous integral. Pauli and Wentzel obtained finite value of the self-energy of 
a photon by means of following definition of Cyy which is obtained from (iii) by partial integra- 


tion: 
Cust — day (fis yr+P*) (7) (iv) 
= —d py a(t , (v) 
(J (aayeatyr— Te @)). 
Then, 
2 3 
Wm sare ay wd) 


=-2 a fo 

For 2 we get W= 8.3 Ay? as as mass term of a photon due to electron which is the 
result obtained by Pauli and Wentzel. 

We will show that the postulate (1) make self-energy of a photon zero also in the case of 


the above definition of integration calculation. The zer arged 
ls o-point energy £ of scalar ch fie 
U is represented by ae 4-function as follows: ne ; 


E soro= =-\(ve Lag a tu ‘| Vo)). dx 


ae J 0 . 
Ar Jace -? ieee Ove “at AOR 2) dx 


— amtdg | (D420 but + 22) 
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9° 
=~ Gate 12a] (aA Attorney, 


Similarly, zero-point energies Zyery of various charged fields are as follows: 
—2n 
(27)347 
By partial integration, as in the case of (iv), we get: 


E; zero = 


Jaa (a) éteteCrye +2, (vii) 


Fy aeons Jas OEY Pao teace. 84.48) a 
Szeru = =. ak df ave 2 tH), (viii) 


According to the postulate (1), we obtain Uo as the energy of the vacuum by replacing Zp in (viii) 
by 4y—e4yp. By expanding Up into power series of é, é-term vanishes by integration calculation 
and we get ¢*-term U2 as follows: 


ine? da 
Cis \ db \ aAu2—2a2( Ah) 2Yete (hy? + 22) 
P= — ayaa (1A) | Tiny? —2a( 4A) Moto ‘ 
ine dan. Ee Ae oa , . 
= aejiag | OD J 2at hak, ats 2B Jt 0 dy (i 
By using (viii) and transforming g—> = , we obtain U, as follows; 
# 
ne* ¢ dw 7c%e (7) . And 
O=— z lepeee (20 -—C 2 } ew = 
2 e(w) eit (2Cyv—Cy,dpv) + w2dpvt (4z) 


"Tae? (Ae)? oY 
which agrees perfectly with (ii). Therefore, the self-energies of a photon can be made to vanish 
by postulate (I). On the other hand, Schwinger) and Furry®) obtained the zero self-energy of a 
photon by means of suitable definition of integration calculation. We can show that their defini- 
tion of calculation is very one that leads to vanishing ¢242-term (the mass terms of a photon), 
which comes from the energies of the vacuum defined according to postulate (I). 

From (4) zero-point energies are: 

2 uty 
 (27)8(4n) 4 
According to pestulate (1) we obtain the energies of the vacuum by replacing 4p by 4p—edy in 


F/zer0= 


| (ae) p26 hu? + 2). 


Lrero- 
2 ; 
8((Au—e4yp)?2+ p2) = B(hu2+ pw?) eA yb, (hu? +2) + Ay AuBpv (du? + a?) + stamens , (xi) 
where 


0 > ° 9 
by (hy? + p*) =p, 0 +H"), 
v 
2 


° ° ee Doe ey} 2 
00 (x) 


7a) = ae 


The ¢2-term U_ of Uo is, by means of (xi) as follows: 


2 
U,=— ( ne vi j (42) [24p4, hyd, (Ap? + 1") 
2 9 9 9 ie 
+ ha Ay Avis ha + uP) + Awd lv? + Hw) 
By partial integration, we obtain: 
2 Pea 
Uemra vag (dh) [28 (dn? + 2) Bp yt Mako” a+ HP) “Cy (xii) 
=! 


(xii) agrees perfectly with the mass terms of a photon obtained by Schwinger and Furry. They 
obtained vanishing self-energy of a photon from (xii) by partial integration. 
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However, postulate (I) is the correspondence theoretical definition and so we 


cannot directly proceed this method in present form into the quantized field 


theory. A photon, introduced into the charged particle fields, make increase the 


zero-point energy of every charged particle up to Evero (because Eyer is defined 

by.méans of momentum of charged particle in interaction). So that, we rewrite 

the postulate (I) as follows: 

(I ) The zero-point energy is deviated from the energy of the vacuum on 
account of the existence of interaction. 

This definition of the vacuum is required in the calculation for interacting 
field and (I’) requires to define the vacuum in considering the existence of the 
interaction.* « 

In order to perform the correspondence theoretical investigations in various 
cases, moreover, we set following postulate :~ 
(1l’) The deviation of zero-point energy from the evergy of the vacuum 

is determined by means of the deviation of the equation of motion 
of interacting field from that of free case. 

Postulates (I’) and (II’) are reduced to the postulate (I) in the problem of 
vacuum polarization in quantum electrodynamics.** 


§ 3. The problem of gauge-invariancy 


We will study the more general problem of gauge-invariancy than that of 
the self-energy of a photon. 
a) System of charged field and electromagnetic field 

We obtain following terms as the energy of the vacuum by the st:ict 
integration of (10) ; 


2 ° 9 $ 
T= evel (Ped) 1} —¢A,\dp 


n wren oT 6 .\ ¢ A? n 
See Ph AY ane 3 he) Beh OE = 
al” er aa ms ea 2 x 60 


-“(44)(4A).  (10’) 


The second term is of aourse the mass term of a photon (11). 
The third is the term which contributes to the process of photon by photon 


* This fact seems to show that we cannot define the absolute vacuum independent of the interac- 
tion, but must consider about the concept of vacuum in considering the interaction, Since the 
relativistic covariance and the gauge-invariance of the theory of the system consisting of 
interacting fields is established in considering the interaction, customary definition of the vacuum 


independent of the interaction may destroy the various invariancies of results of calculations in 


problems for interacting fields, 
** Mr. Machida proposed more quantum theoretical method of treating these postulates (I’) and 


y i 
(II’). In his argument, the field operator in vacuum is the Stationary-state-solution with fixed 
7» (momentum and energy) of its wave equation including interaction. 
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scattering. This term perfectly agrees with the non-gauge-invariant term appealing 
in the case of the photon by photon scattering through electron pairs (7=—2) 
in perturbation theory which Euler found. Therefore, in this case, we can 
obtain gauge-invariant result by means of postulates (1’) and (II’). 

No terms of higher order than the fourth in charge (terms coe”, m>4) appear 
ime 10’) .* 
b) Decay of a neutral meson into photons through the charged spinor particle 

pairs 

As the interaction between charged spinor particle ¢ (mass /) and nevtral 

scalar meson g (mass x) we take the form: 


SE*BYE. (12) 
Then, according to the ordinary perturbation calculation in the momentum space, 
we obtain the following term as the non-gauge-invariant term :? 


He fp A’. (13) 
67? 


But in this case, we define the energy of the vacuum by using (I’) and. (II’) 
instead of Eyero. According to (II’) we write the equation of motion of ¢ as. 
( 


follows : 

{ (a, p—eA) + Bp +fBy}p=0. (14) 
This is just the free field equation in which # is replaced by p+/g. According 
to (I’), as the energy of vacuum we take 


Fear ek V (P—cA)?+ (Utfey—Aiap. (15) 


Expanding into power series of e and f and taking the term of the order é*f, 
this term agrees completely with (13). Therefore, defining the energy of vacuum 
by (15) according to postulates (I’) and (II’), the gauge-invariant result can be 
obtained also in this case. 

When a neutral vector meson U, decays through the nucleon pairs by the 
interaction f¢*7,~U,: it is suitable to define the energy of vacuum as Evero in 
which ~, is replaced by 

st Ga, SHOR! (15) 


Expanding this energy of vacuum, éf-terms vanish by performing the integration 
over angle, so that, terms corresponding to the case of neutral vector meson 


ef _(A%) (AU) which is 


decaying into two photons disappear. Term of order e*f is bn? 


* Accordingly, in the case of the mixture theory (3) or Pauli’s regularization method, we may 
use ordinary zero-point-energy “zero as the definition of the energy of the vacuum, 
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just non-gauge-invariant term in the case of neutral vector meson decaying into 
three photons. Thus, it becemes obvious that the problem of gauge-invariancy 


connects with the definition of the vacuum. 


§4, The self-energies of elementary elaine 


The postulates (I’) and (II’) would lead to different results from usual ones 
in problems connecting with the subtraction of vacuum. One of such a problem 
is that of self-energies of elementary particles. In this case (1) and (iI’) have 
following physical meaning: the self-energy of elementary particle a due to the 
field 6 is the additional energy to the vacuum that arises when a is brought into 
the vacuum of 4. This consists of the zero-point-energy Ey... Which comes from 
the zero-point oscillation of 4 forced by a@, and of the increment of energy of @ 
due to 4. In order to subtract the deviation of energy of @ due to the forced 
zeio-point oscillation of 4, one must use the postulates (1’) and (II’). It was the 
most simple example that the self-energy of a photon vanished by means of (1’) 
and (II’). Now, we reexamine the self-energies of the neutral mesons and of 
the charged paiticles by using (1’) and (II’). 

(a) The self-energy of neutral scalar meson (C-meson) 

(i) Self-cnergy due to the charged spinor field (mass p)* 

Adopting (12) as the interaction, the mass term of neutral scalar meson 
(mass x) is 


° 


Wor =— ee p +£( um —3e ‘log =¢ { + finite term} ¥ (16) 
pre (47) 


in usual perturbation calculation.” On the other hand, the energy of vacuum in 
this case is defined by (15) according to (1’) and Cl’). Expanding it into 
power series of ¢ and f/, we obtain as the f*-term 


i Pk See 2p ©) »,2 ‘ 
-{£ Lov pte —tL hog * +22 fl ae (17) 


This agrees completely with Wo». (x=0) including the finite term. 


The condition x=0 seems to show that Wop, (x=0) is (g.n.¢.) term. 


Thus, when definixg the energy of vacuum also in the case of x0 by Eq. (15), 


the mass term Woy of neutral scalar 1 
é neson due to tl 
by subtraction (17) from (16) 1e charged spinor field is 


2 £8 9 2 
Woy, = — | “f log =P + finite term. Me AR (18) 
r veo (4r) 


* Jn following, we designate the mass of charged particle as # and that of neutral particle as x, 
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W'cx, (x=0)=0, (18’)* 


which has only logarithmic divergence.** 
(ii) Self-energy due to the charged scalar field U(mass p) 
We take following form as interaction : 


aah lt. U*Uo — fu Ue’, (19) 


then the mass term Wes. of neutral scalar meson is” 


9 


Wos.= {-S ov Pre +#(£_F) log =e ih bo eee (20) 


according to usual perturbation calculation. While, the equation of motion of 
the U-field is 


{O-#+ rete" }U=0, (21) 


ie., one that #2 is replaced by (#—pfy—g’¢g’) in free U-field equation. Then, 
according to (I’) and (II’), we define the energy of vacuum of U-field as 


ally aaa eae Pe eae 2) 


Expanding it into power series of ¢, f and g, and picking up f*- and g®-terms, 
the result is 

PA OMY Seay #(£-L)1 g A nVae 93 

on? bape hes 2 4 ter pro (470) 2 
In the case of x=0, this agrees perfectly Wes, in (20). In the case x#0, when 
we define the energy of the vacuum by (22) according to (1’) and (II’), the 
mass term W'g5, of neutral scalar meson due to the charged scalar field becomes 
finite by means of subtracting (23) from (22): 


i ee 
Wros,= (f-t) Gin)’ (24) 
W'os,(x=0) =0. (24’) 


* Prof. Sakata pointed out that we may possibly select (g.n. g.)-terms which appear in the 


ordinary calculation as those which are independent of the mass of the real (not vacuum) field 


at rest. 
** Recently, Y. Katayama,””) J. Rayski!!) and P. T. Matthew’s®) independently obtained the self- 


energy of a neutral scalar meson of which order of divergence is reduced to the logarithmic 


by means of suitable definition of ambiguous calculation, But this self-energy contains stil] 


(g.n.g.) term, 
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(b) The self-cnergies of charged particles (mass [4) 
(i) The self-.nergy of charged scalar particle U duc to the 1. eutral scalar 
y ficld 
Adopting (19) as the iateraction, the equation of motion fer n2utral scalar 


field ¢ is 
{T]—-# 422°U *U }yp=fu *U. (25) 


The self-energy HH’, of charged scalar paiticle 7 due to the interaction (19) is, 
according to O. Hara,” as follows: 


1 


U*U 
Wa 


loa (27) ; 


| oP + (tft—xe%) log 24 + finite term} (26) 
x 


According to (I’) and (II’). however, one must take into account of the altera- 
tion of the energy of the vacuum due to the deviatioa of equation of motion (25) 
from free one. In Eq. (25), the term fuU*U is added to equatioa for fiee case, 
besides the replacement x°->(x?—2¢°U/ *U). Since this term does not involve ¢, 
it cannot be treated as the change of the factor of ¢ in the equation of motion 
of g. It is obscure how to introduce this term into this correspondence theoretical 
argument. For a moment, we will, in following, adopt the new postulate as 
follows. 

(IlI’?) Terms that cannot be regarded as the change of the factor of vacuum 
| field quantity in its equation of motion are left out of consideration. 
and we will postpone the investigation of sech a term in future. Then, taking 
into account of the presence of the term 2¢°(7*U in (25), one must, according 


to (1’) and (II’), define the energy of the vacuim of neutral scalar field ¢ as 
follows : 


1 eS ae 
fal VE 422080 Mt ) 
Expanding it into power series of g, and taking out the g*-term, we obtain, 
it eye: 45 ee G*U 
——— | o 2 x2 e-? low —o°x° cali 
ga SEOs og 4 tel ~~ (28) 


Using the postulate’ (I’), (1I’) and (III’?), the self-energy IIs of charged scalar 
particle due to the neutral scalar field has only the logarithmic divergence : 


U*U 


Wis={}—EP" ely 
iso (27) 


2/ : 
/. log + finite term} 
z x 


(29) 
(ti) The. electromag netic self-energy of a charged scalar particle U 


The equation of motion of the vector potential em 


: ; of the electiomagneti¢ 
field interacting with the charged scalar particle is; ; 
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(T2220 *U ) A=ie(U* grad U—U grad U*), 
(14,=4aic(U*M—IT*V), (30) 


( = : 
T= Oe AAUP); 
4x\ dt 
where J] is the canonical conjugate quantity to U. According to O. Hara, the 
electromagnetic self-energy of charged scalar particle is 
3¢ 9 9 coe 58 
WS aelhs + 6 log (4) +finite term} hath (31) 
An pL i>o (27) 
From (30) we see that the alteration of the equation of motion due to this 
interaction appears in the three components of transverse and longitudinal waves 
and no change in the time component. Therefore, the energy of vacuum is 


F Gal VF eT Na (32) 
from the postulates (I’), (II’) and (III/?). In the case of ordinary zero-point 
energy, the contributions from the longitudinal part and from the time component 
cancel out each other. But in Eq. (32), this is not the case because only the 
oscillation energy of longitudinal wave changes by the interaction with the 
charged paiticle. Expanding (32) into power seiies of ¢ and taking out the 
e’-term, 


tae ve 8) 


Accoringly, the electromagaetic self-energy W’'.... of charged scalar meson based 


on (I’), CI’) and (III’?) is by subtracting (33) from (31) as follows: 


Wa. = SEF log ! + finite term) ooas (34) 
An aa pt leo (277) 

Formerly, Weisskopf*” showed that the divergence of the electromagnetic self- 
energy of the charged spinor puticle could be reduced from the quadratic to the 
logarithmic by means of the subtraction method based on Pauli’s exclusion 
principle. It is vety interesting that the divergence of self-energies become 
logarithmic not only for the spinor puticle but also for the scalar paiticles by 
means of our subtraction method based on the postulates (1’), (II) and (III’?). 
The first order interactions for ¢ and f, which are left out of consideratioa ia the 
definition of vacuum according to (III’?), contribute only the logarithmic diver- 
gences in (26) and (31). Therefore, even if tlte postulate (III’?) should be 
altered in future, it would be possible to conclude that the quadratically divergent 
terms in the self-energies of above scalar mesons disappear. It is also noteworthy 
that g-interactions does not contiibute to divergent terms in (27) 
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(iii) The self-energy of a charged spinor particle due to the neutral scalar 
field 


The alteration of the equation of motion of ¢ due to the interaction (12) is 
{—-#}p=/o* Be. (35) 


While, the self-energy W, of a charged spinor particle due to the interaction 
G2) is? 


W= sats log 2b + finite term} oh* Bob. (36) 
nr \ p I> 


According to (III’?), (85) does not give rise to the change of zero-point energy 
so that, (36) is again the self-energy of a charged spinor particle due to the 
neutral scalar field in the argument based on (1’), (II’) and (III’?), 

(iv) The electromagnetic self-energy of a charged spinor particle 

The equation of motion of the electromagnetic potential A, interacting with 
the charged spinor field is 


D4, =ep*a,g. 


According to (III’?), the electromagnetic-self-energy W... of a charged spinor 


patticle may be the same one as in usual calculation ;™ i.e. 


Wye ke log rat + finite term} P*By. (38) 
L isa 


The results for the self-energies of a charged spinor 
£ particle remain unchanged in the case of using the 
ep _ postulates (1’), (II’) and (III’?). But these features would, 
of course, depend on the validity of postulate (III’?). 


§5. Transitions of elementary particles through 
(4, Ho) (e.v) the vacuum particles 


aie We have shown in §3 that the postulates (1’) and 


(II’) has influence upon the process of y-decay of neutral meson. This fact 
means more generally that (I’) and (II’) would play an important role for the 
transitions of elementary particles through the intermediation of the vacuum 
puticles. In all problems which connect with the vacuum particles, the definition 
of vacuum always comes into question. We take m-4 decay as the simplest 
example of such problems. L. I. Schiff® investigated the case of Fig. 1 in which 
. %mesons, nucleons, #-mesons and light particles interact each other. This is the 
important problem for determining the model of two mesons, 7 and #. Since the 
m-4t decay or m-e decay based on the model in Fig. 1 is intermediated through 
the vacuum nucleon pair, the divergence occurs and Schiff’s . conclusion in- 
timately depends on the way of cut off, In fact, Steinberger was led to the 
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different conclusion from Schiff by using Pauli’s regularization method. Such a 
divergence depends on the physical character of the vacuum and on the structure 
of elementary particles. We analyse it from the view point of (I’) and (II’). 
Accoding to Schiff, we take m-meson gy (mass x) as scalar, and # meson $y. 
neutral meson ¢, and nucleon ¢ (mass #) as spinor. By means of following 
interactions : 


z-meson-nucleon. fY*rPyy (scalar coupling), (28) 
nucleon-#-neutral meson: g(¢*cP¢) (¢,*8¢,) (scalar coupling), (39) 
the matrix element H,,, of the 7+» decay through nucleon pairs, is 
© Ay ‘ ‘dp 1 
iiss | I Ed i +2 ? tT) OPO, 40 
Ls zw \WJo (22+ pw)? (P+ 8)” (4p? +42—22) ( )¢ PO. ( ) 


by. using ordinary perturbation calculation. i.e. it has the quadratic and logarithmic 
divergences. On the other hand, the equation of motion for nucleon is 


{(4, p) + But fe) +7sLl}p=0, (41) 
L=(6,*Bb,). (41) 
Accordingly, the energy of vacuum nucleon is, according to (1’) and (II’) 


briex Se 
(27)° 


|\iov Pe GA Ter ely (42) 


Expanding the above expression into power series of g and f, the fg-term is as 
follows : 


_ te _ Pdp _(.2)4* (43) 
fale). Bhp 
Therefore, according to (I’) and (II’), by subtracting (43) from (40), the matrix 
element A’,,, of 7p decay becomes 


: ee fee pap . ] cr) d,* : A4 
sare tetra arene} CATA (44) 


i.e. which is of the logarithmic divergence, and 


H’g2y(*%=0) =0. (44’) 


§ 6. C-meson theory 


There is Sakata’s®”® C-meson theory as the substantialistic method ex- 
tricating from the diverging self-energy of charged particles. Since, according 
to above argument, the self-energies of elementary particles are, however, altered 
by the definition of vacuum, C-meson theory must also be modified slightly. In 
this case, it is sufficient to regard the neutral scalar field in §4 as C-meson. 
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Since the self-energy of an electron due to the C-meson and the electromagnetic 


are not altered by using our postulates, the self-energy 


sing our postulates, the self-energy of an 
1)3) 


self-energy of an electron 
of an electron are not altered by t 
electron becomes finite by setting a well known relation 


fare (45) 


from (25) and (38). 
‘The condition for cancelling the diverging part of the electromagnetic self- 


energy of charged scalar meson by means of that due to C-meson becomes very 
simple, because their quadratic divergence in usual calculation was reduced to 
the logarithmic one by using (I’) and (II), that is, from (29) and (34), the 
condition for obtaining the finite self-energy of charged scalar particle is 


f= Se, (46) 


g° in interaction (19) may be taken to be arbitrary, and so the most simple case 
is =U, 

The problem of the self-energy of C-meson is of less difficulty than usual. 
Its self-energy (24) due to the charged scalar meson is finite, and its self-energy 
due to charged spinor particle is of the logarithmic divergence. (cf. Eq. (18)). 
But the self-energy of C-meson cannot be finite by means of the mixture of 
only the charged scalar mesons and the charged spinor particles. However, (1) 
the vector mesons are now left out of consideration. (2) The nature of C-meson 
is ambiguous in many points and it is obscure with how many particles the 
C-mecon interacts. From these reasons we cannot unambiguously conclude about 
the self-energy of C-meson. These divergencies which still remain after using 
postulate (I’) and (II’) seems to connect with the problem of structure of 
elementary particles, because they depend on the mass of real field (not on that 
of vacuum field). 


§ 7. Concluding remarks and summary 


Determining the energy of vacuum by using the postulates (I’) and (II’) 
(and (IIIT’ ?)), usual non-gauge-invariant results appearing in various problems 
can be made gauge-invariant. The internal conflict, that the non-gauge-invariant 
results have ever been obtained in the theory in which the gauge-invariance had 
been verified formally, seems to consist in the incompleteness of usual definition 
of vacuum. The postulates (I’), (II’) and (ITI’?) have an impoitant effect not 
oa upon the problem of gauge-invariancy, but also upon all processes connected 
aye particles. By means of this method, divergencies of self-energies 
; ry particle can be reduced to lower order. Especially, in the same 

manner as Weisskopf’s subtraction of vacuum, by which the order of divergence 
of self-energy of charged spinor particle could be reduced to the bearithaas 
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order, the divergence of self-energy of scalar mesons is also reduced to the 
logarithmic order by means of the subtraction method based on the postulates 
Gi’); and (Il). 

Moreover, these postulates reduce the order of divergence for the various 
divergencies appearing in usual calculation of transition between mesons through 
vacuum particle pairs. 

Hitherto, it is known that the direct interaction between spinor fields make 
divergencies of higher ozder. This problem must be investigated from the point 
of view based on above arguments of the vacuum. 

The definition of vacuum based on (I’) and (II’) is, however, of only the 
correspondence theoretical character and it is the purpose of this paper to point 
out that to fomulate such a correspendence theoretical argument into complete 
quantum field theoretical formalism will be one of the important works to be 
accomplished. 

The authors wish to thank Professor S. Sakata for advice and encouragement 
throughout this work. They are greatly indebted to Dr. M. Taketani for his 
stimulating discussion. They also wish to thank Dr. T. Inoue and Mr. S. Machida 
for their helpful discussion. 
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of the vacuum but also the change of the transition matrix, because the vacuum-state from 
which the transition begins is different from free negative-energy-state. As the transition 
through the vacuum is at least second order process (/*) in perturbation calculation, the effect 
of above alteration of the transition matrix appears firstly in the third order process. So, for 
the problems of self-energies and of z-y decay in the second order approximation in this paper, 
it is sufficient to consider only the alteration of the energy of the vacuum. But for the third 
order processes such as the neutral meson decay it may be neccesary to consider not only the 
alteration of the energy of the vacuum but also the alteration of the transition matrix. 
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Introduction 


We need the detailed knowledges on evaporation process for the analysis of 
various expeiiments on nuclear stars, which have been rapidly increasing in last 
several years. Previously, we studied general characters of nuclear stars, but the 
evaporation process was discussed too little there”. 

Main defects of our previous treatment are as follows: i) neglect of «-particle 
emission, ii) rough estimation of fluctuation. We improve these points and extend 
our Calculation in this paper. 

The problem to be considered here is ‘‘ How many and what kind of patticles 
are evaporated from the certain mother nucleus with the given excitation energy?” 
Applications of our results to various high energy nuclear phenomena will be 


done in a separate paper. 


1. Preliminary remarks 


Adopting the degenerate Fermi gas model for a nucleus”, we have the 
following relations between the entropy S, the excitation energy , and the 
temperature 7 for the nucleus (WX, Z), ie., the nucleus of MW neutrons and Z 
protons (V+ Z=A=mass number). 

1(/N-ZY) 7 2 
X~{1—— pier Hs AT? 
9\N+Z/) 4H, 
1/N—Z\)\( AX \3 
seam (S52) (AE 
e8) 18\ V+Z 


E,=Fermi energy. 


In these expressions, each second term in braces is the small correction depending 
on V—Z. We may, therefore, safely neglect these correction terms. (However 
we can not neglect the contribution from the term depending on V—Z for the 
binding energy of a nucleon; see below and ref. 1). Inserting Serber’s value 


22 MeV for Fermi energy E,”, they become 


788 Y. Fuymoro and Y. YAMAGUCHI 


(1.1) 
S,(X) ~0.63, (AX )?, 
X, T measured in MeV. 


Here we do not take into account so-called Bardeen’s remark®, for (1.1) gives 
the energy distribution of evaporated protons which is in hettes agree with 
the experiment of Perkins ct a. In other words, we are inclined to regard 
(1.1), so to speak, as the empirical formula. 

Bohr and Wheelei® gave the following expression for the main part of the 
total binding energy of a nucleus (V, Z)*; 


B(N, Z) =AL,— > b,(Z-Z,)", (1.2) 


where H, means the average binding energy of a nucleon for the most stable 
nucleus (V4, Z4) with mass number 4. &, and 4, are the slowly varying quanti- 
ties of 4, and 4, was tabulated in Bohr-Wheeler’s original papei®. 

_ Then the binding energy 4, and 4, of a proton and a neutron for the 
nucleus (4, Z) are given by 


= wht — Bal eet 
CA N=const « 
(1.3) 
apes =F,—b, (W—N,), 
OV Z=const 


respectively. The differentials of the slowly varying quantities are neglected. 


Next, we will consider the nucleus (V—x, Z—p) (a=u+/, assuming a<A), 


which is produced from the “ stable” mother nucleus (VV, Z) by carrying away 


z neutrons and # protons. The most probable charge Z,., for a nucleus Of 
(A—a) nucleons is approximately . 


PPE TF PRS WEN FON Eas 4 
dA 2 
where 
dN 4-24) dM A=2Z) _ 
dA aA as 


It is easily calculated from the Bethe-Weizsacker’s semi-empirical formula for the 


binding energies”. Substituting this expression into (1.3), we obtain the follow- 
ing formula for Z, and &, of the nucleus (V—x, Z—p): 


* (1.2) is derived from eq. (6) of Bohr-Wheeler’s paper (ref. 6). 
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(1.4) 


These results were used in I. The line 


1+ 
1 cP 
at 


1epresents the ‘‘ Heisenberg valley” near the point (4, Z) in the nuclear chat. 
The binding energy of an «-particle in a nucleus is given by . 


xu= 


E,~2E, +2E,— (total binding energy of an u-particle) 
w4F,—28.2 MeV. 
For heavy nuclei (4100), this is 
E.~4 MeV. 


We may neglect the change of the Coulomb barrier height V throughout the 
evaporation process. It is justified by experiments on energy distribution of 
emitted protons’, and by estimation of thermal expansion of a nucleus”. 


2. The evaporation theory 


According to Weisskopf™, the piobability per unit time that the nucleus 
UA 2X) 1.635 “oncleus a, Z) with excitation energy X, evaporates a “7’- 
paiticle with the energy ¢, +d? is given as follows 


WX, ede) Ms PG —Kj-B,) exp [Sa-a(X—8) SX) 48. (2.1) 


4“ 


We write mass, charge and statistical weight (due to its spin) of this ‘7'”-particle 
as M,;=ayn, 2¢ and g;, respectively ; where # means nucleon mass. Furthermore, 


we use the following notations. 
R=radius of mother nucleus 
=y7,A"8, 4,=1.37 x 10-% cm” 
V,=Coulomb barrier height for j-patticle, 
£,=binding energy of j-particle, 
S4(X) or Sa-a, (X) means the entropy of the nucleus with excitation energy X 
before or after ‘the evaporation. The kinetic energy of an emitted j-particle is 


c—£,. 
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For X¥>¢, we may use the following approximation, 


ay, a 47k, 22 
Se) S4(X) G7 = (2.2) 


T=(dS,(X)/dX)~* means the temperature (measured in MeV) of this excited 
nucleus with excitation energy X. Then the energy distribution of an evaporated 
j-particle is 
i 
(€¢—£,—V,) exp (“1 *), 
a ' 
The mean kinetic energy of an evaporated j-particle is V,+27, and the 
mean decrease in excitation energy of the mother nucleus due to emission of 
j-particle is 


V;+E,+2T. (2.3) 
The “j-particle width” [,(X) is obtained by the integration, 


x 
Boo =], We 2)aé 
! 


Let a,T V;+E£; 


= 


(2.4) 


where K = #/2mrZ2 = 11.1 MeV. For instance, proton width 1,(.X), (g, =2;, 
M,=m, V,=V) neutron width [,(X) (g.=2, m= m, V,=0) and a-particie 
width [,.(.X)(ge=1, me=4m, Ve=2V) are as follows: 


1 Aez? T EMV 

io eee ey ee ae 

hdl ag exp ( 0 T ). (2.59) 
ices kis ek teak a le 
D(X) 0 — pe (= (2.52) 

and 

l ABT? 47 E 2 
ext vae’ (47 _ Bet 2V 

eer exp ( 10 cae, om 


For later use, we write the ratio [¥(X)/U,(X) as, 


P(X) E,+V—E, T 
<- sw kex ee eet 6. 
7G) 7 oot ecto Rene 


The approximate expression of the fission width is®™ 
“ 1 &, 
yr) shes (- f 
W(X) + Pexp (—2), 25f) 


where £, means the fission threshold, and we assume that the small deformation 
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of the mother nucleus causes the fission process (cf. ref. 6). 

Formerly Tomonaga™* calculated the coefficient of heat conduction ana of 
viscosity of nuclear matter and showed that they were very large. Comparison’ 
of the heat conduction with the particle evaporation leads us to the next two 
conclusions : 

(i) It is very rare that a pait of a nucleus is locally superheated and 
evaporates any paiticle. 

(ii) “After evaporating a particle, the residual nucleus will reach the thermal 
equilibrium state before next evaporation takes place. 

The former (i) was previously used to conclude the possibility of fission 
process in very large cosmic ray stars. We see from (ii) that a particle is’ 
always emitted from a nucleus in thermal equilibrium. 


3. The diffusion equation for evaporation process 


Let us suppose the mother nucleus (, 2; X,) at the biginning (time 
t=0). Various paiticles, will be evaporated from this nucleus. The probability 
F(N, Z, X, t)dX of finding a nucleus (4, Z) with excitation energy between X, 
X+dX at time ¢ satisfies the following simultaneous diffusion-type equations : 


BF, Z, X, 1) _ i" 


F(N+a;—-2 Z+2; X', t) 
ot j 


X+V 5+ F; 


x WX', X'—X)dX! 


SITE) Ze Xv) (" W,(X, €) a. (3.1) 
j V 5+ ¥; 


The initial condition is, of course, as follows,, 
FN LeN bh, dead Haak) (3.2) 


‘Solving (3.1) under the initial condition (3.2), we get the desired answer to 
our problem of nuclear evaporation : i.e., 


min (V 5+ ¥;) 


FN, 2) = | P(N, Z, X, +0) aX. (3.3) 
This is the yield probability of a residual nucleus (JV, Z) from the mother nucleus 
(Ny Z; Xo). Other required quantities are derived from F(J, Z, EN, ft): & 
It may, however, be necessaty to perform very long numerical rae 
to solve our problem in this way, although this was done by Horning ay 
Baumhoff for some simplified case”. (Also interesting formal considerations © 
eq. (3.1) were done by Baumhoff.) Avoiding this direct attack to our ey 
we prefer to examine in the first place the average behaviour of evaporating 


nuclei, and then take into account the effect of fluctuations. 
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4. Mean behaviour of evaporation process 


Now we discuss the mean behaviour of the evaparating nuclei. Mean energy 


E® carried out of the nucleus: (V, 7; X) by emission of j-particle is given by 

E®=V,+ £,+27, 

T=3.17(X/A)}. (4.1) 

The nucleus (WV, Z) is the intermediate one in the evaporation process, after 
emission of 7 neutrons and / protons, i.e. 

N=N,—1, Z=Z,—}p; 

(N+Z=A, N+2,=Ay n+p=a) 

and we. assume Spee 


When the nucleus evaporates, on the average, one nucleon, the mean decrease 
of its excitation energy is equal to 


While the probability or the width for emitting j-patticle is 


> - £545 ATF a V,tE;> 
(MN, 2; X) = exp(—7g7- “iy 


(4.3) 


This set of equations (4.2) and (4.3) determines the mean number of 7-particles 
y emitted before the nucleus is cooled down to the excitation energy X. 


= _ Xo T,daX 
v,(X, Xx ) =). SPE” (4.4) 
k 


Thus the average number (¥,(.X,)) of 7-particies evaporated throughout the whole 
process is equal to »,(.X, V;+£;). Of course, number of neutrons and protons 
mz and p carried away from the mother nucleus are given by 


p=da 55% 5) n==>) (@,—2,)» . (y=n, Py Gyrreree ). 


The curve wy=Y,(X', X ) (Sey Py dyrrrres ) in »,-space represents the mean 
bzhaviour of the evaporation process and we call this “ evaporation path’’. This 
is a generalization of Bagge’s mean cooling law. 

For practical purpose, it is sufficient only to take into account neutrons, 
potons and «-particles, since other particles have much smaller widths. 

(A) Case (zp) Consideration of only neutron and proton emission (case 
(7~)) was discussed briefly in I. Let us examine this simplest case more closely. 


In this case, we must solve th 3 
’ e follow : i 
v,=2) ing set of equations (of course »,= / and 
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aX _ (En +27) +,(V 4H, +27) 


O(n+p) Latls 
ee . (4.5) 
n 
= me sa =p=0 r N= NX). 
IB u=p for X= XX, 


Explicit expression for each quantity is as follows. 


eke ), 
1—p 


pe as(= =) ; (4.6) 


For 4,=100, they are 


E,~8 Mev, V~6 Mev, 


4.7 
p~0.9, ote NL oY 


eke 
Various properties of the solution of these simultaneous equations can be obtained 
from consideration on following simplified models. 
(i) Evaporation process with constant temperature. 
This neglection of cooling allows us to treat the only one equation 


_7tt 
oh ae i ae =?) (4.8) 


Its solution (under the condition »=0 for p=0) is, 


23n 


1+yp er? —l v 


lpg on! oF Fam gal (4.9) 
For large 2 and #, this tends 
26(n—= disk tt p)-V=0, (4.10) 
—ye 


which is parallel to the “ Heisenberg valley”. '(See Fig. 1.) 
(ii) We assume that the ratio T/T Ep and &, are constant through 
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evaporation process.* In this case it 
needs only to solve the equation 


“" = N,4+27, 
da 
re 
X,= £,+—2—V, n+ p=a, 
we x\ 
E,=E.=Ey r=3.17(*_/. 


The initial condition is a2=O for X 


=X,. The solution can easily be 


found as 
n xvas 
Ag—a a 
OE “ea aneers 


a 


| vf —atl, 


Y¥ —a+7, Vv +a-—T, 


Heisenberg ce (4.12) 
NUN, oe 9 where 
Y =v 
Asymptotic 10 a(u+X) 
/ ee X= Ae—4 7+ ar510 


—— curve (4.9) with (4.7) (for 7=3 Mev) 
wees solution of (4.5), (4.6) and (4.7) T, denotes the initial temperature of 
(4o=100, Xo=100 Mev) the nucleus and 7, the temperature 
Fig. 1. of the nucleus just after evaporating 
/ a nucleons. This solution (4.12) 
was also obtained by Perkins ct al.” It is woith while to note that a/A, depends 
not on A, but only on 7, and 7). We think that the evaporation process Stops 
when the excitation energy decreases down to X= X,=£,4+(5V/1.4+6). 

(See Fig. 2.) 

We see from the above models that the evaporation path determined by the 
exact equation (4.5) has its end point close the asymptotic line (4.10) and its 
path length is nearly equal to that of (4.12). The results of numerical integra- 
tion of (4.5) are given in Table 1 and in’ Fig... 3, a 


ee 


* Y. Fujimoto and Y. Yamaguchi, Lecture at t 


Dec. 12, 1948), he Cosmic-Ray Symposium (Nagoya University 5 


High Energy Nuclear Evaporation 795 


9 

8 
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MeV. Mev. 
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3 3 

2 2 

1 1 

: "ite 0.1 0.2 0.3 % 0.1 0.2 0.3 0.4 
al Ao al Ao 
xs=8 MeV. xs=10 MeV. 


Fig. 2. Curves of (4.12) 


(B) Case (pa). Next we take into account the a-patticle emission (case 
(upa)).* Then the cooling law becomes, 


AX _Tn(Eqt2T) +0 p(Ept V42T) +1 e(Zat2V +27) 


Ba T, $1, 44K, 
aie Ung Ot 6 (4.13) 
dp T,+2l, 


a=n+hP, N=VY, +2v,, pHYMy+2Nq. 


If we adopt the expression (2.ba) for a-particle width, the number of evaporated 
a-particles becomes much smaller than that observed in cosmic-ray stars (Table 2). 
This is readily seen as follows: firstly, 


and 


bin iQ He TTT (4.14) 


ap Tp 
was used and the results stated there are considerably different from the ones of this paper, 
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Table I. Numbers of evaporated 
particles ; case (724) 
(obtained from (4.5)) 


Ag=100 
(i) E,=8 Mev, 1’ =6 Mev 
p=0.9, Ata iy7 
_—- 
Xo } < f - < 7 > 
50 Mev| 0.27 | 72 
100 | 1.19 4.84 
150 | 2.20 6.54 
200 3.15 8.15 
400 | 6.4 13.6 
a0 | 121 3 
Ay= 100 


.. )£;=8 Mev, 7=6 Mev 
(ii) 
pe=00, It#217 
l-yz 


yh 800 Xo | < ? a < ” > 
ae 15 100 Mey | 0.48 6.52 
pail ) | 
(n= 98.2) 
Agp= 200 
Fig. 3. Solution of (4.5), (4.6), (4.7) (iii) E,=8 Mev, 7 =9 Mev 
p=0.5, i+ #20 
This expression has its maximum at ¢ 
- (Gatv-B) - Xo | <s> | <#> 
~ . ee 
g ae Vn Be 200 Mev | 153 | 10.82 


Thus, for the numerical values Table II. Numbers of evaporated 


Ey~4 Mev, E,~8 Mev, l’~6 bat 414 ( paleg tin of (4.13) | 
‘ ; case (fa), (solutions . 
(A~100), | ( ) (3.5 a), (4.14)) 
this is epreaticese scenic 
A X |< >|<vp>|<re> 
(22)  ~26*=0.27 es boric a 
| iE oes 200 Mev 7.9 2.6 0.67 
\ . : 400 18.5 5.0 13 
While, the observed ratio v,/¥, in cosmic ray 800 22.0 9.4 1.9 


stars is somewhat larger, i.e.~0.5.™ 


This discrepancy may he possibly due to wrong A-dependence of assumed 
entropy of a nucleus. It is our regret that we can not estimate the above possible 


effect. Tentatively, we assume three constant values, 1/2, 1, and 2 for & (4.14) 


and sol ati i ic i i 
olve the equations numerically. »,/v, ratio observed in cosmic ray stars 
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Table IMI. Numbers of evaporated particles, case (7/a). 


Ay=100 ens Ay 
£,=8 Mev ine Se = exp (225e=*) 
Ey=4 Mev Pat.? i d 
V=6 Mev 
k=2 k=l | k=1/2 
X 
| Rva> <yp> <ya> || <ie><vprivern | <ra> <vp>’<ve> 
100 Mev 4.7 0.8,  O.4 47 10, ‘02 J 
200 77 (ge) 1., 78 23, 0.9 78, 26, 0.4 
400 1g0 Shalt, Seas 
800 oan eae 20.8 6.7%, 5.48 a ae 


seems to be most favorable to the 
value~1 for %. (See Table 3.) 

(C) It is woith while to point 
out that the value (¥,(.X%)) + (¥.(X)), 
i.e, average number (P(X,)) of 
prongs in one star, for given 
(and also A,) is not sensitive to 
assumed values of 4 (Table 3). 
Thus we may hope that the average 
number of prongs (for given 4) 
will not much change even if & is 
a rather complicated function of 7. 

Therefore it may be allowed 
to divide the value (P(X,)) thus 
obtained into two parts, (»,(X)) 
and (y,(X)\, such that their ratio 
{y,)/{¥,) is in agreement with the 
experimental one. We plot in 
Fig. 4 the values (»,) and (»,) 
versus the initial excitation energy 
X, (for Aj=100) determined in 
this way. 


5. Fluctuation 


As Weisskopf’s distribution has. 


a broad maximum, actual kinetic 
energies of evaporated particles are 
not always equal to their mean 


values, so that the number of eva- 


800 600 


(Ya) GERM Jani 


400 300 200 100 
a, MeV. 
Fig. 4 
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porated particles show some probability distributions. We discuss the effect of 


such fluctuations 
(A) Fluctuation of @ (number of nucleons carried away fiom the mother 


nucleus). We may regard the evaporation of nucleons ane of a-particles ae 
independent to each other, for the evaporation of an a-paiticle does not result in 
change of binding energies. Therefore we treat them separately. 

We discussed in I the fluctuation of number of evaporated nucleons Y=», +¥,. 
The probability #(’) for evaporating just a nucleon is well approximated by 
Gaussian function as follows: 


(= y= exp [-“-e (5.1) 


V/ 7 dv 2(dv/)2 
where 
At=V2T 
£, 


The probability g(¥,) for evaporating », a-particles will be Poisson’s type, 
fA < Ss (¥,)~%a 
sv Jer. ars 


since the evaporation of a-particles will have little correlations. 

We get the distribution of total number »=»! +», of evaporated particles by 
combining these two distributions. Compaiing 4 with 4y,= (¥.)’, we see that 
4v> Av, in most cases. Thus this distribution is mainly governed by (5.4). 

On the contrary, the distribution of total number @a=v +4v, of nucleons 
carried away is mainly determined by that of «-particles, since 4 <44v, in most 
cases. For some special cases, such as low energy evaporation process, 44», 
becomes smaller than 4)’, and the distribution of »/ plays an important role. 


(B) Fluctuation on x or ~ (Case (7p)). Formerly we assumed binary 
distribution for the number / or » of evaporated protons or neutrons, i.e. 


(m+p)! __ (n)"(p)? 

m\ pl ((n) + (py)? 

’ This neglection of correlation in successive evaporation of a nucleon is only valid 
for the case of small {~) and (x ss 


When (#) and (%) becomes considerably larger, i.e. larger than the 


breadth of Heisenberg valley, we can not neglect the effect of change of binding 
energies, which imposes serious correlation between them. 
treat their fluctuation as Markoff’s chain. 


For simplicity we put (i+e)/(— z)=1, then the ratio of the widths is 
Fon exp {P=2BO=A) | 


(5.2) 


In this case, we may 


(5.3) 


Pp 


Assuming the constant temperature 7, we get the diffusion equation for the 


—— 
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probability F(x, ~) for finding the nucleus (V—7, 7—p). We get its steady 
solution /(7, ~), which depends only on 2—A, as follows : 


23(n—p) —¥ 4 4 
EG@.f\=Glichem Medieasl civ (5.4) 
where C is determined by the normalization condition, 
+0 
3 Cals dee 
n—p=-© 


Of course, we should use the average temperature 7 suitably chosen. 


6. Summary and acknowledgement 


We calculated the prong number distribution and its relation to excitation 
energy on the basis of Weisskopf’s evaporation theory. These results are applicable 
to high energy nuclear reactions, meson stars, cosmic ray stars, etc. 

The authors express their hearty thanks to members of Tokyo University 
and St. Paul University for their kind discussions on this work. 


It seems to be necessary some additional remark to avoid misunderstanding. We can not give the 
complete explanation of nuclear stars only using the evaporation theory considered here. But we consider 
that particles emitting in the star are roughly classified in two classes: the one is so-called “ Tleisenberg 
nucleons” and the other is evaporated particles. This classification is justified by the recent experiments.* 
Thus we considered the id2a/ evaporation processes somewhat in details in this paper. To make a 
discussion of nuclear stars more complete, we must consider the Tleisenberg nucleons etc. 


* Cf. ref. 5), also see C. F. Powell et al, Phil. Mag. 40 (1949), 862. 
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§ 1. Introduction 


In connection with recent experiments of aitificial meson production by 
nucleon-nucleon collisions, several authors” have calculated the cross-sections 
for these phenomena. However, their results are more or less unreliable because 
of the arbitrary choice of model for z-meson on the one hand and of the doubtful 
validity of the perturbation method on the other hand. As for the first point, 
we shall discuss all the four possible models for z-meson to eliminate the arbitra- 
riness. But, as for the second point, it is admitted to neglect damping and 
reactive corrections, because near the threshold energy for z-meson production 
mesons are emitted from loosely bounded parts of their clouds about the nucleons. 
In §2 we shall give a general formula for collision cross-sections in a completely 
relativistic form, and in §3 explanations for the approximation of our calculation 
will be given. From §4 to §7 cross-sections for all the four models are calculated 
and in the last section (§8) we shall compare them with each other. 


§2. General formula for cross-sections 


Since the calculations are much simplified in a completely relativistic form,* 
we combine Méller’s formula for collision cross-sections and Feynman-Dyson’s 
diagram method for calculating S-matrices. We consider the most general collision 
process as follows: two particles collide with each other and then # particles 


including the original ones are produced. From M@ller’s formula the cross- 
section is 


= 7 ee § KK) 3( WW) | (Kegs, Kel Z [Ay Hy)? (2.1) 


t=1 J 2; 
with 


S= Vv | kw, —kew,?—|k x ke. (2.2) 


0 0 0 
(k,", w,°) and (k,°, w,") are energy momentum four vectors of incident particles, 


* ‘ Se508 
C. Morette has already computed jin a relativistically invariant manner the cross-section fer 


meson production by nucleon-nucleon collisions and we are much indebted to her manuscript. 
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respectively, and (#*‘, w*) is the same for the z-th particle produced. (K°, W) 
and (K, W) are total energy momentum four vectors before and after the 
collision, respectively. 6 results from the relativistic generalization of relative 
velocity of incident particles. The relation between a transition matrix (K,,--, 
K,,| I | h,°, k°) and a corresponding S-matrix in momentum space is 


(K,, os k,,| Si k,, k,°) =0(K—K)d(Ww— W?) (kK, “ot, Kn | T| f a Ke") 


0.) 0 
V Wy *Wy =A Wy We 


(2.3) 


From (2.3), we can readily find a way to calculate /-matrices relativisti- 
cally with slight modifications of Feynman-Dyson’s diagram method for S-matrices. 
As an illustration we assume particles to be either nucleons or mesons. 

The wave function of a nucleon and the adjoint operator for it are as follows: 


ax) = VT/ZQAY D [dv / Vn ep (GE (D1 «Dyes a+ Gz (DY o( Det), 
(2.4) 

$.(#) = V1/2Q2)= vo {ap/ Vn + pp Pt (Diy, (pera +, (DP) ur, a( pera). 
ae (2.5) 


Here m is the nucleon mass. 7+(=1,2) and +, — distinguish the direction of 
spin and the sign of energy, respectively. (We use the natural unit system 
4=c=1 and the metric tensor g¢4=£»=£3=£u= —1 throughout.) Commutation 
relations for these decomposed fields are 


{bt (p), Gh (D’) ba={ 95 (DP), Ge (B") bs =O er Vin + OD — DP’), (2.6) 


remaining anti-commutators=0. 
Spinor wave functions satisfy the following relations : 


>a Ut g( Pp) Us, 6(P)= (Tr Put Mt) a3 and pa (DP); s(P)= (7. Pp —) os (2-7) 


with 7,=f8 and y=fa. 


Similarly for symmetrical meson wave functions, 
$(2) = V12Cay [dhe / V p84 BE (Wea HWM) E=1,2,3), 28) 


(for mesons with spin 0) 


or 
: ft — x yee 
G(x) = V1/2(22)? fan/ V4 Ke vi, (Ke) (Get (Keenta +4; (Kee M%a), (= 1, 2, 3) 
: (2.9) 
(for mesons with spin 1) . 
where 7,(£) is a polarization vector of a meson with spin 1 and satisfies 
3 7 . 
oD v7,(k) Un (k) =—Ouy + CAPT) ‘ (2.10) 


r=1 
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commutation relations are 
[si(he), * (el) a0" Ve he O(k—K’), 


remaining commutators=0 


(2.11) 


for mesons with spin 0, and similar expressions for mesons with spin 1. 
Then for each figure of Feynman-Dyson’s diagram method,* there corresponds 
one transition matrix (|/|) which is obtained by the following procedure. 
(I) For each free nucleon line ~“1/2(27)* «*(p) or expressions like that 
is substituted. 
(II) For each free meson line “1/2(2z)* or “1/2(2z)* u%(F) is substituted 
according as the spin of the free meson is O or 1. 
(III) About other elements of figures as internal lines and world points there 
is no difference between S- and /-matrices. 
Thus we have obtained a completely relativistic procedure for calculating cross- 
sections and this will be used in the following sections. 


§ 3. Approximation of calculation 


Ae shall calculate the transition probability from a state in which we have 
two nucleons with momentum p® and gq to a state in which we have two nucleons 
with momentum p and q and a meson with momentum &.** Hereafter, we assume 
that quantities with a symbol ~ refer to the laboratory system and those with no 
symbol refer to the center of gravity system. The threshold energy &, for z- 
meson production and the maximum energy €,, of emitted mesons when the 
incident energy £& is given are as follows: 


€ m= L/2 + 2/4in + a E/(2m+E) A (E/2—2/4in)— pe (for the laboratory system) 
and 


F =2u~280 Mev. 


F =2u4 °/2m~290 Mev. 


€,,= (ml + 4°/2)/,/2in(2m-+B) (for the center of gravity system). 


Z x 2 a ; j 
As €<m near the threshold energy, We may neglect £ compared with g’+p 
°F come: fi One 
so that p°=p"—me.*** Tn what follows we calculate cross-sections in powers of 
e/m and leave only the first non-vanishing term. 


spt St Fuither we shall carry on our 
calculation in the center of gravity system, bec 


ause it helps us very much. 


* A detailed account of Fe 


in the following paper: 
ere we assume the nucl 
these points we 


ynman-Dyson’s, diagram method for calculating J-matrices is described 
F, J. Dyson, Phys. Rev. 75, (1949) 1736. 


eons as free and neglect’ effects of the target nucleus upon them.. On 
shall have some discussions in §8s. 


See Foldy and Marshak, reference 1, 


* 


ORK 
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§4. Symmetrical pseudoscalar meson 


(A) The interaction Hamiltonian density for symmetrical pseudoscalar 
mesons is 
F(x) =fh(a)rste (x) b'(4), (4.1) 
where 7;=717e737 f is the coupling constant* and 7,’s are charge operators. The 
Feynman-Dyson’s diagrams for meson production now considered are given by 
Fig. (A), (B), (C), (D) and those which can be obtained from them by inter- 
changing # with g. 


qa gq Pp q p o 


fy 


> Ie 


k 


The transition matrices which correspond to these figures are summarized 
into the following expression : 


(p, a, &|Z|p", a) =[1/2(22)*P? (27) 79 (ut (p) reee* (P") )[1/ (Ap) I) 
(%* (Q) {7575 (Qavat Axvr—™) Vet t ste (Qn Ta— Pata — ™) Ts u*(Q°)) 
— (£9) —(P°9°) + (0a, 2°9°) (4.2) 


where (pg) or (f°g°) are expressions which can be obtained from the previous 
one by interchanging ~, with g, or p," with gx, respectively, and (pg, "g') by 
interchanging ~, with g, and ,° with 9," simultaneously. 

Taking the absolute square of (p,q, & [| p’, @), summing over the final 
spin states of the nucleons and averaging over the initial spin states, we find in 


the first approximation 


pane S| 7 ~ me ; 
pakitie a=? (1) ate bee} 8 
(yp, a kT p', a) "=> (57). (p—p) 2 (P—P)'(P" +p) 


where 7 and 7’ are numbers which are listed in Table 1. Mesons are emitted 
spherical-symmetrically and their energy distribution is given by the following 


differential cross-section 


* Cross-sections for symnietric pseudoscalar mesons with pseudovector coupling are by a factor of 
order €/ smaller than those for symmetric pseudoscalar mesons with pseudoscalar coupling, so 


we only consider the latter case. By the same reason we shall neglect interaction I{amiltonian 


densities which contain derivatives of meson field’s operator tor other types of mesons. 
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da=(i/m’) (f?/4%)* P(e) — (7 /2)8 (©) Jee, (4.4) 

with 

AV 2 yp wee 
EP arn m 
f(¢) TRL aah (4.5) 
§& (€) = (€/2€,,—€)f(€). 
: a 37. 
In integrating the differ- Table I. Values of a, B, 7 and 7” 
ential cross-section (4.4) analyt- ans ae RE 
ically, we make the following ap- = a 
proximations: we replace tanh™* its ¥ fF y , © Yee 
——<——— ee ! », * £ a aa) - - 
V (En— €) /Em by v Ca eh ea hs 5. ph A x— 1 1 1 1 
and W(e€ +p)/€, by a suitable | am, ae ; 8 2 ee 
average value (V(€+))/€n) ay - Be ea 1, ae 
Then the total cross-section nm, n | aw = 2 —%/ 3 V3 
is a, nm, wo 0 0 1 2 
a, 2 2 3 (/ and » indicate a proton and a neutron, respectively, 
ane Sli ) i) and z+, z— and 7° indicate a positive, negative and 
x [yr (7//2) G], (4.6) neutral meson, respectively.) 
where 


a wl ‘é —- fE (<a) 
sa 2 (y/ Em pas 


MLE», 


ee We. Cid my (47) 
dE ~(2€,.— 
G =8r(/ ae (=F € Gs ft) )-F. 


Qm a 


(B) Foldy nd Marshak have treated meson production as a second order 
process in which one step consists of the creation of a meson by one of the 
nucleons and the other step consists of the mutual scattering v 
nologically determined nuclear potential between them. 
we take a nuclear potential 7 which is given by 


ia the phenome- 
Following their procedure, 


4 
V=SV, 


(4.8) 


~ 
- 


with 

Vi=(er'/4r) § dade’ (x) f(x) exp (—pr)/rpi(x') o(2’), 

Vy= (gy'/4m) § dacdach (x) y7uh(*) exp (—pr)/rP (a rr“), 

Ve= (ge /4n) § dardac' (x) zh (x) exp (—pr)/rp(a')e(2"), 

Vis (ge/4r) | dxdal$ (a) care (x) exp (—4r)/rh(x')eray(x"), (4.11) 
where r=|x”—a’'| and x,=2x,. 


Though we write spin matrices in a relativistically invariant form and assume 


the pseudoscalar coupling between mesons and nucleons instead of the pseudovector 


one, it will bring no difference to the cross-section, because in our approximation 
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we treat nucleons non-relativistically and the equivaleace theorem between the 
pseudoscalar and the psudovector coupling holds.* 


The angular distributian of mesons is spherically symmetric and differential 
_cioss-section becomes 


tom EMG GE) Yo CEC Ze fee 


(4.10) 
with 
KOH YY Ged les (e2— pe) 92(1/e)4. 
7 (<? — pe) — =f Grae: : 
g(é) =2 Se ee = 7 (4.11) 


The total cross-section is 
- LMC CFCC) 
with 


5 Ves 
Fatt [54h _ eee (<. aah ae bacilli 


97 e+p 3€,—f /Qe,—p a( bb Je V 1(2€,,— 
B vel EB) fo ate feet) (4.13) 


If we assume instead of the expression (4.8), that 


V=(f?/4n) § dacdul f(x) 7(%) exp (—pr)/rb (2) ree"), (4-14) 
which can be deduced from the symmetiical pseudoscalar meson theory. Then 
the cross-section is just the same with (4.4) in (A). 


§ 5. Symmetrical pseudovector meson 


The interaction Hamiltonian density is 


F(x) =f P(x) ror ntl (x) 9. (4): (5.1) 
The method of calculation which has just been used for symmetric pseudoscalar 
mesons can be used in the case of symmetric pseudovector mesons, and in the 
first approximation in which the nucleons are treated non-relativistically the 
angular distribution of mesons emitted is almost spherically symmetric but not 
exactly so. Here we write down only the differential cross-section. 


* The equivalence theorem between the pseudoscalar and the pseudovector coupling only holds 
for the case (B) and not for the case (A). 
** This expression differs from the expression (23) of reference (1) in the sign of g-term. 
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son (LY (4414 yo-(S Sr oof 


_ 


with 
im 249) Sf 2 
=2',/ €m : ( ’ 
Ke) es 3 
pe 
(<9) as 
g(é) €,,(2€,,— €) 


where a, y and 7’ are listed in Table 1. 


§ 6. Symmetrical scalar meson 


(A) The interaction Hamiltonian density is 


H(#) =f («) 7:9 (4)9,(4)- 


In the first order of €/m the cross-section vanishes by the interference between 
four matrix-elements which come from Fig. (A), (B), (C) and (D), respectively. 
So we must, proceed to the next approximation, in which relativistic effects is 


expected to be important. 


10-*/Mev 
10-*/Mev 


0 20 40 60 80 100 
€—yu Mev 
Fig. 1. Energy spectrum for pseudoscalar mesons Fig. 2. Energy spectram for pseudoscalar mesons 
(A) (Eq. (4.4)). (B) (Eq. (4.10)). 


f+(g/2): 2+ (or z-) production by nucleon- 


(5.2) 


(5.3) 


(6.1) 


nucleon collisions. F+(a1) iin* Cone) Oe ee te 
nuclecn collisions via Serber potential. 


Production of x-Mesons in Nucleon-Nucleon Collisions 807 


The differential cross-section is 


da= (1/m*) ( f/4n)'[af(e) —g(e) Jae (6.2) 


with 
He) = SMO. 8) = Nee). (6.3) 


That is, f(€) and g(€) are composed of three parts, each of which has a 
different origin. 

(I) As a momentum transfer between the nucleons is not exactly the same 
in the process (A) (or (B)) as in the process (C) (or (D)), the remaining 
fraction which evades cancellation gives cross-section of the order of e/m. f,(€) 
and g,(€) come from this effect and are given by 


aa) A Piet m(2€,,—€) (€2— p2) 9? 
(3 m m A 
Ale) 3 Ene €6 


(6.4) 
E1(€) = (A (€)/4) (€/2€,,—€)*§(1+ K(€)), 


where 


USGS. 


tanh (= cw 
(2e,, = €\) eee im 


1/Mev 


(P.V)(V) 


vector 
0.8 1.6 


0.6 1,2 


GORA tT hie ETE e 


pseudovector 


vector 
Ral ~~ — —— 
ee gE OPE 9 2 +40 60 80 100 
¢—yp Mev ¢—p Mev 
Fig. 3. Energy spectrum for scalar meson (A) Fig. 4. Energy spectra for vector (Eq. (7.2)) 
(Eq. (6.2)) and pseudovector mesons (Eq. (5.2)). 
nd : (ee f: x* (or x~) production for the vector theory. 
Lo Eis een Sea brerclon bya Lon L. (0 (Because g can be neglected compared with /) 
N—N) collision. (21/8) f—(13/32)g: 2+ (or x-) production for 


f—&+ x* (or 27) production by P—WV collision. the 7.v. theory. 
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me a 


x10-” in ¢m*/Mev. ster. a= (a/m*) 


32 
5 nucleus, 


0 40 =g0 120 160 200 


€-p Mev €..—f Mev 
Fig. 5. The experimental energy spectrum Fig. 6. The total cross-section for production 
at 90°. of #.s. mesons (B). 


(Here we assume Serber potential.) 


(II) The second part of the cross-section is obtained by considering the 
nucleons recoil effect : 


fr(6)=(A(©)/2)L1— (K(€) /8)], 
&2(€) = (A(€)/4) (€/2€m—€) [1 + (5K(€) /8)]. (6.5) 


(III) The interference effect between (1) and (11) makes the third con- 
tribution : 


fs(e) =A(e)[1 + 1/16) (€/2€,,—€)*(1+K(€) +2Z(€))], 
Sal €) = (7:(€)/2) (€/2€m—€) [1+ (1/4) (€/2€,—€) (14+ K(€) +2Z(€))], (6.6) 


where 


© a 
L(6) = joer tan / Sa € 
(<) V €m(Em— €) be e. @ 
The exchange effects,* however, are still of higher order in €/m by the 


interference between four matrix-elements and may be neglected. The angular 
distribution of mesons differs from spherically symmetric one. 


* Exchange terms are always smaller than ordina 
Scalar or pseudovector theory. So 
theory and calculate the cross-sect 
have the same results, 


ry terms for the symmetric scalar, vector, pseudo- 
if we assume a potential which is deduced from the meson 
ion for meson production as a second order process, we shall 
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x 10? 


RS 1 PSA) 


0241 10° SCA) 


0.1 


~ (21/8) F-(13/32) G 


0 0 
0 40° 80 120 160 200 0 40 80 120 160 200 
€n—P Mev €,,—# Mev 
Fig. 7. The total cross-sections for Fig. 8. The total cross-sections for production 
production of S(A) and P.S. (A) of Vector and P.V. mesons. 
mesons, 


(In Fig. 6, 7 and 8 we plot only the average total cross-sections for the incident proton over 
the nucleons in the target nucleus having the equal numbers of protons and neutrons.) 


The energy distribution curve of mesons and the total cross-section for 7- 
meson production are plotted in Fig. 3 and Fig. 7, respectively. 

(B) Next, as in the pseudoscalar case, we assume the phenomenological 
nuclear potential 7 ((4.8)) and treat meson production as a second order process. 
The differential cross-section is composed of three parts, each of them corresponding 
to effect the (I), (II) or (III) in the case of (A), and is given by 


iene Py th (6.7) 
with 
d= (Flac fi() —BG81(0) He (68) 
d= (F\-5 OA (1-2) 
+4 (G-2OK(OK LOO (6.0) 


-#(6- 2) AMOK Ose | 


Qén— € 
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€ 


ado,= —(f)ecs, (€) + 2OKO 


2€ 54 ate 


a 


~(6¢ +£C,)F,() 14 K(6) +2L(0)) He, (6.10) 


where 
(42)°CQ=g;' +324, 
(42)’C,=2,3'+ 6g3¢7—32,', 
(41)°C,= (g"—g5')’+3(e"—-£4)", 
(42)"C.= Gr —gs)gs +3 (Seg) gs +3 (br —Es)si—3(gr—ge)ge, (6.11) 
(40)°C,=g3 (6s —g1) + 3¢7(ei—a). 
(42)’Co= (g/—gs)"+ 6(e°—g5") (g¢—27) —3(e3—g4)*- 
With g;=7f and g;=0 (2¢*=3), the expression (6.7) agrees with the expression 


(6.2) as is expected. 


§ 7. Symmetrical vector meson 


* The interaction Hamiltonian density is 


A(x) =f 9 (x) rz (4) Gi. (2), (7.1) 


and in the first approximation of treating nucleons non-relativistically the cross 


-sec- 
tion vanishes. 


The differential cross-section is 
do= (1/m’) (2/42) af (©) —Re(€)] (7.2) 
with f(E)= SYA (6) and g(e)— > e(€). 


Functions 7,(€) and g,(€) are rel 


ated to those of symmetrical scalar mesons as 
follows : 


RO=(/OA), fl) =Be/e—ppn(6), 
A) ==), 


and similar relations between g'(€) and g*(€), using suffices » and 5 which indicate 
vector and scalar mesons, respectively. 

The angular distribution of mesons which are emitted by the effects of the 
nucleon’s recoil (f," and §:') is spherically symmetric unlike that part of scalar 
mesons. Since this effect is the most predominent among other contributions in 
the low energy region on account of a factor (3¢°/e?—42), we shall have in 
Practice a sherically symmetric angular di 
those obtained artificially. 


(7.3) 


stribution for low energy mesons as 
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§8 Discussion 


Comparison with experiments of the predicted angular distribution, energy 
spectrum and total cross-section for production of 7-mesons will throw light on 
certain properties of the z-mesons and the strength of its coupling with nucleons. 
But at present experimental deta are very few and inaccurate, so detailed analysis 
must be put off and here we satisfy ourselves only with summarizing and comparing 
our results each other. 

(1) Angular distribution 

The angular distribution curve of 7-mesons produced can be expressed in 
our approximation as follows : 


D(0) =1+x(€)P,(cos 9), (8.1) 


where z(€)=O for pseudoscalar mesons, x(€)~O for pseudovector and vector 
mesons and x(€) is a quantity of order 1 for scalar mesons. 

(2) Energy spectrum 

If we assume collided paiticles are at rest initially, then €,,~163 Mev for 
incident protons of energy~345 Mev. But as it is not the case, we make another 
assumption that nucleons in the target nucleus are moving in Fermi distribution. 
In consequence of these movements we shall have larger €,, than~163 Mev, for 
example, ~260 Mev for the nucleon which moves in the opposite direction to 
the incident proton with energy~25 Mev (maximum Fermi energy in carbon 
nucleus), i.e., With momentum~220 Mev. To compare directly our results with 
experimental data we must take the average over nucleons in the target nucleus. 
But here we only plot the predicted energy distribution curves** with €,,~240 
Mev together with the experimental curve*** (Fig. 1, 2, 3, 4 and 5). 

Near the lower end of the curve, ie., when €—p#<€,,—#, it increases in 
powers of (€—y)1” for pseudoscalar (A), pseudovector and vector mesons and 
in powers of (€—y)*” for pseudoscalar (B) and scalar mesons. Near the upper 
end it decreases in powers of {€,,—€)'? for any case. 

(3) Total cross-section 

The total cross-sections near threshold energy are plotted in Fig. 6, 7 and 8. 
They increase in powers of (€,—/4)* for pseudoscalar (A), pseudovector and 
vector mesons ard in powers of (€,,—/4)* for pseudoscalar (2) and scalar mesons. 


* If we consider other nucleons in the target nucleus, some transitions might be forbidden by 
Pauli principle. But as was shown by Feynman we have only to take account of this effect in 
intermediate states. For incident protons with energy~345 Mev recoil nucleons have energies 
larger than the maximum Fermi energy, so it will be admitted to neglect this effect. 

** For pseudoscalar (B) and for scalar mesons (B) we take Serber potential, i.e., gy Ar= (8¢ 
+ 3c) /16, £0°/4n= (—3c, +¢)/16 92 /An= (¢, —3cg)/16 and gy/4n= (— 4 —¢3)/16, where c,=0.280 
and c,=0 404. See referenc 

*** In this experiment meson. nich were emitted in the direction perpendicular to that of incident 


protons were observed, and in this case we may put ¢€~é, approximately. 
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Table II. Numerical values of total cross-sections 


Meson type Pseudoscalar Pseudovector Scalar Vector 


(72/47) 0.05 (2172/2)? 0.05 0.05 0.05 
Incident energy Total cross-section (cm*) 
(A) (4) 


345 Mev 7.0x10-27 4.7x10-*1 4.6x10-9 2.2x10-3! 1.1 10-% 
500 Mev 3.5x10-°6 1.3x10-29 4.7x10- 6.9x10-39 1.1 10-*% 


(These values are the average total cross-sections for production of charged 
mesons over the nucleons in the target nucleus having the equal numbers of 
protons and neutrons.) 


As a preliminary estimation we give some results ia Table II, assuming f*/4r~ 
0.05 (2m/f)* for pseudoscalar mesons and £?/4z7~0.05 for other mesons. 

(III) z+ and x- ratio 

If we assume that the incident nucleon is a proton and target nucleus contains 
equal numbers of protons and neutroas, then the 2~/z* ratio is given by 


k=1/3 for pseudoscalar and pseudovector mesons, 


R<1/3 for other cases which we have considered, because of the presence 
of f-term. 


Of cource, it is true that we must take into consideration the effects of Coulomb 
potential of the target nucleus for low energy mesons, but it is to be noticed 
that R may differ from 1/3 by the presence of §-term. 

In conclusion, we should like to express our sincerest thanks to Mr. Haya- 
kawa, Mr. Yamaguchi and Mr. Fujimoto for their kind interest in this work. 
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§ 1. Introduction 


Recently Koba® has proposed a new method of calculating S-matrix, using 
only Lagrangian but not Hamiltonian, i.e., in the conventional field theory the 
law of motion is described by Tomonaga-Schwinger equation based on the 
Hamiltonian formalism, while S-matrix which is the integral of Tomonaga- 
Schwinger equation and calculated by Dyson’s method requires only Lagrangian. 
Therefore, so long as we employ the integral formalism as the basis of field 
theory, we need only Lagrangian and the role of Hamiltonian is virtual. The 
proof will be given in the following sections. 


§ 2. The chronological operator P* 


In the following discussion we always use the interaction representation, the 
field quantities satisfy the same equations and commutation relations with these 
of free fields in the Heisenberg picture. 

Now the total Lagrangian density Lj... is decomposed into two parts: 


Liova= Tite “ Loins : (2. 1) 


We discuss only about Zin, so we simply write it as Z. Here Z is a Lorentz 
invariant Hermite operator allowed to involve the derivatives of the fie'd quanti- 
ties of any order, i.e., let ~,, uf be the field quantities, then Z is a polynomial 


or even a power series of the following quantities : 


Ou, OUg we ous Vous 
; » vets Uns — Hy peers 
OX, O4,0%,y Ot, A402, 


Uns 


Using such Lagrangian, we make an operator : 
U [o, oo|=DI0,[¢; dul, (2.2) 


where 


Uo, 6]= “(<) (itr. f dy PR(L(4))+-L(4))- (2.3) 


n! 
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Now we must define the notation P* which is the modification of Dyson’s P. 
But as it is hard to define it directly, we define all the matrix elements ot 
RAUL (4 )ice ed (Aq Os 

{n a quite similar way to the Dyson’s diagram method, the matrix element 
of any process is given by the combination of the following terms : 


Pi ad el oe eae Nie (2.4) 


’ 
, 
OF agen OP yy OF" yieasOE Vk 


so that the general definition of P* is reduced to the definition of (2.4), which 
is given by 


< P* 4 ibialneh med ee 8) F 
OF nseen tps. | Oh panet wy 
gtk 


z OX ys-+-O% ps OF" y4-..02' yy <P(iter #5.) > o- (2.5) 


By this definition, the operator P* commutes with the differential operator in 


contrast with the ordinary F, this is just the same with that proposed by Koba™ 
and Matthews.” 


P is determined by giving a family of space-like surfaces, and so long as 


o(x)=-0(2"), e(%, x') commutes with the differential operator, so that in this 
case we know 


APG “ppt Se Jeu 


More generally speaking, if all the surfaces o(x,) differ from each other, for 
instance 


a(%;)>o(4.)>...>0(x,), 


(future> past) 
then P and P* coincide with each other: 
PE(L (41) + L(4n)) = PUL (4) LZ (an). (2.6) 
In the same way, if 
a,>¢,, (s=1,--+,4; t=£+1,-++,n). 
then P* is decomposed into 
PRL) +L (Hn) = P*(L (ay) LZ (a4) P* (L(t) LZ (ay)). (2.7) 


What we must notice here is the fact that the quantities such as <Pu Was 


>» are represented by 4,y, Sy, or their derivatives only which are independent 
of the choice of the family, so we can conclude: 
Ole, a], consequently U 


, [¢, a] are determined uniquely by giving two space- 
like surfaces o and o. query. Oy etving pace 
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§ 3. The generalized transformation function 


Now suppose that the surface o is deformed into o’ by an infinitesimal volume 
dw at a point # on o, then 


U,{0’, a,|—U,|o, oy | 


nN ik > \n for 
= (Dar) lite [ee fete PE) LG) 2) 
+ (5) =(4) \ae, [dr. [aes : {tr,P* (Le) LG.) LG,)) 
M: ce o Go Go 
If we employ [o,, o,) as the integration domain of le then [s,, ¢,) [o, o’) have 
O1 


a 
no common point, and so P* is separated by (2.7), thus we get 


U,[o’, a|—U,[¢, 9%] 
= [4 [2ed¢n | oy; (=) [aee--[ar,P* (L(x) Law) | 
+{5-(4) (ze, [ae Pr(La)L(2)) | oy, (4) 


fitzy-(dry PPL) LZ a) | 


| 

-[ifdteren pote ote[ 2 (G) es fb mcen ee Poder 

ae (3.1) 
We put 


Lim. {2 [Z@)ax,J=-4 Tie), 


tim, aL a) Ji fe BEDE) [= gH, 


eee eee seseseese 


lim, — oH ree zy ae as, PCL ss L(2,)) [=—4H, C2 ele eae) 


In general, H,,:-- do not vanish,” because of the existence of 4-dimensional 
d-functions in the integrands, and as the integrands of (3.2) are family indc- 
pendent, A’s are independent of the method of limiting process o’—a and com- 


pletely determined by giving x and a oaly. 
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Inserting (3.2) into (3.1), one obtains 


U,l¢, a]=—4 {A Ualo, o.|+ H(4)Un-s[, o,|+ es }. (3.3) 
da(x) h 
Now put 
H[x; o]=3) Ha(2)- (3.4) 
Then summing (3.3) up with respect to ”, we get 
6 ee : 7 35 
———U [¢, |= ——Al*z; € a, 9%. (3.5) 
sig! le l= GF HLe; 6] Ue % 
Next in the same way with (3.3), one can derive 
= A al=5{Uasle, o.|H, (4) + U6, %]H2(40) + of. (3.6) 
d0,(%») h 
Comparing (3.2) with (2.3), we see 
H,,( x) =th lim. — U, [o’, o], (r=1, 2,-++) (3.7) 
Sela af 
especially 
H,(4)=—L(2), (3.8) 


but the definition (3.7) is not symmetric with respect to future and past, so that 
we can equivalently define it as follows: 


Hx) == tim + | Ue, o]+Ufo, a”). (3.9) 
oll>o 


Then its functional derivative is given by 


: -F, (% mes fim. | Ue a] , dU fo, o”] 
(2’) 2 a 20 dw a(x’) 3a(2’) 


ba 


inseiting (3.3) and (3.6) into above equation, one easily obtains 


Toy ig C2) ota 2+ (2) the (2) 4 PR \. (3.10) 


This is the fundamental equation of our theory. 


If fy Hy...,H,-1 are Hermitian, then from (3.10) H; is also Hermitian, while 
tive — J, ic Hermitian, so that all H’’s are Hermitian by the mathematical induction. 
Next summing (3.10) up with respect to 7, we get 


7) a : 
doa HUE l=) | Ales 0), He"; 2] (3-11) 
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This directly guarantees that the Hamiltonian (3.4) actually satisfies the integra- 
bility condition. 
As A[z; o] is Hermitian, UV [o, a] is unitary from (3.5). 
Thus translating all the above results into the p.operty of U [o’, a], we get: 
The property of U[o’ ,o]. 
(1) U[e’, oe] is an unitary operator uniquely deteimined by giving two 
space-like surfaces o and a’. 
@)* 7 la,.o|=1. 
Gs) Ula", a | U0 feF'e)|. =o", “c]. 
Clearly these results show that U [o’, a] is the generalized transformation function, 
and S-matrix is readily given by 


S=U([o,—o]. (3.12) 
§4. Density matrix 


In this section, we will discuss the property of the density matrix in the 
interaction representation. 
Consider the generalized Schroedinger equation : 


go lela) Bio), (4.1) 
0a (x) 
and expand the state vector ¥[o], 
P[ol=S'a,[7] Pr; (4.2) 


where @,, refer to some eigenvectors which are independent of ¢, and satisfy 
Coz, ®,.) =9nns 


and the coefficients a,.[o] are c-numbers. 
Inserting (4.2) into (4.1), we get immediately 


ik ee Sun (2) anf], 


(4.3) 


. Oda t| Fak ae 
| ele) 2 Aye lawl) 


where 
Ae (Oa ea) 
Now we define the density matrix as follows: 
0(a, 0) =(a,lo]. axlo’])- (4.4) 
Then from (4.3), p(e, a’) satisfies the following equations : 
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nad ten 30 p(o, o)=H(z)p(a, a), 


ooh (4.5) 
dah o(a, o)=—p(o, o) H(z’). 
The expectation value of a physical quantity Q is given by 
(Q)=(¥Flo], OF[o])=Sazle] aslo] Onn 
=S,(p(<, 2) 9). 
We write simply 
p(s, 7)=p(9), 
then 
(Q)=5,(p(2)Q). (4.6) 
Next from the normalization condition : 
(¥[e], Plo})=Sax{e] a,fo]=1, 
one easily obtains 
(eC, (2, 2”) w= Dale] a[o’] afo’] ai[o"] 
=a,o] af[o’’], 
mt (a, 0") plo’, o”)=p(a, o”). (4.7) 
especially 
p(a)°=e(a). (4.8) 
Solving (4.5), p(¢, 0’) is given by 
(a, &)=U[a, a] p(a) Uap, 0’), (4.9) 


which describes the law of motion for the density matrix by means of the 
generalized transformation function. 


If we write for simplicity as 


p(a) = Po U[e, a |= S; 


then (4.6) is: rewritten as follows : 


(2)=S,(SpoS“* Q). (4.10) 
The relation (4.8) means that the Hermit 


0 and 1, thus we know the connection between the generalized transformation 


function and the density matrix, but here the concept of the density matrix is 
differently used from the ordinary one, 


€ operator p(¢) has the eigenvalues 
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§5. Integral formalism 


As has been proved in the previous section, only Lagrangian is required to 
calculate the generalized transformation function, and equivalent Hamiltonian always 
exists in principle, but we need only the existence of it, and not the concrete 
form. 

Therefore we can start from the integral formalism, so we will translate the 
differential laws into the integral laws. 

(1) The law of motion. 

[zo+4 on Pll= 0, (14.) 
z da(%) 
P[a]=VU[o, a]? [a]. (1i.) 

(2) The conservation law. 

Let Qfe] be a conseivative quantity, then the following condition must be satisfied : 


[z@+4 7 Bola awe’ a) |-0, (2d.) 


Qia]U\e, o,J=VU[e, a] Ola]. (2i.) 
(3) The supplementary condition. 
In order that the supplementary condition : 
Ela] V[o]=0 


does not contradict fo the law of motion, the following condition must be satisfied ; 


[ Roan tol l= - 
| 2) +— dale) 5 1 0, (3d.) 
E[o]U[e, 9,]=U[e, 0] =[o]. (3i.) 


An example of the interaction which is easier to introduce in the integral 
formalism than in the differential formalism is of a type as follows: 


L(z)={ fMO)(4)} Ar) 
=\7.(e) 7a 2) A,(*) as. 


where 


(4) = {A- — ke ye*er dk. 


Gr (27)! 


Thus even the non-local interaction can be treated in a covariant way in this 


integral formalism. 
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Appendix 
In the text, we see that the Lagrangian method is more convenient than 


the Hamiltonian method for the calculation of S matrix, and that practically the 
former is the only one when canonical description is impossible. So it seems 
necessary to clarify the role of Lagrangian more deeply. 

In the canonical formalism, Lagrangian is used to derive the field equations 
from the variation principle and also to quantize the field variables. While in our 
theory, the unitary use of total Lagrangian is impossible. For instance, the free 
part of the Lagrangian is utilized to quantize the field variables and to derive 
the equations that they satisfy in the same way with the canonical case. 

Such a procedure is always possible since the field quantities in the interaction 
representation coincide with those of free fields in the Heisenberg representation 
whiche an be described in the canonical form. 

Concretely speaking, we employ the representation in which field operators 
satisfy the following variation principle, 


2| Ze dx=0. (A.1) 


Thus the role of the free part ia the interaction representation is completely 
equivalent to the total Lagrangian of the canonical formalism in the Heisenberg 
representation. 

On the other hand, the interaction part of Lagrangian is considered to be 
the generating function of S matrix or the generalized transformation function 
that rule the causality. 

The relation between the free and interaction parts of Lagrangian resembles 
to that between the tram and the rail. 


On this stand point of view, S matrix is a functional of the interaction pait 
of Lagrangian 


py = >i. -) fr. [dr P*(L(2,)...L.(2,)). (A.2) 


The employment of the interaction representation is essential in our generalized 
theory. 


In the canonical formalism, field equations in the Heisenberg representation 
can be derived in two ways. 


The one is the variation principle : 
r) Baer dtg==). (A.3) 


which is independent of the canonical formalism. 
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While the other is based on the field quantization, after which we can cons- 
truct the total Hamiltonian by the common rule of the canonical formalism. Then 
the equation of motion for an arbitrary field quantity Q is given by 

dO 00. 3 
= + 
at Ot Ok | 


FLrorat » Q). (A.4) 


As the characteristic feature of the canonical formalism, field equations derived 
in the both methods coincide. 

But it is not true in our general case, as has been noticed by Koba.” Our 
method is equivalent to the second one essentially, the quantization, however, is 
referred only to the free part, but not the total Lagrangian. 

Also the interaction part of Hamiltonian cannot be derived automatically from 
Ligrangian in contrast to the canonical case, for it should be constructed by the 
requirement of the integrability condition. 

Corresponding to (A.4), we have 


faa) OE) a a. <+{4S,! Hud’; 2, Q(z)>, (A) 
ax, OF, i 

where < > means the expectation value on a surface a, and field variable Q(x) 
is expressed in the interaction representation. Ajm,[%; o] is non-local, but the 
spread is confined within the surface o@, so that the initial condition on a surface 
% is enough to determine the future behaviour of the fields system in agreement 
with the causality principle. 

Whlie the field equations derived from the variation principle (A.3) in the 
non-canonical case have the non-local interaction which spreads not only space- 
likely on a, but also time-likely,® so that the knowledge on the initial surface 
7 is no more enough, and the future behaviour depends also on the past history 
in contrast to the basic causality principle of quantum mechanics. 

Thus the unitary description by the total Lagrangian is broken and the 
Variation principle does no more hold in our non-canonical theory. 

This is a very gloomy conclusion, for theorems which can be derived from 
the variation principle have no validity and they should be demonstrated by another 
cirect method. 

For instance, gauge invariance is no more clear from the form of the total 
Lagrangian, since it requires an intimate connection between the free and interac- 
tion parts of Lagrangian. 
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§ 1. Introduction 


We may define distribution functions for given number, %, of molecules as 
being proportional to the probability density of finding 2 molecules in a certain 
configurational position, averaged over all positions of the other molecules of the 
system. These functions are the exponent of minus the potential of average 
force between molecules, devided by £7. Mayer” derived a set of integral equa- 
tions which relate variations in the potentials of average force at two different 
fugacities. On the other hand, Kirkwood” evolved an integral equation for the 
radial distribution function, and Born and Green® obtained a similar equation. But 
the relations between these integral equations and Mayer’s general variational 
method are apparently not known. The purpose of this paper is to attempt to 
clarify these relations, at the same time, to derive new integral equation which 
is a generalization of the Borh-Green equation. We will deal with a multicom- 
ponent system, because the final equations are in no way more complicated if a 


system of many components is considered than if a single component system is 
handled. 


§ 2. Symbolism 


In treating a problem involving the coordinates of an extremely large number 
of molecules it will be necessary to use a system of symbols originally introduced 
by Maye1” in order to attain some brevity of expression. Let us consider an 
infinite system containing o types of molecular species, and use the subscript s 


(or ¢) to refer to a special species of type s(or ¢#). The symbols without subscripts 
refer to a set of quantities 


P=Pu Poy tenes »Ps» seeeee sPos (1) 
s=21, Soy teens yFagrrree® sot (2) 


The symbol p, indicates the density of molecules of species s, in molecules per 
unit volume, and 2, is the fugacity normalized in density unit. 
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The letters 2, m, m, etc., when written without subscripts, will be used to 
indicate sets of integers: 


R= hy, Boyrreee Eisai Fos (3) 
of #, molecules of type s. The short hand notation: 
A=! hl---Byl---Ze!, (4) 
f=)" 0.5 -+ 058s +++ Pero (5) 
gaz," ™ seeg ks 25% (6) 
will be adopted. 
The symbol 
(2,) =Aisy Visy Ztsy Qitsr Jriss’** (7) 


indicates the configurational coordinates of the i-th molecule of type s, and in 
general will include internal coordinates, Jiisy Yoers***,aS Well as the Cartesian coor- 
dinates, %4, Js, 2%, of the center of mass. The coordinates of 4, molecules of type 
s will be represented by 


(A$=(1), (2),+(4). (8) 


And we will denote the coordinates of all molecules of such a set & by 


(Aj = {Ay}, CAT LAP Ce (9) 


The volume elements are similarly expressed by d(é4), d{%,}, and d {fh}. 

A function, %,{2{, of the coordinates of the set x of molecules can be ex- 
pressed as a sum of the component functions ¢,{»}, which are functions of the 
coordinates {¥}, of subsets of molecules out of the set ”. These component 
functions are defined by either of two reciprocal sets of equations 


Y,= (Sin) Av (10) 

or 
f= Caiz},) (—)"™ fv. (11) 

where the summation symbol, (3}{»{,), indicates summation of functions of such 


subsets, {v},, over all subsets, for all subsets, for all values of »v, 0 <», < x,, 
including a single term for which all »,=0. 


§ 3. Variational method 


With this notation we now define functions F,(p, 7, {x}), of the coordinates 
of a set 2 of molecules, by the statement that in a system of infinite extent. 
(V>n,/p, for all s) at the temperature 7, and density p, the probability that there 
are an appzopriate set 2 of molecules at the coordinates {”} within the configura- 
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tion element djz} is equal to p"Fn(p,T,\2}) a {x}. It is to be noted that 
with this definition the functions F,, are so normalized that, 


in| ff fee ats [=v m™ 


Mayer? and McMiilan” have shown tht for two sets of densities ps and Py, for 
which the pressures are fy and #£,, and for which the fugacities are 2, and 43, 
respectively, the distribution functioas are related by two sets of equations : 


Gr=T (| [Chemin (13) 
m>0 yn'\ Vv 
or 
Grae "||. (eaem djmt, (14) 
m>o yy! V 
where 
V 
Gefen tel [22] Alen Tie), 
exp Hee || 2 | Fale Tim) (15) 
Cra lex AV | & | R 
[exp 47 BT Flom Ts im), (16) 
and 


-— sa ~ 


We will consider two sets of functions G,(y) and Gi(y) which depend upon a 
parameter y. And we now assume that Gi(%))=G%s for a certain value of the 
parameter, 7, and that functions G,(yv) depend on y only through the equetion 


fod 
~ 


G.())=% 


m>0 772! 


Ye: [oRen (9) fom). (17) 
Then, from (15) and (16), we have G,(%)=Ga, and the reciprocal of (17) is 
CD=DCO” [f-- [Carn O) a it. (18) 


It is necessary to assume that G,,(y) and GX (y) are symmetrical in permuta- 
tions of the coordinates of identical molecules except the coordinates of a certain 


molecule z,. Then we cannot defi i 
: ne unsymmetrical G and Gx 
the set 7) contains the molecule ty a) 7) 


We now define the functions ¥, and ¥X by 


y= In GaCy) 


wy (19) 


Integral Equations between Ditribution Functions of Molecules 
ay 
and the component functions ¢, and ¢¥ by (11) and 
y= (2{2},)(—)"" Fe, 
so that (10), and, 
Pr= (21a oy 


are valid. From these equations Mayer obtained 


du Di (Diahs) [fo [btm Kemmd fre), 


m>0 771\ Vv 


in which the kernels, Ay.x,m, are defined by the equation: 


Ky snnm= 2 {x i) Ge a ees ee 


By the same operation he obtained 


gt= DK)” (Sria},) || [ born KEnm dtm, 


>! m!\ 


with 


Ge AS (@2, {x}x) Oe an 
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(19’) 


(20) 


(21) 


(22) 


(23) 


(24) 


(25) 


§ 4. Differentiation with respect to the coordinate of a molecule 


Let us consider an arbitrary function, 2,{7}, of the coordinates of the set 
2 of molecules. And for the sake of covenience we will represent %,{} by 
P,(2,), P.,( (4s) (7)) and Puy (C45) (12) (%,)) when n= (i,), w= (2) (7) and n= A) 
(7) (2) respectively. Kz+nn,m Will be denoted by K(jk+n}; {2}; {m}). The 


potentials of average force between molecules 7 are defined by 
W,§nt{=—AT In F, (fa, T; {nt), 
W*in{=—AT In Fa(s, T,}2})- 
The component potentials, w,}7}, are defined by the equations 
wy {y3=(oint(—)"™" Mains 
w *{vj=(Din}v)(—)*"Waiah. 
If differentiation is performed with respect to the coordinates (¢,) 


molecule of type s, we have, according to (15), (19) and (26), 


; 1 aw, 
7 i) = aye amit, 


(26) 


(26) 


(27) 


(27) 
of the z-th 


(28) 
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A awh) GD), (29) 


P ,,( (4) ODS OX te 


since p and p do not depend on (Z,). 


From (28) and (29), it will be easily seen the functions ¥ 4" are the 
averave forces. Y,’s and ¥Y*’s are zero unless z contains the 7-th molecule of 
c € > we n ‘ 
type s. And from (10) and (27) we obtain 


5 (i) —-— 1 9u, (4) 30) 
J.(74) erage ( 
(i) (G))=— F-Pt G) GD) (31) 
And similarly we obtain 

WK (4 ee ir" Ow F(z.) 39 
YOO PET ae, (2) 
eine tal Bel 

bi() W))=—s7 iy (33) 


The function W,,}} for three molecules approach the sum of the single molecule 
functions when all distances between molecules are large. And Kirkwood?” 
assumed that for small distances they are given to good approximation by the 
sum of the functions for all the pairs. This assumptioa is certainly not exact 
but used by Born, Green® and Rodriguez” And it is equivalent to the assum- 
ption that wai}, wax}a}, gia} and Pin} are zero except for the case where 
the set 2 contains only one or two molecules. 

If we introduce the Kirkwood approximation mentioned above in (22), and 
if the asterisks refer to the fluid or glassy system, for which w*(i%) and $X(4) 
are zero for all s, we obtain the following integral equations : 


HU=Vee [EGE Ks Gi HAE). BA) 


PG) (ID) = Fal Ge) G)) # EGE) RRC: GMA) A(E). (85) 
Here the kernels are given by 


K( (7) ; (25) :(4,) ) = G.( (4s) (&,)) 


G, (2,) : Ce 


and 


K((4) (W034) ;(4,)) = Sr LED DAD) G(s) (2) 
CCOLGAHCARCS)) Gai) Gd) nic aan (37) 
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Inserting (36) in (34) and (87) in (35) and subtracting (34) from (35), we 


(i) Ey) ePID ACE). 
(38) 


have 


Yul (is) (J) — Glin) = GE (@)(W)) +4 fg 


(34) and (38) can be, according to (15), rewritten in the forms 
i) = Zl Cee Foe) las) (4) eC) a4), (89) 
and 
Yat (te) (F1)) — Poe) = Pe (4%) (Ze) 
>) Para Fu) *((4,) (& Fete (4s) (He) Ar) ) ab 40 
FE OO may 2 
It will be easily seen from (26) and (27) that the Kirkwood approximation 
(Ci) 4) Fa) GD), (41) 


can be expressed in the form 
; — Fn (Gs) (Ki) ) Fer 
| pee e t k,, — 
BD Ee) FG EG)F AE) 


Substitution of (41) into (40) leads to 
(TAC CALC Ay bea b (25) = pe ((z,) (7:)) 
P*((, ) (k,)) Lh CHS Aas Ce) ) ad ) 


(2k, 
(42) 


ce Par( Zar— Zar) 
y aeaail 


at Zar 


§5. The Born-Green equation 


Let us consider the case where the both values of the fugacity, 7, and 2 
. Then ¢,(2,) and w,(z,) are 


correspond to the gaseous, liquid or glassy states 
In this case (42) can be written in the form 


Zero. 
OF, ((2,) C7.) Fh @OD) D2 ( (24) a). 
OX 
Fa hie -) (a, FG 1 ) 4, )) al Wye (4) Gy d(kr), 
OX py 
(43) 


OXts 
= 5) Peeler 200) £4((4) (4) | 
where w,((z,)(7,)) is assumed to be a function of (z,) and (7) only through 


i] 


the distance between molecules 7, and 7, 
For the system composed of pure substance (43) becomes 
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aF,((z) (7)) + FOO) aWF¥(()(7)) 
“a EE Ox; 


: (FCA A)ACAA)) IWS COA) 1 py (44 
2 ere Fp) [FXO OFA) oi Pr tae U(h).(44) 


Here we denote F,, and w,, for a single component system by F; aan WW, TeS- 
pectively. If the asterisks refer to infinite dilution, z.((z)(7)) is equal to 
u(z, 7), the potential between the pair of molecules zand 7. Then (44) becomes 


AF) (7)) 4 COCA) ou, 1) 


On; kT OX; 
, 4 2 Bs k 5 
=pF,( (2) (7) 1 (7) Pz (J)(Z)) 9 a5 ) dk), (4 5) 


in which p=ps— fs; ard is equal to p. if ps=O. This integral equatioa was fast 
derived by Born and Gieen®. Hence we conclude that (43) is the generalization 
of the Bron-Green equation, and a statistical-mechanical theory of fluid miature 
will be evolved, based on the integral equation (43). 


If g,’s refer to infinite dilution and ¢,’s refer to solid state, we can write 
(39) in the form 


OF, (4,) a= = F,,.((@,) (4,)) Ou,,.( 255 hy hb 46 
eet: [22\Ge ie ath), (46) 


in which ~,,(z,, 7,) is the potential energy of a pair of molecules of types s and 
t, 4, and 7,; and for single component system 


OF (i) _ [Fel (i) (@)) Bui, &) yp 
ay, =e LT Pals 


(47) 


Eq. (47) resembles the equation used by Kirkwood and Monroe” in their theory 


of fusion, and is essentially equivalent to the equation used by Osawa and Sato” 
to discuss the multimolecular adsorption. 


& 


§ 6. The Kirkwood equation for pair distribution function 


In the presert section we asseme that the asterisks refer to infiaite dilution, 
and the functions GX(y) have the form 


Ge (fat, 7) =exp(— Un( na y) (48) 
where 


Yalta}, N= RMN ly ji) + (I= LDtalin Js (49) 


p=lt= 


the potential energy of the system being assumed to be expressible zs a sum of 
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terms corresponding to the potentials of pairs of molecules. Since U,({x}, 1) is 
equal to the potential energy of molecules x, U,{n\}, we have 


GE(jn}, 1) =exp(—Z!4), (50) 
kT 
which is equal to G*}x} for ¢=0. And from (33), we obtain 
$8(G) W))=— ee 0. (51) 
k7 
W((4,) (Z,))’s are zero unless either 4, or |, is ¢. And in this case Gx {n’s 
aye zero unless 2 contains two and only two molecules: Then we obtain, from 


(15), (22) and (51), 


PA (>) C7,)) 4 ai ts Gs ji) 


2 Gy. By tel Ge) Ge) Ee) 9) _ Fr Gn) hd» ID Lad 52 
* ni eet A] FAG) G) I) FAG) 1) feo ve 


by the similar procedure as used in the derivation of (34). According to (11) 
(15), (19) and (26), ¢4((4)(j)) is given as 


. wars. 1 B00 ( (2) (i) IY) De 
Ps ( (25) (f)) = LT ay i 2) 


Inserting (53) in the left-hand side of (52), and replacing /, by tn», we have 


O74 ( (7) ( fe)» ¥) 


ay =Ug(iss cB, 
2 4, i? qi he) —W, hs P ’ 
Sip fats [exo {— aE). = WEG). D)} 


Wr bs En ’ — W,((is), ) rs 
—exp{- (GG) DATA DY) \ate,). (54) 


Here 2({2},y)’s are defined by (15) and (26), if one replaces Gy by Ga(y)- 
Since W,,( (4) (jz) (Ar) 91) — W,,((is) (fr) Y) is the potential of average force upon 
h,, when 7, and 7, held fixed and W,,( (4s) (2), Vy - WG), 7) is the potential of 
average force upon /, when 2, held fixed. Eq. (54) is the one which Kirkwood 
has already derived. 

Here we will consider a fluid or glassy system for which 7w,((%,), 7) is zero 
at least for y=1. Bearing this in mind, expressing W,({2},7) in terms of 
Wn({2}, 7) and using the Kirkwood approximation mentioned in § 4, we integrate 


(54) between zero and y with the result: 
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2054 (45) (Je) » I) =I (Ze) (Fe) ) 
_ Weel (45) Ar), PP exp| _ Worl We) ) a: 1 
kT 
a(k,) ay. (55) 
Here we assume that p, and w,((7)(4,)) for 747, ts, do not depend on y, 


since the variations of these quantities with y are of completely negligible order, 
O(1/n,). And for a single component system (55) becomes 


+3." Won fae h,.) exp{ 
r=1 0 


a0,((2) (7), ¥) =I2 (4, 7) 
Wo((2) (2%), v))P 7((2)(& a 
+/\" fasG, 4) exp{— D2) exp{—- 2X OO) Yt—1|a(eyay. (66) 
Eq.(56) has been utilized by Kirkwood, Monroe® and the present author in 


order to calculate the pair distribution function of liquid argon from the interatomic 
force. 


§7. The Kirkwood equation for single-molecule function 
Let us consider a solid system, for which '.(7,) and w,(z#,) are not zero. 
And, in this section, we will confine ourselves to the case of single Component 
system to avoid complicated arguments. If U,({2},7) depends on y in such a 


manner that given by (47), and if the asterisks refer to infinite dilution, we ob- 
tain, from (22) and (23), 


- 2)¥2((2) (2), 7) 
h+ W@=pfuG, JZOON ag, o7 
v Ay(@) yw) a 
According to (15), (16) and (26), ¢,(2) is expressed by 
JN=v y= 970%), 9Ilnp_ 1 Bw,(%) 
Got, (=a? Ga of ws 5 SP Sete ; 
in AL ey ee ay kT ay 
in which p(y) is defined by 
pLOV = {I--fex {- Umm, 9) 
exp yi a aa Na JexP i bd fan}. (59) 
We introduce (58) into (57) and integrate it between zero and y, obtaining 


Vi £0)—P0)! 4 nC) _ (Ci), 9) 


kT pte) kT 
Pad 8 Bey p) £2 (2) (4), 9), we 
f fe “(z, &) F(),9) W(h)ay, (60) 


in which the const 
vanishes for y=0, 


According to McMillan and Mayer”, we have 


ant of integration is determined by the fact that z,((2), 7”) 
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(ea) Sy yp ie. (61) 
Os / ny g 
Differentiating (59) with respect to z, and making use of (61), we obtain 


y fexp ? OY} Pam # | J--fexp{—ZenGZt2 Otatmjac) (62) 


fm 
gs m>0 772! rv 


Then we differentiate (62) with repect to y, obtaining 


| expe Zz oe aL py) V+4T In p(y)] 
kT bile ay 


ce {f--{[-4 oa OI) Jexp|— anna Ora 
dimid(z). (63) 


“M 


m>0 770! 


I 


And if we use (49), (63) can be rewritten in the form 


Voor: 0 [expt | al py) V+AT In p(y)] 


AT Oy 
== 2 Ll. Lara in z,exp|— Toshi 404 I) aim sa@a(e). (64) 
AT m>0 mm) JIV k1 
Using (15), (17) and (48), we now obtain 
YL ag) VAT n ON) — BN ne, HELO, NAOH). (65) 
AT ay TY 


Now let us define U(y) by 


Tin= Effi HOW Aad, (66) 
Then we obtain from (65) 


Here p(1)— p(O) may be negligible compared with p(1) or p(0). Bearing the 
above mentioned in mind and using (26) and (66), we can rewrite (68) in the 


following form: 


rw Es =a) (2) )F, (2) dh). (69) 


y 


where 
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= SA) Peas: 7) (£), )d(k)dy, (70) 
In aire (i, BF, (() (2), yaad (Bay 
\(k))=—L ni, b z ; 71) 
M((G)(2))=— ali, BL) (2) ( 
LA) ))=( OQ Hvay, (72) 
and 

(2), 9p) =< FLOM”) (73) 

(OOD-FO, DRO 


The non-linear integral equation (69) has been used by Kirkwood and Monroe® 
in the theory of fusion. 


§ 8. Conclusion 


According to Mayer’s general varitional method, we consider possible varia- 
tions of the: potential of average force. Then we obtain the Born-Green integral 
equation and the Kirkwood integral equation for the pair distribution function 
if the Kirkwood approximation is introduced. 

As shown in § 5, this variational method shows that the generalization of the 
Born-Green equation to multicomponent systems and the one to the case where 
both of the fugacity sets for which two sets of the pair distribution functions 
are related are not zero are possible. This method will play an impor- 
tant role in the theory of liquid solutions. And this will be used to make clear 
the origin of the discrepancy between Mayer's statement and Rodriguez’s one, 


as to the relation between the singularity in an isothermal and the equilibrium 
of liquid and gas. 
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1. It is well-established both theoretically and experimentally for the crystal 
- electrons that there exist the rather predominant interactions between the lattice 
vibrations and the electron in crystal, which play an important role in elucidating 
the temperature-dependent properties of the crystals. For example, the absorbed 
energy by crystal electrons from the incident radiation cannot usually be reemitted 
except for some of the specially treated crystals, i.e., crystalphosphors, coloured 
alkali halides etc., but may be transformed into the lattice vibration energy, 
thus resulting in the appearance of the characteristic absorption alone. From 
the same reasoning the mentioned interaction should be expected clearly to 
exhibit itself in the phenomena of the temperature-shift of electronic energy band 
in crystals, which will be experimentally observable, owing to the approximately 
continuous distributions within each of the energy band, at both ends of an energy 
band separated from each other by the forbidden region of finite breadth charac- 
teristic for both insulators and semiconductors. From the theoretical view point, 
since the temperature shifts of the electronic energy band are clearly ascribed 
to the reactive effect of the interactions with lattice vibrations of crystal electrons, 
the mentioned shifts of the edges and of energy band may be considered to be 
consisted of two kinds of terms, i.c., the “ self-energy” term and the “ radiative 
damping ” one, similarly to the case of the electromagnetic interaction with free 
or bound electrons.” The former is found to give rise to the finite amount of 
energy shift due to the finite number of degrees of freedom of crystal lattice, 
while the latter to the level-broadening specified by the half-value breadth of the 
energy level, which states of affairs may be represented through the wave 
function (1, ¢) of the crystal electrons in the. foliowing manner ; 


g(r, i) =exp{ —i(Z, + 0E) t/4—P'thexp (kr), 
in which &, expresses the electronic energy without taking into account the 
interactions with lattice vibrations, 62 the level shift of the “ self energy” type, 
and &l’ the half-value breadth due to the “ radiative damping.” Further, the wave 


number vector & specifies each of the electronic levels in crystal. Since the energy 
distribution of lattice quanta is well known to be governed by Planck’s law at 


834 T. Muto and S. Oyama 


thermodynamical equilibrium, both quantities of d# and 4/” should be dependent, 
to a more or less amount, upon the temperature of crystal. 

It is already observed®’® experimentally that in many crystalline insulators 
there is a displacement of absorption toward longer wave lengths with increasing 
temperature, which-phenomena is simply manifested in many cases by the deep- 
ening of the colour of crystal. On the other hand, the temperature-shift of the 
absorption toward the red side has been theoretically discussed in detail by 
Radkowsky® by taking into account the level broadening due to the radiative 
damping effect through the interaction with lattice vibrations, which effect seems 
to give a successful explanation for the mentioned shift in polar crystals since 
the thermal expansion effect is too small by a factor of at least four as shown by 
Méglich and Rompe.® Recently the temperature-shift in non-polar crystals (Si) 
has been observed experimentally by Pearson and Bardeen™ thiough the detailed 
analysis of the changes with temperature of the electrical resistance and the Hall 
constant of the mentioned crystals, which results are shown easily, by a simple 
calculation, to become too large to be explained by the level broadening in 
Radkowsky’s theory. In the present paper we shall show that the observed 
temperature-shift of the energy level in non-polar crystal may be successfully 
attributable to the “self-energy” effect through the interactions with lattice 
vibrations and, further, the level broadening in polar crystals becomes of com- 
parable oider of magnitudes with the “ self-energy” effect, thus the total con- 
tribution of both effects being responsible for the level shift of the mentioned 
crystals. 


2. Hamiltonian of an electron under the combined influence of both the 
electrostatic field with crystal period and the lattice vibration may be written as: 
follows. 


For the semimetal (non-polar monoatomic crystal), we have 


k 1 ° 9° 9 4 
“iia a @ AE (te) +42°%(w) of (w)}—R grad V(r), (1) 
a Ne a ee peta 
Aah H, 0H 
R=—. > pm €.{a,(w)exp (iter) +.a,* (w)exp(—iwr) }, (2) 
in which 7, =Hamiltonian of an electron in the crystal field, 
A= Is of the lattice vibrations, 
6H= - of the interaction with the lattice vibrations of an 
electron, 
(1, 2,3) =Suffix specifying each of three modes of vibrations with different 
polarization, 


€,=Unit vector of polarization of elastic waves denoted by 4, 
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fu 4 =Canonical momentum and coordinate of a normal mode of vibration 
corresponding to the wave number vector w and polarization 4, 
y,(w)=Frequency of a normal mode of vibration belonging to w and 2, 
M=Mass of an atom situated at each lattice point, 
N=Number of atoms involved in the crystal. 
Further, a,(w) and a,*(w) express, as usual, the operators responsible for 
the absorption and the emission of a phonon, whose non-vanishing matrix 
elements. are 


(W,(w) +1ar*(w) |W.) =/ 4 _(H,() +1), 


dnv,(w)* 


eee (3) 
N,(w0) —1 Ja, (w) |7,(e0)) =,/ 4 : 
(Vw) la, (w)| V,(w)) tay, (BY ee 
N,(w) represents a phonon-number, the average value of which is given, at 
thermodynamical equilibrium, by Planck’s law 


NV, (w) = {exp (Zy,(w)/h7) —1}", (4) 


k, and T being Boltzmann constant and absolute temperature. 

For the ionic crystal (diatomic polar crystal),”® as shown by Born and 
v. Karman, the normal modes of the elastic vibrations of the mentioned crystals 
can be divided up into two branches, the acoustic and the optical one. The 
number of normal modes in each branch is 3V when 2 is the number of ions. 
The oscillations of the optical branch are those in which two neighbouring ions 
of opposite sign vibrate in opposite directions, while, in the acoustical branch, 
neighbouring ions of opposite signs have almost the same displacement. Since 
the optical branch gives much the greater polatization of the lattice, it gives a 
much stronger interaction with the electron. So we shall be allowed, in the 
following, to take into account only the optical vibration in the mutual interaction 
with crystal electrons except for the case of extremely low temperature where 
the contribution of the acoustical branch to the interaction energy may be rather 
appreciable, owing to the increasing phonon number, in comparison with that of 
the optical branch. 


Now we have 
C Keats 
Vice ake : Ine yh { in(wr) +22 *cos(wr , 
Hair BL eB PUL + serpy Lig IO * ant 
i oe 


OH 


in which the real canonical coordinate and momentum , Xw, Y,, ate taken to 
designate each normal mode of lattice vibration of the optical Pe taned ass 
lattice constant and v the frequency. Since the vibrations of optical be ie se 
7, nearly the same frequency » it will involve only a small error if we give 


the Reststrahlen frequency » to all waves. 


L 
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The matrix elements of X, and VY responsible for a phonon-absorption and 


-emission may be given by 


1|X,|V,) = gle 1), 
(Ma 1|XelNe) = Fi] 57 (Met) 


(6) 


= h > me 
(NM, £1| Y, |v.) = tea (Moto): 

3. Now we shall evaluate the correction to the electronic energy in crystal 
owing to the mutual interaction with electrons of the lattice vibrations described 
by (1) and (5). Since the first order energy correction is easily shown to 
vanish due to (3) and (6), we have to manipulate the second order correction 
of the electronic energy, i.e., the correction of the “ self-energy” type, which 
shall be found to become of a finite magnitude owing to a finite degree of 
freedom of crystal lattice in contradiction to the case of radiation field (infinite 
degree of freedom). 

Thus, it follows, for the second order energy correction d£, 


Of= yt | (k, N,(w) \dA\k’, NV, (w") |? a 
kiwelt E(k, N,(w)) —E,(k’, Ni(w’)) 


in which £,(k, V,(w)) represents the energy of the system H,4+H, without the 
mutual interaction and & the wave number vector of a crystal electron specifying 
each of the electronic states. 

Now, for brevity, we shall confine ourselves to a simple cubic lattice, which 
model will be sufficient for our present calculation, since our results obtained 


will easily be extended to the more complex lattice structure. In such case, we 
have, as usual, 


EL, =k /2m,*, (8) 

for the energy near the bottom of the upper band, and 
4= —E— Fo (4n-2) 9 
. 2m,* \ a ; @) 


for that near the top of the lower band respectively, in which m,* and m,* 
express the effective masses, 
numbers, 


without tl 


”% the unit vector whose components are integral 
a the lattice constant, and &, the energy gap between the bands 1 and 2 


re energy correction under consideration, the bottom of the upper band 
being taken as an energy oligin. The energy expression (9) may easily be 
derived from the energy formula obtained 


Boor sate © Ty in Bloch’s tightly bound electron 
oximation. 1e eigenfunctions corresponding t 8 
written by P g to (8) and (9) may be 


Ye” (7) =Nexp (kr )u,® (vr), (10) 
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and 
$2 (0) =N-Pexp (kn) u,2 (r), a) 


respectively, in which ~,” (2) and w,°(#) are periodic functions with a crystal 


period and may be supposed to be nearly independent of & in the energy regions 
under consideration. 


Semimetal (Nonpolar crystal) : 


For the non-vanishing matrix elements of 0H, we have, as already shown 
in the theory of electrical conductivity, 


(KN, (0) |S |K—w, Nw) +1) = 9) A (wo) #1) Cn 
pense CL SEP OZ) 
(Ie, N,(w) 8H |ke-+w, N(w)—1)=,/ av (wy ey, 
MN 4nv,(w) 3 


respectively, in which the longitudinal waves (¢=1) of lattice vibrations only. 
are shown to make contributions to the scattering matrix elements approximately. 
For the case of the electron or holes near the top of band 2, the similar expres- 
sion for the scattering matrix to (12) will be valid except for the replacement 
of C, by C,, corresponding to (10) and (11), where the average kinetic energies 
C, and C, are given by 


3¢ 
— | 
© om#ds 


in which r, and s denote the volume of elementary cell in crystal and the dis- 
placement along the direction of polarization of longitudinal wave respectively. 


ee (7) re (i=1,2) (13) 
s 


0 


Tonic crystal (Diatomic polar crystal) : 


: oe eye ane 4 i 1 ae 14 
(ke, V(w) dH |k w, V(w) 4) =somyit tN Od RG Jane (14) 
where C/ represents the constant quantity being nearly equal to onc, and 7=1,2 
specifies each of the bands 1 and 2 respectively. 

Taking into account (12) and (14), the energy correction dE in (7) is 
easily seen to take the following simplified form, 


|(k—w, V,(w) + 1)0A|k, V,(w)) |? 
E(k, N,(w)) —ZB(k—w, Ni(w) +) 
\(k-+w, V,(w) — 1/6 |e, Ni (w)) | , (15) 


sE= 3 | 


w 


0) 


in which the summation over W may be evaluated approximately, as usual, by 


replacing by the following integration, 
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3 (n/a (ala ae | a 
so Ae | do] do,| dw: 


w (27)? —x/a —nla —nla 
aes ["w*deo|"sindtd [ae (16) 
(27)? 0 0 0 


3/ po 
w,=VW 62° /a, 
taking into account the range of values of each component of Ww, namely, 
—1/a < Wy Wy, W; < 7/a. (17) 


Thus, it follows, for the case of semimetal (nonpolar crystal), 


‘Pik (aia) pele deo sinddb ag |— NCO) +" 5 — 1, (1s) 
9MNe (27) “w+7 —awt+y 


Q 


for the bottom of the upper band and 
D2 3 ‘ is 
0B, (ean) = — ACF Ne ( ean sind d0dp | Nie) +2, Mie) _} 
a 9MNc (27)* CU+7 —aw+7 
(18)’ 
for the top of the lower band respectively, where ¢ expresses the sound velocity, 
a,=h?/2m,* (t=1,2) and y=he. 


T> 0: Niwa 2 ST (4’) 
hy rw 
Then, we get 
_ Tha’ ame 
ser WG [Lae Joao) 
\ = On Me \u, : Tew+(Z) og(“ mi dh 
te Vio g men Be fs 7}, (19) 
wy ay ~ Uy, a 
lia? ey Te), 
dE (k= n {2 ee - ( z 
i( “n) = On Me * +—2# +(Z) log Wy 1) 
rm. [emt Ty og eT |e, 7. (19)’ 
Oy dy yy +7 


For the breadth of the forbidden region separating the bands 1 and 2, 
therefore, we obtain 


4E=dE,(k=0) —{—E, + 0E,(kK=n/a-n)}=E,—AT, (20) 


where 


E =F, lka’® Sour £1] ws 
G nee See a Fens (1 log (“204 +1)| 
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‘er lbgr emo AG) en 


_ Vikya® ce C ‘T 20, 5 rae patel + Gl] 20,+ IF log 20 }. (22) 
~ 9x2 Mey ay ay * UW +7 is Cait aeaok, 


The energy correction independent of 7 in (21) is due to the interaction 
with zero point vibrations of lattice, and, further, each term of (21) cannot be 
measured separately but only as a whole (&,) as easily seen, 


Pete: N,(w) = exp(—kew/k,T). (4)” 


Substituting (4)” into (18) and (18)’, we obtain, after a similar process of 
manipulation to the case of 7 > 9, 


37,2 
4E=Eo— “ Ry fmt? {7 *_ exp (- =wiah ge hwyc pS T? ) 


97 MKC 
2m,*c* at ie Inc Mae 2m *C 
2 wo On -m(2a) 0 Bag tn?) 
ke P hols 
*2/4 * 2 * */2 
ia enee cary ae 2m, ie [-2(= 2m, e+e, Te cons fe |} 
he kyl 


reey, 2 {T?—exp (- eS) ca —— 


apts shige amie ao (2 2m,* ae) + z(—™ fe _ tee | 
Re ; 
mee (2 ca | -2( Pie ao) +E, Saat Q2m*c iii 
age 


-8 
(io ea 

_—T 40° ky hie \ (hw ¢ 2 
re fea Cr + ie Ce Pee (-3 2 aah “r+T)), 22)! 

eee ’) p (s (22) 
since the terms invo:ving Z, may be considered to become negligibly small in 
comparison with the remaining terms especially except for very low temperatures, 
as shown numerically. 

For the case of the ionic (diatomic polar) crystal, we have, after performing 
the similar calculation to the above case of semimetal, 


ip te 
4E=E',—f'T, (23) 


where 
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ch { Stee ealP eed a 2 oe V p/h? Wy —1 |, 


Eg =Ly— 7 aap V a,j V ahv > V a,//v-W,+1 
(24) 
and 
Bic: é {[ tan“}( Vu,//v +7) 4 I log Vu,/hv-w,—1 Cr? 
87°a° MY V ahy 2V a,j V u,/hY-wy+ 1 
tan-"(W/a,/ fv -2) Loo e/v-M—1 Jer}. (25) 
+ V Unhv * 2V a/v te V y/hv-twy +14 ~ 
HOTA 
“(AS a,/hy-u V a,/hv-W—1 
oad ae eth {| tan7!(Va,/4¥- wy) vs u log a, / hv + Wp Cc” 
aia Ara®My V a,jv 2Vunv ”° Vu,j/-wt+1t"* 
+| tan71(“a./hv+ wo) 4 ati", log V dy / hv + Wy—1 C a 
V Ugh 2 V aj V d,/ V+ Ws +1 
ma ) 
x exp (— : 26) 
xp ( LT ( 


4. Generally speaking, the combined effect of both self-energy-and radiation- 
damping types will give rise to the observed temperature-shift at the edges of 
the energy bands, but, for the real crystalline insulators and semiconductors, one 
or the other of the mentioned effects may be considered to predominate. In 
order to see such situation and, further, to proceed the comparison with the ex- 
perimental results, we shall carry out the numerical estimates of 4E for the 
cases of NaCl and Si as follows. 


NaCl: a=2.81 x 10- cm, 
M= MyaMa/ (Mya+ Me) =14.0 x proton mass, 
y= h,O/h=0.583 x 10" (0=280°K), 
a,=6.05 x 10~*/e,,  m,*=e,m (m=electron mass), 
a,=6.09 x 10-*/e,, 1m *, == Sat, 
7=3.86 x 107%", 
C/~C/~1. 
4E= E,—0.303(“e, —0.037e,) 


—10.8 x 10-"(We, + ve )(1-0.037 we T (eV). 
& Ey 
(27) 


Since we have at present no reliable data regarding €, and ¢,, the values of 
€;=€,=1 shall be tentatively adopted in the following evaluation. 
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Then, we have 
4E=£,—0.291—2.08 x 107 (eV), (28) 


the temperature coefficient of which is found to become of a comparable order 
of magnitude with that of the temperature-broadening due to the radiative damping 
effect as shown in Radkowsky’s theory and, therefore, the self-energy correction 
here computed should be taken into consideration, together with the level-broad- 
ening one, for the quantitative comparison with the experiments. 


Sk: a=2.7 x 10-*cm 
M=28 proton mass (=46.8 x 10~-*'gr) 


As for, G,, Seitz’s calculated value is 104 (eV)? at 7=295°K, which leads 
to the electron mobility of 243 cm?/volt sec, while, the experimental value is 
300 cm?/volt sec,” which gives rise to C2=84.2 (eV)* according to the well 
known relation 


ple 3 V Qrehc Mn, (% =the density of unit (29) 
VBA) 2 (m*) "Ce cells in the lattice) 


In the numerical evaluation of 0, the mentioned two values for C,’ have been 
taken respectively. Regarding thie value of C2 (average kinetic energy of the 
hole at the top of the lower band), we have made use of the observed value 
f,/P2=3 and the effective mass ratio mentioned below (Pearson and Bardeen) ® 
together with the relation of 


fy Ge La (30) 
fo Ci \m,* 


Pearson and Bardeen’s experiment™ shows @=1 and e,~0.33 or 0.65 according 
to whether we may estimate them by the observed depth of the donor level or 
by the mobility deduced from the changes with temperature of both Hall constant 
and electrical conductivity, but e,~0.65 seems to be a more plausible one. The 
sound velocity in Si has not yet been measured directly as far as we know, and, 
therefore, we have estimated it by the use of one or the other of the following 


relations, 
1 Am kn Cr hol (Rn = W 62? /a, p=density) (31) 
t On hoc" 
ane (Sr wise. V=volume of the crystal Bo 
or CH ) 31V ’ ( ys ye ( ) 


together with the calculated value for ¢ of diamond (Seitz). The results are 
€a=1 : ¢=10800 m/sec. 


from (31) 
6,=0.65 : c= 8720 m/sec. 
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and 
c= 8900 m/sec. from (32). 


Taking into account the uncertainties involved in the estimated values of ¢ 
mentioned above, we have attempted, for the nimerical evaluation of JZ, to take 
the values of c=0.5 x 10° cm/sec, 1.0 x 10° cm/sec and 1.5 x 10° cm/sec respectively. 

Substituting the mentioned numerical values suitable for Si into (20) ~ (22), 
except for c, €,, € C, and C,, we obtain 


Table 1 


r 2-104(er)2 | C2=84.2¢er)2 
e1 C2/C2 | ‘ 2 in Nees ses of 
| x 10®cm/sec a B a B 
penny eat 3 | os 0.113 8.0x10-4 | 0.092 65x 10-4 
ee | 1.0 87 2.1 46 1.7 
iad , 1s | 38 0.9 31 0.75 
0.65 1.02 0.5 0.043 3.4x10-4 | 0.088 28x 10-4 
38. 1.0 24 0.95 19 0.70 
. 1.5 16 0.38 13 0.31 
a : 0.33 0.189 0.5 0.015 1.1x10-4 | 0.012 0.85x 10-4 
earson an 
or 
tei 1.0 | 8 0.27 6 0.22 
15 5 0.12 4 0.10 
4E=L£y\—a—BT ev, 
my*=em; mo*=em (e.=1) 
£-dkE 


(Y) 


cV 
0.5 
&=6=1; C2/C°=3, Cfi=104(ev)2. 
(1) ¢=0.5 x 10® cmjsec ; 
/ (2) c=1.0x 106 cm/sec ; 
0.25 / 


(2) (3) ¢=1.5x 108 cmjsec. 
/ ara 
Fig. 1. 
- isis 35 hearticale (3) 


500 
1000 °K 
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c)) 


4.93 x 106" (14202) (ev), (C, in (ev), ¢ cm/sec) 
‘eae | 


Ale nd 85 x foe 1 2Cr 94x 10", 1% 
€.C; hg 


é 11 


=F,—a—fT (cv), (33) 


and the numerical data of 4E for each of the appropriate values of the quantities 
involved in (33) are shown in Table 1. The general trend of JZ with 
temperatures is described in Fig. 1, including the low temperature region 
according to (26). 

From the detailed analysis of the changes with temperature of the concentra- 
tions of both electrons in the upper band and holes in the lower band, Pearson 
and Bardeen® have observed, for the temperature coefficient B of 4, to become 
of 3x 10-4 ev/degree for Si in the higher temperature range, which value is in 
fair agreement with our theoretical one for the suitable values of the physical 
quantities, i.e., e,~0.65, c~05.x 10° cm/sec, or e,~1, c~0.8 x 10® cm/sec, as seen 
in Table 1. The contribution from the level broadening discussed by Radkowsky 
is easily seen to become negligibly small, for the case of Si, in comparison with 
that from the energy correction mentioned above. 


One of the authors is indebted to Educational Ministry for financial aid. 
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Note added in proof: After sending our manuscript of the present paper to the editor of 
Progress of Theoretical Physics, Fan’s letter in Phys. Rev. 78, 1950, 808 has come to our notice in 
which the temperature-shift of the energy gap in monoatomic semiconductors has been worked on a 


quite similar line as ours and, further, the observed shifts with temperatures have been analyzed 


theoretically. Our criticism of Fan’s theoretical analysis together with our calculated shifs for Ge and 


diamond’ shall be published soon in the forthcoming issue of the journal. 
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§ 1. Introduction 


Investigations concerning the cosmic-ray phenomena underground provide 
impoitant clues to infer the nature of the cosmic-ray mesons of moderately high 
energies. 

The directional intensity and the temperature effect of cosmic-ray under- 
ground are computed numerically in order to give further suppoits to the theory 
of Hayakawa and Tomonaga” which was proposed to explain the well-known 
bend in the intensity vs. depth curve of the cosmic-ray underground on the basis 
of the fact that z-meson decays into “#-meson and the latter, having no strong inter- 
action with matter except an electromagnetic one, can penetrate far underground. 

The result of the computation of the diréctional intensity shows fairly good 
agreement with the experiments of Greisen” and_ thus, though the calculated 
temperature effect is not yet quantitatively in good accord with the experimental 
data of Barnéothy and Forré®”, we may conclude that the rays existing deep in 
the ground consist of s-mesons and they lose their energy according to the 
equation of Hayakawa and Tomonaga” which took only electromagnetic processes 


of #-mesons into account, and consequently that the quantum electrodynamics 
proves its validity up to the high energy of some 10" eV or more. 


§2. The directional intensity 


As was stated in I, most of soft ™-mesons which were produced by the 
primary rays high in the atmosphere decay into #-mesons before reaching the 
sea-level and hard z-mesons that arrive sea-level are immediately absorbed after 
penetrating into the ground, and #-mesons play the leading idle in the ground. 
The directional intensities at various depths may be obtained, following the same 


method as of I, by means of the energy-range relation of #-mesons in the eaith 
and the integral energy spectrum of obliquely incident f-mesons on sea-level. 


Taking the absorption vf a-mesons in the atmosphere into account the diffusion 


th Zand the zenith angle of incidence 


an » I, by the substitution of 7 //cos@. ,: 
results in the transformations of /cos And this 
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and moe (1) 


B- B/cos 0 (1’) 


for the results obtained in I. Explicitly, we have next equation: 


Og(£,4,9) BB /cosf Tee ; 
d//cos 0 are El/cos 8 g(4y 4 0) See 1,0) + ¢(&, 4, 8). (2) 


Here g(Z,,/,0) means the differential energy spectrum in question, and 1/A 
the absorption coefficient of z-mesons for which we assume the same value as 
for nucleons: A=125gcm~*. And 


B= (ptxl/tn) (U/p) =3.4x 10° BeV, 


where /4,, Tx are the mass and the mean life-time of the z-meson respectively 
and p means the average density of air. The numerical value of B is derived 
by putting #,=286m (a: election mriss), t,22 1 x10~ sec. 

By adopting for the source function, go, of m-mesons 


£o( Fort, 9) =const. (Bo°*—2,2E6*")exp(— (Z/cos 0) /A), 
determined by Hayakawa and Nishimura®, Eq. (1) can be solved and we have 
g(Z,, 4, 0) =const. J/cos 6-exp(—//A cos 9) (Li9**° 2.240") 1+ CE SERN CE EA 
(3) 
Then the differential energy-spectrum of /-mesons, f(Z,4, 8), can be written 


down, making some alterations for the integrand’ and the upper limit of. the 
integral” appeared in Eq. (22), I, as follows: 


ete i dE, (' B/cos 6 “go (Ey fl, O)al 


f(L, l, 0)\dE=dE 


E M4./Ha)? Edo Eyl! 
1 (Hey [Bey PE (B;°*—2.285*") i 
Ly ee ee ee | (Zo 0) dk. 
const. d@Z-A- (1—exp(—Z/A cos Dari Tle w +E] (B/cos ay P 


(4) 
As this equation is not correct for the lower energy region because of the 
neglect of all sorts of energy losses of both z- and #-mesons, we fix the value 
of const. in this equation using the data of Rossi” for f(Z, 10, 0°), the spectium 
at sea-level for vertical incidence, at the energy of H=3x10"eV. Accordingly 
FE, 10, 4), becomes 


the integral energy spectrum of #-mesons on the sea-level, 
F(E, 10, 8) =| 7% 10, 0) aE’. (5) 


z 
We integrated (9) numerically for the values of 6=0°, 15°, 30°, 45°, 609, 75°, 
and 90° taking My, the mass of the #-meson, to be 217 m. 

On the other hand the minimum energy E required for #-mesons coming 
from the direction @ to penetrate to the depth 2 gcm™ from the suiface of the 


eaith is given, according to E54. (12), I, by 
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E(x, 0) =—* _[exp{ (p+7r) -x/cos 6-1], (6) 
Gish eg 
which takes account of the energy losses of #-mesons due to ionization, pair 
creation, and radiation expressed by the following parameters: 
a=2.5%10° eV pete em~, 
p=16x10% pergem™”, 
r=10x104" . pergcar. 


For given x and @ we have the directional intensity, /(+, @), by substituting the 
value of & determined by (6) into (9): 


I(x, 0) =F(E(4, 8), 10, @). (7) 
When 7/cos @ is sufficiently small Eq. (6) may be approximated by 
E(#,9)=a-x/cos 0. (6’) 


As the energy is fairly large for + >200mw.e. (water equivalent) the second 
term of the spectrum 2.24>*° has no contribution. This results in the cos'’ @ 
law reflecting the energy spectrum of z-mesons for the angular distribution 
irrespective of whether the 


ee ae the 0 neers pay depth under consideration is 

--—-+-- cos" : 
09 —<--— cos? @ above or below the bend. This 
i 7 (Peagacoge 


is because the effect of the 
transformation (1), appearing 
only in the term exp (—d/ 
Acos@) of (4) which being 
far small as compared with 
unity, is negligible for every 
direction on one hand, and the 


¢ culculated value 


transformation (1’), appearing 
in the denominator of Eq. (4), 
just effects the cancellation of 
the factor cos @ together with 
Eq. (6’) on the other. There- 
fore, if any deviation from the 
cos'* @ law may take place, it 
must be attributed to the fact 
that, at high energies, obliquely 
incident paiticles suffer more 
Solid curves represent downwards the calculated. results energy losses than tkose of 
for x= 200, 730, 1700, 3000 m w.e., respectively, veitical due to radiation and 


i 2 . . . . . 
pair creation, and consequently, we may expect a shaiper angular distributicn 


RELATIVE INTENSITY 
=) 
te 


(0) 10 20 30 40 60 60 70 80 90 
ZENITH ANGLE (DEGREE) 
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at a greater depth. The results of calculations illustrated in the figure show the 
correctness of our expectation. For r=200m w.e. the angular distribution is 
exactly cos'* @ reflecting the energy spectrum of the primary rays; for +=730 m 
w.e. it may be approximited as cos*@ and is a little steeper than the experimental 
result cos’’@ of Barnéthy and Forio®; for +=1700m w.e. it is about cos”*@ 
and this seems to be in fairly good agreement with the Greisen’s experimental 
result of about cos’ @ at the same depth? when we take the experimental error 
into account. Although we have so far only one available experimental result 
of Greisen at 1700m wee. to testify crucially the deviation from the cos'* 6 
law, it my seive, if this expeiiment is really true, es a strongest evidence for 
the validity of our equation (6) for the energy loss. 

Thus it may be concluded that those rays which penetrate far underground 
are #-mesons and that they interact with matter according to the ordinary 
quantum electrodynamics. Fuither the Greisen’s opinion? that the phenomenon 
under consideration can be well exp!ained by the neutrino hypothesis of Barnothy 
and Forré seems highly improbable. 


§ 3. The temperature effect 


The temperature effect of the hard component of the cosmic-ray with regard 
to the atmospheric temperature at the earth’s surface has negative siga and is 
—0.15% per °C, whereas according to the Forré’s experiment which was 
performed at the depth of 1000m w.e. the same effect with regard to the outer- 
air temperature is found to be +0.74% per °C. The positive sign of the 
coefficient under the ground might be regarded as a suppport for 7-p-decay 
because the higher the temperature, the more 7-mesons decay into #-mesons 
before reaching the sea-level due to the upward shift of the height at which 
m-meson production takes place, and consequently the more rays are produced 
which are penetrable into the ground. However, according to the simple calcula- 
tions of Forré” and one of the authors” the underground temperature effect can 
not exceed the value +0.16% per °C. Although this is right as far as the 
qualitative argument is concerned, there seems to be a considerable numerical 
discrepancy between the theory and the experiment. We will examine here 
whether the result which agrees quantitatively well with the experiment can be 
derived from the idea of I. 

Change of the temperature affects the cosmic-ray intensity through the factor 
B. When the temperature rises the air density p decreases, and accordingly 
increases B. The influence of increasing B on the cosmic-ray iostensity can 
easily be seen from Eq. (4). The temperature coefficient has, of course, positive 
siga, since we are neglecting the decay of #-mesons. 

Hereafter we consider only the vertical intensity. At. a point far deeper 
than the bend inequality B < Z, holds and therefore, the denominator of Eq. (4) 
can b2 expanded as 
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Thus we see that the intensity is propoitional to B at the most. The intensity- 


temperature relation may be written down as 
41/[=4B/B=—4p/p=4T/T. (8) 


Accordingly, if the temperature at the production layer of 7-mesons is put for 7, 
the temperature coefficient becomes 


a lara /T~1/250=0.4% per °C. (9) 


Generally, when the above-mentioned approximation breaks down, we can never 
have the temperature coefficient larger than this value. Is there any reason to 
explain the existing small gap between the theory and the experiment? 

First of all the temperature used by Forr6 for the correlation was that of 
the earth’s surface. On the other hand the temperature used in (8) is that at 
the height in the atmosphere where almost all the 7-mesons decay into “-mesons. 
Next, the temperature used in the approximation of Z in which we neglected 
the absorption of z-mesons in the atmesphere meant the average taken for the 
region lower than the height of 7-meson production. But, as regression coefficients 
among these three. kinds of temperatures are found practically unity, the result 
does not depend upon whichever of these temperatures is taken for the correla- 
tion. Consequently, if the temperature coefficient obtained by Forré is really 
significant the cause of the gap can not be hit upon for the time being. 
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Note added in proof: Our Statement, that the angular distribution at great depth can be explained 


by the electromagnetic interaction of #t-mesons accountin 


g for the bremsstrahlung and pair creation process, 
1s not correct, 


The recent experiment has established much steeper distribution, 2==3.140.1 for small @, 
wheras our theory predicts 7 ~ 2.4. This discrepancy will be explained by the n 
#H-mesons with the cross-section about 10-2%m? per nucleon, 
Professor Creisen for his helpful discussion. 


uclear interaction of 
We should like to express our gratitude to 
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§ 1. Introduction 


The covariant form of the S-matrix has been derived by Dyson’ using the 
Tomonaga equation. The elements of this matrix can be calculated, by a con- 
sistent use of perturbation theory, to any desired order in the fine structure cons- 
tant. Feynman? and Dyson’ gave detailed rules for carrying out such calculations, 
and it was shown that the divergencies arising from higher order radiative cor- 
rections could be removed from the S-matrix by a consistent use of the ideas of 
mass and charge renormalization. This formalism was adopted by Karplus ana 
Kroll? in order to calculate the fourth order radiative correction to the magnetic 
moment of electron, and it was verified that the application of perturbation theory 
is 1eally consistent in quantum electrodynamics. 

The success of this renormalization program has given impetus to a re-examina- 
tion of the difficulties in the meson theory. There are three questions to be 
asked in this connection: (1) Can the renormalization program of electrodynamics 
be successfully applied also to the meson theory? (2) In view of the largeness 
of the coupling constant describing the interaction of meson and nuclear fields, 
can one obtain valid solutions to problems concerning with these interactions ? 
(3) Will the cross section for meson scattering and creation remain finite with 
increasing momentum, on account of the reactive and damping influence ? 

If we restrict ourselves to the (pseudo) scalar meson with (pseudo) scalar 
coupling, we can always obtain finite results, to any desired order, by the renor- 
malization program. It has been pointed out, however, by many authors’ that 
the renormalization of mass and mesonic charge is not sufficient to remove all 
divergences in any theory in which a gradient coupling is used. The second 
question above poses formidable analytical difficulties, since we have at present 
no other covariant method than the perturbation one in solving the Schrodinger equa- 
tion. Case’s® treatment of the anomalous magnetic moments of the nucleons, and 
the results obtained by Watson and Lepore® concerning the radiative corrections 
to nuclear forces, suggest the necessity for research about this second question. 

Algo the third question cannot be solved, so long as we apply the perturba- 
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tion method ‘to the S-matrix. In this connection, we remember the Heitler and 
Wilson’s theory of radiation damping’. As far as the meson problem is concerned, 
this theory gives 1easonable results except that the scattering cross section be- 
comes too small for high energies. But because in this theory no account is 
taken of virtual processes, not only the infra-red catastrophe remains unsolved*, 
but also the anomalous magnetic moment of electron cannot be explained. It 
was shown by Pauli? that the Heitler’s foimalism could be generalized so as to 
include all the radiative corrections. If we could introduce the renormalization 
program into this theory, we shall get a reasonable formalism for treating all the 
scattering, emission and absorption processes”. 

The present paper deals with the general theory along this line. As the 
covaiiant formalism of Tomonaga and Schwiager has been used, the relativistic 
invariance of our theory is self-evident. If we apply the perturbation method 
directly to the S-matiix, the most impoitant property, i-e., the unitary character 
of this matiix will be destroyed. On the contrary, this unitarity is strictly main- 
tained in our theory in any approximation. This pcint must be stressed especially 
in the meson theory, because in this case the coupling constant may not be so 
small as ia the electrodynamics. Many app'ications and discussions will be publi- 
shed in the subsequent papers in this journal. 


§ 2. Formalism 


We shall stat from the Tomonaga equation 


ite . 
1 —__—¥|e|—=A(x)®lel, 2.1 
sy PIR HOPE] @.1) 
where A(2) is an invariant Hamiltonian density usually dependent also on the 
surface o. (We use ‘natural units throughout this work; 4=c=1.) The eq. 
(2.1) can be solved by the unitary matrix U[e] ia the form 


PY [o]=U[e]¥[-0o]. (2.2) 
U[eo] is the so-called S-matrix". If the perturbation theory can be applied, S 
is given by 


Ss =>i( ra a i Roa ic. dn FREE Shaves H(%a)), (2.3) 


where the symbol / is the chronological operator defined by Dyson’. As was 
shown by Tomonaga and othe1s, the self-energy matrix must be subtracted from 


the Hamiitonian density in order to carry out the mass renoimalization consis- 
tently. 


If we put 


S=1-ik, (2.4) 
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Ene elements of the mattix R represent the amplitude of the scattered wave. 
Since the S-matrix is not unitary if we content ourselves with the proper approxi- 
mation in the power series (2.3), we shall perform the Cayley transformation 


ofS 
ge feel Nee 
145k 
S is always kept unitary, even if we replace K by its approximate hermitian 


matrix. From the eqs. (2.4) and (2.5) we obtain the Heitler and Pauli’s equa- 
tion of radiation damping 


(2.5) 


HK KR, (2.6) 
which is the fundamental equation we are going to treat. 


The integral equation defining K was established by Schwinger’, assuming 
the first o1der real process does not occur: 


— | aeer(2) V[o}, 
Vio]+ 2)" ax! e[o, 0’) W(x’) Vfo']=1. (2.7) 


If the successive approximation for V’[¢] is convergent, K is given by 


ae |" ae =e Wie eC WG lle ee 


n=1 Biss 
Denoting the z-th order term of (2.3) and (2.8) by S” and AK™, respectively, 
and the anti-hermitian part of S” by S{”, we get the relations 
Reais)? (far 2s. 5). (2.9) 
Since the total energy and momentum are conseived, the matrix clements 
of S, R and XK must have the form 


Ci BAL Prise «stn tel) SO Pe Pa) Poecsno PL Ad Pa sns-sons P02) oe Oh adO) 


P=Sip' and P=3} pf (2.11) 
i=l i=1 


are the energy-momentum four vectors of the initial and final state, respectively, 
and A is a submatrix of A belonging to a fixed value of P=/. Denoting sim- 


ply (~| 4|f) for the matrix elements of A, the fundamental equation (2.5) 
is written down in the more convenient form 


(AIR A) = (IR) 2 FIR IRA) GK Ay Har O-P), 


852 N. Fuxupa and T. Miyazima 


where the integration covers all the final states (including the number of particles) 
2 
f initial state. 
that can be reached from the initia ; : 
We shall often consider a collision between particle L (~)') and 2 (*,)- 
c : . ' 
After the collision, these particles will take the energy-momentum and # res 


pectively, creating a number of new particles (#',.----- ,k™). The invariant differen- 
tial cross section for this process was defined by Méller™. 
doz — 4b | (Poyeereesh™ | R | Por Pre)? OCP sons)» an) 
Q(por Po) 
where 
p= Lp xp) (2x2) (2.12) 


is the flux of incident particles, and 


as oe 2 - a (p)2-+ (APD?) 0-0 (")2+ (ee) "1 2a) 19" (P—P,) 2.13) 


is the so-called statistical factor. aXé6 is a tensor whose components are given by 
(a,b,—46,)- 


§3. The Feynman and Dyson method in the meson theory 


To describe the meson field we choose the simplest possibility which may be 
expected to yield interesting results—namely, the pseudoscalar theory. In the 


case of interacting meson and nuclear fields, the invaiiant Hamiltonian density 
FHT is given by 


H= if, $x) ers PC) be(2) ieee O(z) “aks 7) 
ar Fe 
= SE GR) teks PO) tu) O.G (a) (2) 878 (2) 82), 8.1) 


where the dummy indices 4 and f mean the summation from 1 to 4. 7, zt, 7; 


are the conventional isotopic spin matrices and t, is merely the unit matrix in 
the isotopic spin space. 


The vacuum expectation values are 
<P($,(x), $,(%’)) > = —0, 4,, (x—+’), (3.2) 
<P(a(x), Yala’) >"=—i Saree (4—2') (3.2) 
where the index be refers to both spinor and isotopic spin components, and 4 (4) 
and S»(4) are - times those given by Feynman, except that S, is a unit matrix 


with respect to isotopic spin indices ; namely, 
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Ay (x ask =e ae. ke ’ RS R ees hy ys oS 

= fa ” =a @ ( ) ( ) Ba i.2 ie, ( ) 

Ke: d'pe- a 5 se aie pas 3.4 
Vs = yi | pe (2) Ce i (3.4) 


where € is a positive, small constant which will not be written baat dn in the 
following. As was demonstrated by Matthews”, the third term in the Hamiltonian 
may be omitted if we put effectively 


<P), HEY) > = 2 PEL), ble) >e CS) 
ax, ° 

According to Feynman and eee} the matrix element of S in question is 
given by a sum of term each of which may be described uniquely by a graphically 
represented transition scheme. For any given transition diagram G, the corres- 
ponding matrix element in the momentum space is obtained by the following 
prescriptions. 

(i) The coefficient is 


(=T)s-"-t (— 2) (az) O'(P—-B,), (3.5) 


where 2, m, J, and p indicate the number of vertices, of open polygonal arcs o 
nucleons, of loops, and of internal nucleon lines, respectively. 
(ii) Corresponding to each external line of G, we put a constant spinor or 


pseudoscalar 
(p), BP), G(2) oF G,(2), (3.6) 
whose normalizations are given by 
P (Dp) re (2) =E(P), (3.7) 
$* (PME =S (for each 4), (3.8) 


E(f) being the energy of nucleon. 
(iii) Corresponding to each internal line of G, we put a factor 


Oi, 4 x, (# Ve Srp’). (3.9) 


(iv) Corresponding to each vertex of G, we put a factor 


F=fig, (for pseudoscalar coupling), (3.10) 
Gata, (417 2) (for pseudovector coupling), (3.11) 
lai 


where the + or — sign is used according as the meson annihilates or creates 
at the vertex, #’ being its energy-momentum four vector. When the self-enerey 


term of (3.1) is operating, we have to put 
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iM oor it? (3.12) 


The virtual energy and momentum should of course be conserved at each 
vertex. After the above calculations, we must perform the integrations over some 
energy-momentum variables which may take arbitrary values. If these integrals 
are divergent, however, the result may depend on which energy-momentums we 
shall take as the integral variables. When we want to get the matrix element 
of K™, it is necessary to calculate the anti-hermitian part of S”. It will generally 
be shown that this anti-hermitian part is obtained by taking the principal values 
of the integrals over energy-momentum variables at the poles’. 


§4. The renormalization of mass and mesonic charge 


This section is concerned with the evaluation of modified functions to second 
order of Sy, Dy and F; in the case of pure pseudoscalar coupling. The calcula- 
tions including both the pseudoscalar and pseudovector couplings are so much 
complicated that they will be published in the subsequent papers. After we have 
carried out mass and mesonic charge modification in these modified functions, all 
the Feynman diagrams, which contribute to the matrix element in question con- 
taining second order corrections, will have direct physical significance. 

Since the evaluation is almost the same as that performed by Karplus and 

Kroll’, we shall mention here only the results. 
(a) The function S’,(f), which describes the property of a virtual nucleon as 
modified by its interaction with the meson field, must be calculated to second 
order. The leading term, of course, is S,(~). The second order correction S@( p) 
to this function arises both from the ability of the virtual nucleon to create 
and absorb a meson, and from the nucleon self-energy term (3.1), as described 
by the Feynman diagram Fig. 1, 


AT 


Fig. 1. 


After the mass and mesonic charge have been renormalized, SP(~) becomes. 
Sy(f) which has direct physical significance. The result of calculation shows 


(2M én } 
P l sey Mu + (i 7 ~—M) (1—n) }1———_** 
Sr(P) =q2()I, du w(1—w), .s trp ) d—x) 1 A, } 


A,+ (f° + M*)u(1—u)v 


(4.1) 


where A, is given by 


A= Mt + #i(1—u). (4.2) 
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(b) Quite similarly, the second order correction J¢)(R) to J,y;(R) is described 
by the Feynman diagram Fig. 2. The physically significant part 4»,(4) of 49} (4) 


a Hoy 

ee ee eee iN a ee ~ 

b eid k ie Se 
2 pis ~ 


is given by 


Apt) =— (FE) au u°(1—x)? [ao 


M* i 1 
x eo Wig Cae (4.3 
Ties (2+ p)ul—u)v] Jo Ap+ aoe ae 
where 4; is defined by 


A= M*?— pi u(A—x). (4.4) 


(c) The second order correction Z/(“, 2°) to the vertex matrix F; is described 


by the Feynman diagram Fig. 3. 


Fig. 3. 
L& (¢, #) can be written in the form 
LOE Pye Ar yt Lt, F), (4.5) 


where A, is the divergent integral which may be interpreted as the mesonic charge 
renormalization term. This separation is not unique however, because we have 
no general principle to do so. Watson and Lepore® carried out this separation 
in such a way that 7; .(¢', 2°) becomes) zero: if the nucleons Z, ¢? and the meson / 
should be free simultaneously—which, in fact, is impossible, for £= + (¢’—#°) on 
account of the energy-momentum conservation at the vertices. In this case, Z; 


(A277) “is given by 


1 1 
LA, *)=— (FL) Tare [dee | ao 


¢3 
k 
4n EE \ 40 
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K piel hee 
x | Ae 4 [Mt— we o(1—v) Ts (4.6) 
5) 4liz +*iko 
af fasten. 
i RE I 


where 
A= M? + p2A—u)—4i # v—2) 
+ (B+ upd vl—v) +[(A)24+ e+ (424+) 2), ED 
Mo = + pg A—u)—pe @ v—2), (4.8) 
K=(t7 04+ M)7,0G 7 24+) 1-2) 
+(704M)7,Mut+r(¢ 7 + UM) Mu 
— (B+ pe)tv(1—2)7; 


—[((A)?+ IP) 0+ (2)? +?) 1-2) Jud) 7s (4.9) 
bap ek for i=4 or £=4 
=20,—1 for i4 and £4. (4.10) 


We notice here that in the charged meson theory L(#, 7°) is identically zero. 


$5. An application to the meson scattering 


As well known, the cross section for the meson scattering increases rapidly 
with the-momentum of the incident meson. Many authors have shown that this 
difficulty could be overcomed by taking into account the reaction and damping 
influence of the meson field. It must be poiated out, however, that the calcula- 
tions so far were either non relativistic neglecting entirely the recoil of the nucleon, 
or introduced the approximation which might be considered not always appropriate™. 
We are, therefore, going to treat the influence of damping on the meson scat- 
tering, as an application of the fundamental equation (1). The radiative correc- 
tions and the multiple processes arising from A”, A, etc., will be neglected 
here, which are to be published soon in the following papers. The charged 
meson theory was chosen, because in this case the above integral equations can 
be solved exactly. 

Denoting the initial and final energy-momentum of the nucleon and meson 
by ~, & and g, &, respectively, the matrix element of K for the meson scat- 
tering (not containing the mult.ple process) has the general form 


(Pr BIA | boy ho) = (27)* ECA) GK) PCP) te 1A (P25 Por Fo) (Pr) 


+h (p)t Tk Bp, &} Dos ko) P (po) je (5.1) 
If we put 
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(p, B|R| Por bo) = (27)* Hf (4) Go(2e) 1 P(P) tat XP, 25 Por Fo) (Po) 
+9(P) te VCD, 23 Po 2o) PC po) } (5.2) 
and introduce them in the equation (1), we shall get a pair of integral equations 
XP, 23 Por &y) =ACZ, 23 Por & —_ 7. fae’ Lala 
Bp Po o) (Pp Po o) (47)?W d A(p, & Pik) 
ACP) XPM; por ko), (5.3) 


[a2 B( ree of Th, 
ACP) VP, B 5 Lov Ao) (5.4) 


for the unknown functions X and Y, where we have referred to the center of 
mass system, / is the absolute value of the momentum of the meson or nucleon, 
W is the total energy of both particles, and A(g’) is an annihilation operator 
defined by 


ul 


Y(2,%3 Po &o =B(p, k; Por &o) ——— > 
(2,45 Po ho) =B(A&; P )—Gaw 


A( pt) =—i7 p+mM. (5.5) 
‘It is to be noted here the eqs. (5.3) and (5.4) have nothing to do with each 
other, which is not the case in the symmetrical theory where the neutral meson 


may appear in the final states. 
If we retain only the second order term K® of &, A and B are given by 


Aen tS le 
(Ath) +, See 


B® _ 5+ doit (Po—*) 


(Ah), ee 
where 
a=f?M + 2f¢ (M24 W?) +4 (IM? +3uMW’), (5.8) 
B ie 
a=f?W+ a *W), (5.9) 


a= Ma ELM (HB) + Fae M0) (5.10) 


9 Af, 3 9 9 
b= +E ms gle — (A.—*)’]. (5.11) 

The eq. (5.3) describes the scattering of Y* by WV (neutron) or Y~ by P 
(proton), while the eq. (5.4) does that of Y* by P or Y- by. “Since, the 
latter integral equation depends on the angle of scattering, it cannot be solved 
exactly. Solving the equation with respect to X in the first process, we get 


the result 
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5 hee (5.12) 
1+77 
where 
QIA es a 
= +M 5.15 
ee ) (5.13) 


represents the influence of damping, # being the energy of nucleon. Note that 
in the center of mass system the energy of nucleon or meson (€) does not change 
after the scattering. 

The differential cross section am becomes in this case 


dQ 
Ww ae 


i= 


X'7,), (5.14) 


where we have taken the average about the initial spin and the summation over 
the final spin of the nucleon. Performing the integration over the angle of scat- 
tering and putting 


ou 
Ll % 3s eT 
Fag(se g), Fe & (5.15) 
we obtain the total cross section 
8x C(4) 
——— Fi ——ee 
MM?” x(x)’ CR 


where 
C(x) = (1+ 2pT (a) )F (2) + Fip'a®| —142r47°(—14 p7(2)) —r' pT) 
— (r*p°/16) (0+ 27%(2))*%, (5.17) 
m(x) =| (0+27(4%))(14+297(2)) + Fixe [—142r4+r(—14 pT(2)) 
—r eT (a) — 74e*/16) (p+ 27° (x))*]}? 
+ FYyta"(142pT(2))*| 27 (2) +r(p427(2)) (4p T(a)) (2, G18) 
E(x) = (142°) + 7%(2) —r (p+ 27(2)) { (1+22) 2[2(1 wa 


— (rp/2) (14+-x*)3—r (1+ pa?) 8] 4 7a) 2—r— (rp/2) T(2) | 
i (7°/4) (0+27(2))*}4(2+4 p*x’) +4rpT (x) —2rp?x* 


+ (rp/2)*(L4+2°+4 7?(x)) —2r(4—r) (14+ p7(x))}, (5.19) 
p= 0.15, et (5.20) 
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and #7 (x) is the energy (€) of the incident meson in the laboratory system, 
namely | 


WT r. sd ee 2 2 4 2x? 4 
(7) pena Voces 4 (1+ 2°)3(1 + p?s*)?. (5.22) 


The values of Z7(4) for a number of x are given in the Table I. 


Table 1. 


When we neglect the influence of damping, the cross section becomes 


Aue eg aeie) (5.22) 
M* 7(2) 


where &(%) is given by (9.19) and 


n(x) = (1+2pT(2)) (p+ 27(2))*. (5.23) 


Table II gives the cross section o® for different values of x and r when aol 
For other values of / multiply only the factor F* to this o®. r=0 corresponds 
to the pure pseudoscalar coupling and r=1 to the pure pseudovector coupling. 
Due to the influence of the recoil of nucleon, the cross section does not increase 
even if x reaches as much as 5, or in the laboratory system, the. energy of the 
incident meson amounts to 1.5 times the nucleon mass. This is contrary to our 


usual expectation. 


1 
0.00 x 10-27 6.13 x 10-27 80.0 x 10-27 
2 2.28 0.04 5.94 48.4 
3 1.58 0.23 5.79 47.5 
4 1.13 ‘~ 0.37 5 89 48.5 
| 0.46 : 


Table II. Values of of). 7=1. 


Table III, and II, give the ratio a/o for r=1 and r=2, respectively. 
r=2 corresponds to the case where f and g have the opposite signs. In the 
other two cases, the influence of damping is negligibly smill even atr=5. As 
expected, the bigger the coupling constant and the momentum, the greater the 
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influence of damping.* It is to be noted that tie cross sections o” and o do 


not approach to zero but to a finite value as / diminishes, as was first pointed 
out by Araki”. 


Table III,. The ratio /6@ for r=1. 


Table Ils. The ratio ¢/¢@) for -=2. 
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§ 1. Introduction 


Using the Schwinger-Tomonaga formalism, corrections to the meson current 
due to interaction with the nucleon field have been calculated, and that from 
motives chiefly of the following considerations. 

(1) The calculation of nuclear forces and the anomalous magnetic moment 
of nucleon, applying the Schwinger-Tomonaga formalism to meson theory, has 
been performed by many authors, the results obtained however being not in 
quantitative or even qualitative agreement with expe1iment. There are indeed 
some attempts to get rid of this difficulty, e.g- by modifying the model or by 
advancing the calculation to higher order approximations. But it may also be 
worth trying to attack the problem from a different aspect, namely to deal with 
real meson problems, instead of being concerned with virtual mesons, which has 
been the essential phase of the the theories in question. 

(2) Recently, Kinoshita and Nambu” gave the corrections to the meson 
current due to electromagnetic field, and pointed out that vector or pseudo- 
vector meson should have, besides the anomalous magnetic moment, an intrinsic 
quadrupole moment. If however such a quantity had any real meaning, it would 
have been necessary to take account of the interaction with the nucleon field, 
which is stronger than that with the electromagnetic field. 


§ 2. Method of calculation 


9 


The calculation proceeds according to the method used by Case” in deriving 


the magnetic moment of nucleon. Thus for the requisite effective Hamiltonian 


limited to the second order approximation, we have 


+0 
Hy a= 2 (GE) [fae tee Med, Ha), He), 
!\ he 
where P denotes the chronological ordering opertor, and 
Het = Het HY, 2) 


Hft=—i1e $i, uPA, (3) 
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HeY= Ate (4) 

jai b— 99 3), (scalar or pseudoscalar meson) (5) 
hic Ory OF 

= 2 (twb.—bite) (vector or pseudovector meson) (5’) 

_ 3b, _ by 

Pe i 3 — 

Sy, Os, 


Here.e¢ is the proton’s charge, ¢ and ¢ are the spinor operators of the nucleon 
field, 7, the usual Dirac 7’s, and tp is the isotopic spin operator for proton. A, 
denotes the four-vector potential of the external electromagnetic field, and ¢, ¢ 
and ¢,, $6, describe the charged scalar (pseudoscalar) and vector (pseudovector) 
meson field, respectively. H; tepresents the interaction Hamiltonian between 
nucleon and meson fields, and, e.g. for the scalar coupling of scalar meson, may 


assume the form 
H,=f $ tpy $+complex conjugate, 


where f denotes the coupling constant. In the following sections, we shall use 
J for scalar and pseudoscalar couplings, ¢ for vector and pseudovector couplings, 
and g’ for tensor and pseudotensor couplings. spy is the isotopic spin operator 
changing neutron into proton. 

Now, our problem is to take the zero nucleon, one meson portion of the 
effective Hamiltonian. When limited to the second approximation, it can readily 
be seen that the meson current part of H™, i.e. H¥*, need not be considered, 


‘ 
| vo x 


(a) (b) 
— nucleon line 
ttteee meson line 
—-— photon line (external) 
Fig. 1 


*) The term of order ¢2 is neglected. 
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49 att a rt! 
{ ‘ A x 
1 ‘ ‘ 
' ‘ \ 
P*' ' P i ony 
N- N- 
N 
“ 1% 
: P v Pa 
' ' 
‘ at ' : 
‘ ' ‘ 
’ amt tae 
a (a) 
(a) (b) 
ee Se 
} l 
aE fe p~ 
PF Pt 
je P- 
: 
\ | z° 
(a) (b) 
Fig. 2. Fig. 3. Fig. 4. 


as this part can contribute only to higher approximations than the second. This 
situation makes the calculation simpler than that for nucleon. The Dyson diagram 
of the nucleon current part H7‘, up to the second approximation, is shown in 
Fig. 1, (a) and (b). They differ from each other only in the direction of nuclear 
lines. Detailed processes for positive (z+), negative (77) and neutral (2°) me- 
sons are illustrated in Fig. 2, 3 and 4. In these figures P*, V*, P™ and V7 
denote proton, neutron, antiproton and antineutron, respectively. Circles mean 
the external electromagnetic field. It appears from these figures that positive 
meson processes belong to the (a) diagram, negative ones to the (b) diagram 
and, finally, neutral mesons have processes which belong to both (a) and (b). 
But the (a) and (b) processes have the matrix element possessing equal magni- 
tude and opposite signs, so that the neutral meson has, after all, no correction 
current. For this reason, we did not take the interaction Hamiltonian //; for the 
neutral meson into account. 

| Thus, the portion of zero nucleon, one fosttive meson, in the case of scalar 
_or pseudoscalar meson with scalar or pseudoscalar coupling, comes out to be 
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Ys Me oe A Affe aX, $(4;) (42) Tr] Sp(40—41) ert eSr (42-70) 


8c? 
Fel) Sela) exe 7) |, (6) 
where S, denotes Feynman’s function for nucleon, whose Fourier transform is 
Sox) =2-t me dk exp(kz) al (7) 
nak, (8) 


in which J/ is the nucleon mass. In the double sign of (6), the upper sign 
corresponds to scalar meson, and the lower to pseudoscalar meson. (7;) is unity 
for scalar meson. 


§3. Scalar and pseudoscalar mesons 
As usual, we consider a typical Fourier component of field operators. Let 
d(x) =4(p)exp(ipx), (x) =4(P))exp(—ip’z). (9) 
py and 7,/. satisfy 
p=p"=—*, (10) 


where 
Saar (11) 


# being the z meson mass. Inserting (7) and (9) in the expression (6), noting 
that 


Tr(tp try Txp)=1 
Tr (hy +0) ERT +) (75) ER +) (75) ] 
= Fiky hel ki!" Tr(7a Ta Yo Te) +410 (hat hy £24), 
(Vd Tu Tv To) =4 (Din Dov—Son Day + Ou, Dre )s 
and integrating with respect to (4%, 4%) and (4, 2’), we finally obtain 
Ha ~f A,\ Gp hee’) hu + Re) &,!— CARY RE Ry thy! + te) 
c*(2z7)4 (2 +2x2 ) (42422) (4’" +x5) 
x9(P)8(p)exp(—idp x), (12) 


*) In place of scalar Products £4 ry , 


ku yp and 2,2, we adopted the conventions kx, ky and &% 
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~ with 
Ry! =hy+ (4p) ws ky =hy + py! (13) 
and 
(42) p=Dy'— Pp (14) 


In the same manner as has been done by Feynman, (12) may be integrated 
over &-space, and be expressed as follows: 


ps i. ii 
saa 1).~ \ 


[ Wenn # (2—x) + I—4r4+7)) yA—+) p/%) 
O(x,7)§O(x, y) + (46/m)*v~A—Z) § 


_ydaatsy— 2) (Ap/0)* 4 (yz) —1)in2 ED) + AP/I IAD N(9-4.9/), x 


P (x,y) + (4p/%)*y A—2) O(xy) 
x$(p’)6(~) exp(—742%), (15) 
where 
O(4,y) =14+4(2-y) (4-y—-)), (16) 


mie Ge) (17) 
Returning to the interpretation of 7,, —ip, and —2(4f), as operating on 


the general wave functions, i.e. in —ip,! ee and —i(4p),—> ie 


(d(x) ¢(#)), the expression (15) reduces to 


ef” uf A) 
(fs era ne 


_A(x—y) (4-9-1) F 129-4) + (1—z#+y)} yd—+) CF 


O(x,9){ (x 9-8-2) va x 
yd—-y+4y—-*#) ce P(x, y)—y (A—*) = 
a3 2° 9 + (3(y—*)—1) In VAEAY 
Ons) —90- ee ™ O(2, 9) 
| (18) 


where [[]? denotes the d’Alembertian operator, and 7, is given by (5). 
It is apparent from the expression (18) that scalar and pseudoscalar mesons 
have no additional magnetic moments. The same kolds true for the vector and 


pseudovector couplings. 
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§4. Vector and pseudovectér mesons. Vector and pseudovector couplings 


In this case, the expression corresponding to (6) becomes 


H= 


mn as 4,{ dx dx, $,(4,)$,(4) Tr [s 'y(£p—4;) Tee Sr (41—%) Tew (1s) Tv 


Sy (4%4—42) tre (7s) 7a | (19) 


where (7,) is unity for vector meson. As before, let 
d(x) =$.(p)exp(ipr), $,(%) =9,(P’)exp(—i’z). (20) 
From the conditions 


ag, (2) _ a¢,(2)_o 


Oty Ox, 


it then follows that 


PB (2) =pi’ $.(p’) =0. (21) 


This being taken into account, the replacement 
D’> (4p) pr>— (4p) (22) 
is justified when they operate on ¢,(~’) ¢,(p). 


na all 
4rve & ax ay x 


(tes) 0-2) (“PY 


x|- (PtP) ug ; 
2 (x9) +() y(1—z) 


Me" 24) F (— 2a —a4 yay PY 


Oe, 1, D{OC, n+(“) xa- —*)| 


_ p+’ wa. 
2 


ée Cpe (4A) 9 nv ~e-9) ((2—x) (¢—y) —1) + (12 (1—2z)) 
P(x,9) +(“2-) 91-2) 


Ox.) +(Ppa-a 


+(2) 90-2)" {PtP oes, ,(l+2—y)—((4), 


P(x, ‘7 
Sur— (4P) ,9ys) (2 -+)} 
4 P+P')n(40)V(4P)a 2dy (x—y) (1—z) d,(2’)b.(P) exp (—i4px,), 
2 O(x,y) +(2 +) 2-3) f 
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—_ in the double signs the upper ones are to be taken for vector meson and 
; os ones for pseudovector meson. Due to the operator interpretation of 
P,’s, this reduces to 


=— Opti . z : ; 3 
cf (55,)40 ax jay | AraC) pee ie (o%b,—P 10) 


1 ie a /a¢t ad 
oC ¥ g,— gk Pa 
# he Ore f— $s x CS 


where 
A= (l+z—y) (l—2) 
Ox, 9) -yA—2) 


2 
% 


= {A(z—y)*—4+y) F (—2(4#—y) F 1-4 +9) } 209). ag 


0 (2,9) {O(29)—9 2) = x5 
l+2-y Oy) yt 
enrol (25) 
B __Ae—yMQ—2)@—/)—1F G4—#) =z) +40 +2-9) (l—x) CF 
O(x,y)—¥ —ayo D(x, y)—J (1-2) 9 
1 fae—y)*—24+y) F (— 2-9) F (l—xr+y))}7A—x) OP 
‘ D(x, 9) {99 weed a-»H} x 
ets (2,9) -y1-n 
di Lek) ee A OTT 2 
+( 9 +2 7)! AeSs (26) 
n__ 2y=2) 9) 
VG, 7) J (1—+) mi (27) 


While the second term of (24) gives, according to Pauli”, the anomalous 
magnetic moment of meson, the third term is. due to an intrinsic quadrupole 
moment as pointed out by Kinoshita and Nambu”. 

Thus, the magnetic moments in meson magnetons ¢/2z are: 


1 (4, (%, ey) (2-2) (2) —D 42-2 
\a le jy P(x, 7) 
| (28) 


the 


Vector meson: - =--(£ 
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Pseudovector meson : 

Cy | (ae fa A(a—y) (2—2) (4-9) -D 4% (29) 
Anhc’ mw Jo Jo P(x, y) 
Performing the integration and inserting the value of z meson mass 278m, and 
therefore 4=0.023 (which is very small and hardly contributes to the evaluation), 


Lov et 


we get 


= +021 (<— =—0.11(£). 30 
by = + =), Pye &) (30) 
On the other hand, for the quantity measuring the magnitude of the quadrupole 
moment 
9 1 1 x 
o=(£_) = |arlae 31 
~ on oa Cal (31) 
we obtain the value 
= — ay pk: eS) —30, i » 1 2 
O= 0.0001 ) 1.9 x 10- (5) em. (32) 


It is to be noted that Q has the same value both for yector and pseudovector 
mesons. 


§5. Tensor and pseudotensor couplings 


While in the case of nucleon, the calculation for tensor coupling gave a 


logarithmic divergent expression, a similar though rather complicated procedure 
of calculation yielded a definite result in our case. 


To quote only the main results attained, the magnetic moment and the quan- 
tity Q are given by 


(EE) fe fp alert hens te, 


; P(x, y), =) 
and 
oh fete (cy SB" 32—y +2) (wy) F122 49) £29) 
Antic or O(x, y) 
(34) 
respectively. The evaluation gives ; 
. = +0.91( 2 tl g” 
‘ eer Eve o.11(4-), (35) 
and 
0.=30,, 3207 ees lee E 
a x ) gcse im em? ; (36) 
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§ 6. Conclusion 


Using the method of Case, corrections to the meson current due to interac- 
tion with the nucleon field have been calculated up to the second order for all 
types of mesons and couplings. Definite results for the additional magnetic 
moments of the vector and pseudovector mesons have been obtained. Their values 
are at most several percent of the normal values. There also appeared for the 
case of vector and pseudovector mesons the terms representing an intrinsic qua- 
drupole moment.* For vector and pseudovector couplings, its value is very small and 
the sign is the same for both types of meson. For tensor and pseudotensor 
couplings, it has the magnitude which is no longer negligible, and, moreover, its 
sign is opposite to each other for the two types of meson. These situations might 
have an role, in future, in deciding the type of = meson. 
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§ 1. Introduction 


Recently, various experiments on the artificial meson production were carried 
out at Berkeley”. It seems possible to derive some conclusion on the meson from 
their data and the theoretical calculations performed by many authors”. But 
when we are going to compare both results, we at once encounter many difficul- 
ties, among which the effect of nuclear binding is the most important. We must 
not forget that the experiments are performed by bombarding high energy par- 
ticles on the graphite target in most cases, whereas the calculations are done only 
for elementary processes. Therefore, it is necessary to estimate at first this effect 
of nuclear binding. Here, we calculate it assuming the Fermi gas model for a 
nucleus, and furthermore, other effects will be also taken into account, if they 
are necessary. 

After some preliminary discussions (§ 2), we treat the case of X-ray bom- 
bardment (§3) and of proton bombardment (§ 4), separately. The energy spec- 
trum, the excitation function and the angular distribution of produced mesons are 
calculated. Also we determine the coupling constant of the nucleon-meson interac- 
tion, by comparing theoretical results with experimental ones. Finally we discuss 
the results, which seem to favour the pseudo-scalar meson theory. If this is the 
case, higher order calculations will be necessitated. 


§ 2. Qualitative discussions 


iM is weit known that the Fermi gas model for a nucleus is very useful to 
explain high energy nuclear reactions®. If we adopt this model, we see that 
there are two kinds of effects of nuclear binding on meson production: 1) the 


spread of available energy for meson production and 2) the restriction of possible 
final states. 


Let us discuss there these effects briefly. 


nae 1) The momentum of a nucleon inside the nucleus (target nucleon) 
8 SP, (Py=Fermi momentum, £y=P}/2M= Fermi energy and M=nucleon mass). 
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Thus the spread of available energies of incident 7-rays or peed due to this 


internal motion of a target nucleon amounts tot 


iE 
S2VEE,, where & or & means the initial energy % incident y-rays or protons, 
respectively. In the special case, where £=F, we see that this effect is more 
important for the proton bombardment than for the 7-ray bombardment. 
Effect 2). This effect reduces the number of final states, since the Pauli 
principle requires that the final momenta. of nucleons in laboratory system must 
be larger than Py. The momentum of incident y-ray or proton is &/e or P=E 


eee respectively. P is considerably larger than £/c for A=. Therefore, we 


must suppose that more final states are excluded in the case of y-ray bombard- 
ment than the proton case. 

As has already been pointed out, the effect 1) reduces the threshold energy 
for the meson production considerably”. In general, we may consider that the 
effect 1) increases the cross-section for meson production. While the effect 2) 
forbids some final states and makes the cross-section smaller. 

The above estimation shows that the effect of nuclear binding may be quite 
different for the different kinds of incident particles. In the case of y-ray bom- 
bardment, the effect 2) will predominate and will make the cross-section smaller. 
On the contrary, the effect 1) will be more important in the proton case, and 
will make the cross-section larger. 


3. Meson produciion by x-rays 
i y, ry 


In the experiment of McMillan et al.°, they bombarded the graphite cylinder 
by X-rays with the maximum energy 335 Mev, and detected mesons which were: 
emitted nearly at the right angle to the incident beam. For analysis of their 
data, we must take into account following corrections besides the above effects 
of nuclear binding. 

3) Energy spectrum of X-rays 

4) Geometry of expetimental arrangements 

5) Scattering from the absorber 

6) Absorption in the carbon sylinder. 

We distinguish the laboratory system (the target nucleus is at rest) from 
the nucleon system, in which a nucleon is initially at rest As is easily seen, 
the changes in momenta of a meson and a y-ray, arising from the transformation 
from the nucleon system to the laboratory system are negligible, while that for 
nucleon is of a significant order of quantity. This allows us to neglect the effect 
1) and to calculate the energy and the intensity of mesons in the nucleon system. 

Then the effect 2), ie. the exclusion of final states, can be easily estimated 
as follows. First we determine the recoil momentum of a nucleon in the nucleon 
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—<—<———— ————— > 


incident y-ray 


Fig. 1. PO is the recoil momentum in the nucleon system. Let P be the 
origin of the momentum space in the laboratory system. Then the initial 
momentum of nucleon is uniformly distributed in the sphere 4, and the final 
recoil momentum in &. Pauli’s principle excludes the final states corresponding 
to the intersection of 4 and &, which is shown as the hatched part. 


No. of mesons 
in arbitrary scale 


+ experimental data of McMillan etal. 


5 


ve 


0 50 100 150 MeV 
: ene 
Fsi. 2. Energy spectrum of mesons mir abo 
a) Scalar meson with scalar coupling. 
Curve 1: We consider only 2) and 3). 
Curve IL: Furthermore, effect 1) and 4) are added. 


Curve mM This final curve is obtained from IL with corrections 5) and 6). 
Theoretical and experimental curves are normalized to fit both best 
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No. of mesons 
in arbitrary scale 


10 


0 e0 1? 150 MeV. 
energy of meson 


Fig. 2. b) Pseudoscalar meson with pseudoscalar coupling. 


system, which is a function of the incident momentum of y-ray and also of the 
angle or energy of the emitted meson. By the transformation into laboratory 
system, this recoil momentum is distributed uniformly in the sphere B (Fig. 1). 


Then we can find the forbidden part of final states as the intersection of the two 
spheres A and B. 


Energy spectrum of produced mesons 


scalar meson with scalar coupling, S(s), and the pseudoscalar meson with pseu- 
doscalar coupling, Ps(~s). The cross-sections of meson production for these 
cases have been calculated by various authors”. 


We consider the case of the 


The curve I in Fig. 2 shows the energy spectrum of mesons emitted at the 
right angle to the incident beam, where only the effect 2) and 3) are taken 
into account. As seen at once, this curve is nearly flat and drops sharply at~ 
110 Mev, whereas the experimental one has a long tail in higher energy regions 
This shows that further corrections are necessary. 
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Heretofore, we have approximated the energy of mesons as uniquely deter- 
mined by the incident energy of y-rays. The effect 1) and 4) make the energy 
of mesons distribute around the above unique value. We have calculated this 
energy distribution of mesons for a fixed energy of 7-rays and given it in the 
appendix. The half-width of this distiibution is an increasing function of the 
incident energy of y-rays. Therefore, it may be sufficient to take it into account 
only in the region of higher energy, rather at the upper end of the cnergy spec- 
trum, where it drops sharply (curve I in Fig. 2). The curve II in Fig. 2 is 


obtained in this way. 


Fuithermore, the corrections 5) and 6), i.e. scattering out and absorption, 
are estimated. For the former, we apply the formula of multiple scattering”, 
which reduces the intensity of mesons in high energy region. For the latter, 
we take into account only the ionization loss®, which is effective in lower energy 
region. 

The curve III in Fig. 2 includes all corrections estimated above. We find 
that there is no appreciable difference between the cases S(s) and Ps(fs), and 
both are nearly in agreement with the experiments. 


Results obtained by the calculation on elementary process 


———— Corrected by the nuclear binding “effect 2) 


cross-section 
in arbitrary scale 


-1 0 +e 
cos § (angle in laboratory system) cos 6 
i) Scalar meson with scalar coupling ii) Pseudoscalar meson with pseudoscalar 
coupling 


Fig. 8. Angular distribution of mesons produced by y-ray with energy 230 Mev. 


; nia) ar distribution McMillan et al. found experimentally that the angular 
distribution of mesons is nearly isotropic, if the geometry of their experi- 
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mental arrangement, i.e. the effect 4), is taken into account. The experimental 
data on mesons with energies 43~59 Mev. are shown in their paper”. 

We calculate the angular distribution of mesons produced by monochromatic 
y-rays with energy 230 Mev. Our results may be compared with their experi- 
mental data, since this y-ray emits a meson with energy~50 Mev. at the right 
angle to its direction. 

The full curve in Fig. 3 is the angular distribution obtained by the calcula- 
tion on the elementary process”. In both cases S(s) and /’s(fs), the pronounced 
maximum is seen in the forward direction. 

The effect 2) considerably reduces this maxima, since the meson emitted in 
the forward direction is accompanied by the small recoil nucleons and large part 
of the final states are excluded. We represent in Fig. 4 this correction factor 
for y-ray with various energies. Applying this factor to the above full cutive in 
Fig. 3, we get the corrected angular distribution for 230 Mev. ;y-rays. This is 
shown by the broken curve in Fig. 3. 


Correction 
factor £=326 Mey, 
1 
280, 
230 
Fig. 4. Correction factor due to 
186: eee 
nuclear binding on the angular distribu- 
tion of mesons prodiced by monochro- 
0.5 matic gamma-rays. 
0 
= 0 ; rt 


‘cos 0 (aagle in laboratry system) 


We may conclude that the angular 
In the case S(s) the strong 
ems to contradict with 


Further corrections are not important. 
distribution is nearly isotropic for the case Ps(-ps). 
maximum in the forward direction still remains and it se 
the experiment. 

Absolute yield Finally, we compare the theoretical and experimental 
‘tted at the right angle to the incident beam, 


absolute intensity of mesons em 
d one nucleon for meson production 


The mean cross-section per one quantum an é ees 
) i roi D i ined from the experiments as 
by X-rays with energies 140~335 Mev. is determine p 
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follows, 
3.8 x 10-” cm*/sterad. 


The scattering out and the absorption reduce the meson yield approximately ay 
90%. Thus the corrected value is 4,8 x 10-® cm?/sterad. While the theoretical 
values are 


OF ¥55x 10-3 cm*/sterad, for S(s), 
hie «he 


and 


£ Ff x32x10- cm?/sterad, for Ps( ps). 

hc he 
The effect 2) makes the cross-section nearly 15% smaller than that in free col- 
lisions. Comparing these figures each other, we can determine the coupling cons- 
tant f of the nucleon-meson interaction as follows : 


f2/hc=1.15 for S(s), 
and f° /kc=16 for Ps(ps). 


§4. Meson production by protons 


Next we consider the case of proton bombardment, in which we need only 
to take account of the effects 1) and 2) described in §2. Here we use the 
center of mass system of the incident proton and one of the target nucleons, 
because theoretical cross-sections are conveniently expressed in this system. 


First we treat the effect 1), for it is more important than the effect 2), which 
will be discussed later. 


Energy spectrum of produced mesons” For 
simplicity, we assume that the momentum FE of a 
nucleon in a target nucleus is distributed uniformly 
on the surface of the occupied sphere with radius 
P,(~+200 Mev./c) in momentum space (see Fig. 5). 
The energy spectrum of produced mesons in the 
center of mass system have been given by Fukuda 
and Takeda®. Transforming their result into the 
laboratory system (the target nucleus is at rest), 


multiplying the suitable weight factor, which is pro- 


ae and finally summing up over 


Fig. 5. The momentum p of a 


nucleon in a target nucleus is poitional to 1— 
assumed to be distributed uni- of’ 


formly on the surface of the all possible orientations of P, we get the spectrum 
sphere with radius py. which can be directly compared with the observed 
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data? (cf Fig. 5). The full curve in Fig. 6 shows the energy spectrum of mesons 
produced by 345 Mev. protons for the case of Ps(fs). On the other hand, if 
we do not take into account this effect, we get the result expressed by the dotted 
curve in Fig. 6, which is not in agreement with the observed results. Further- 
more, this effect leads to a similar spectrum of mesons even if we assume S(s). 
Thus we see that the main feature of the observed energy spectrum can be ex- 
plained by the effect 1) and the assumed type of mesons is not relevant. 


——— Assuming initial nucleon in a nucteous is at rest 


—— Correcting above curve with the effect of internal motion 


$ Observed value in‘ appropriate scale 


relative yield in arbitrary: scale 


140 160 180 200 220 240 20) 280 300 320 340 MeV. 
Meson energy (including rest energy=144 Mev.) 


Fig. 6. Energy spectrum of mesons (r+ and x) (345 Mev. proton bombardment) 
Absolute yield” The total cross-section per one target nucleon for 345 


Mev. protons can be obtained by integrating the above result over the meson 
energy. 


o= sha se 0) 0318 u—O 00113 B ’ Ao 


(2) Assuming the pseudo-scalar meson theory with pseudo-scalar coupling and 
using the phenomenological nuclear force” (quoted as Ps(ps, NV. F.)), 


o=L(*) x 9.5ux10~ red. 
hc \Mc 

Where a and are listed in Table Bi 

y means the ratio of the number of positive mesons to that of negative mesons, 

and the observed value” is r~3. In the above calculations numbers of prctons 
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Table I. 
a 
ji | E 
a es 
mt pn e er 
punt 2 / 4 
a 1 1 
{npn 
pune 2 | 0 
nut 1 | 1 
° = 
mn—>fnn— 2 2 
nnn? 0 | 0 


In this table, /, 7, z°, z+, and z— mean proton, neutron, and neutral, positive, and negative 
meson, respectively. 
and neutrons in the target nucleus have been assumed to be equal. For com- 
parison, we quote here the result for the case where the initial target nucleon 
is at rest and the nuclear binding is not taken into account. 


(3) For Ps(s) ; a=(7) (4) (0.0102a—0,0044 f), rx6. 


Comparing (1) with (8), we see again that the ettect of nuclear binding 1) is 
quite important in the case of proton bombardment. 

Equating these results to the observed total cross-section per one nucleon 
2.7 x 10-* cm’, we can determine the coupling constant as follows ; 


f° /tc=2.0 for Ps (ps), 

I? /aic=270 for Ps (ps, NF), 
and 

f?/&c=0.35 for S(s). 


The results of Ps(~s) and S(s) seem reasonable, while Ps,(ps, NF.) leads to 
quite an unreasonable result. This may be due to the special character of the 
phenomenological force used here. Therefore, we may suppose that this assumed 
phenomenological nuclear force is not adequate in higher energy region. 

Furthermore, if we assume that the momentum distribution of target nucleons 
is uniform inside the occupied sphere, the agreement of the calculated and the 
observed energy spectrum must be better than in Fig. 6. 

In order to examine the effect 2), i.e. exclusion of final states, we consider 
the two typical cases (a) and (b), as given in Fig. 7. In the rormer case (a), 
the final nucleons have relatively large moment» because the resultant momentum 
of the two initial nucleons is large. And the effect 2) is rather small. On the 
contrary, the final nucleons in the latter case (b), may have smaller momenta 
and thus considerable effect is expected. However the contribution of the case 
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(b) to the total cross-section can not be so large. Thus we may neglect the 
effect 2) in our calculations. Strictly speaking, the consideration of ‘hte effect 2) 
results in the decrease of the number of mesons with higher energies, since these 
mesons are produced mainly in the case (b). This also makes the agreement of 
the calculated and experimental energy distribution better than in Fig. 6. 

Excitation function for slow meson production” Let o(Z,¢) denote the 
cross-section for production of meson with energy € by an incident proton with 
energy £ in the nucleon system defined in §3. Then, taking into account the 
effect 1), the yield of a meson with energy ¢ is approximately given by the 
following expression : 


Pr 
‘“ P. 2 : oe alk 
a(Z£, =| tPr a(t 3 ee se dP’, 
F 
2 


where P means the momentum of incident proton. Fig. 8 shows the relative value 
of o(Z, e) for e=6 Mev. in the case Ps(ps), while the observed results represent 
the relative yield of z~-mesons with energies 2—10 Mev.” The agreement of 
both results seems satisfactory. This excitation function depends mainly on the 


available volume in momentum space, for which a is larger than the 


M 
threshold energy. And the assumed type of mesons is not essential. 


—_—__ SEE anni ad 
Py Po Fig. 7 
(a) 


(b) 


§ 5. Conclusions 100 
From the above analysis, we see that 2 
the effect of nuclear binding is quite dif- = _ 
ferent for different kinds of bombarding i 
particles. Furthermore, even if we con- E ips Sheet eet. 
sider the mesons produced by one kind a with 2-10 MeV 


of particles, this plays a significant role 
for some of their behaviour, while only 
a little for the other behaviours. Thus 
it is desirable to compare such experi- 200 = 250 =~ 300 350 Ko 
mental arid theoretical results as are not Fig. 8 

seriously affected by the nuclear binding, in order to determine the type and the 


—- Calculated relative 
yield 
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coupling constant of the mesonic interaction. 

The above analysis are summarized as follows: 

1) We can say only little about the type of mesons from their energy spec- 
trum, becavse the spectrum is largely modified by the nuclear binding from its 
original form, which is obtained through the calculation on the elementary process. 
The modified spectrum is nearly in agreement with the experimental one, whether 
we assume S(s) or Ps(ps). 

2) The angular distribution of mesons produced by X-rays are of special 
interest. The experimental data seem to favour the Ps(fs). 

3) The excitation function for mesons produced by the proton beam depends 
only a little on the assumed type of mesons. The experimental curve can be 
explained for the most part by the available volume in momentum space, and 
the cross-section plays only a little role there. 

4) The ratio of positive to negative mesons with high energies is little 
affected by the nuclear binding. But this ratio depends so slightly on the type 
of mesons that it is hard to deduce any conclusion at all. More precise experi- 
ments are hoped on this point. 

5) At present the absolute yield of mesons are the most important data. 
We can determine the coupling constants from them. Our results are summarized 
in Table II. 


Table I. 
The values of coupling constants 


Determined from 


S(s) 
Ps (gv) 
Ps (75) 
Ps(¢s, N.F) 


meson production by X-rays meson production by protons 


As is seen from this table, Ps(~s) gives inconsistent values for the coupling 
constant, while S(s) and Ps(pv) lead to somewhat satisfactory ones. But it 
may be possible that the 4-th order calculation revises this inconsistency in the 
case of Ps(p~s)”. Furthermore, the phenomenological nuclear forces fitted to the 
lower energy phenomena do not seem adequate in higher energy regions. The 
same conclusion was also derived from analysis of cosmic-ray stars. In connec- 


tion with cgpture processes of negative mesons, further discussion will be given 
in another paper. 
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Appendix 


Distribution in meson energy due to effects 1) and 4). Let k, and P, be 
the incident momentum of y-ray and nucleon and g and P be the final momen- 
tum of meson and nucleon, in the laboratory system. We may consider P, only 
in the direction of K,, for all are symmetric with respect to this direction. The 
kinetic energy Z, of meson is obtained from the conservation of energy and 
momentum. 


k ky +P, i] 
By=E,— fh, Cot aE a cost (1) 


where @ is the angle between K, and p, and &, is the value of £¢ for i=5, 
PaO” (P= (le, 24) / ho) - 
The distribution function of A is given by 
3 P 
ralarad Gita yeiwet), (2) 
_ where 
y=Po/ Pr. 


From the experimental arrangement, the detection probability of mesons emitted 
with an angle @ is approximately given by, 


S(-a)dx (l>x2>-1), @) 


_ where 
zx=cos 8. 


E, can be expressed by new variables x and y as follows, 


E,=E,+ Ax—By+ Cxy, | (4) 


A= hod Ba lobe. and ca Pt. 
M M 
Fae cae 7 
g in B and C may be approximated by VEP+2E yp, ie. the value at OFese 


P,=0. Then we have following numerical values of 4, B and C for 4=340 Mev 


A="1 Mev., B=75 Mev. and C=45 Mev. 


We may neglect the term Cxy in (4), as we are interested in the case for which 
ay <1, and as C is somewhat smaller than A and B. Then it becomes simply 


E,=E,+a(x-y), (5) 
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aw75 Mev. 


The probability F(Z) for Z, to be larger than Z is obtained by the integral, 


PE) = NI de dy -2) (1-9). (6) 


The integration is carried out through the region 5), as shown in Fig. 9. Neglec- 


. L— EF, : 
ting the higher order terms in °, we obtain 
a 


— E— fs 
FE) = 3+ = 2—*140(=—*) . (7) 
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§1. Introduction and Summary 


Yukawa’s meson theory received brilliant support in the recent experiments 
concerning mesons artificially produced by the Berkley cyclotron. 

It is, however, well known that in spite of the advance of the quantum 
electrodynamics and qualitative agreements between theoretical results and ex- 
perimental data, the problem of nuclear forces and. nuclear magnetic moments 
have not been solved completely. These facts show that the present model of z 
meson is inadequate to describe those phenomena. Now let us consider the 
history of the meson types. After the prediction of Yukawa, Kemmer” found 
‘that there exist four different kinds of meson fields. The development of the 
meson theory based on his classification seemed to be successful ‘in the qualitative 
interpretation of almost all phenomena. 

Especially, the pseudoscalar meson field has been studied by several authors. 
- Further, the mixed field theories came out in order to remove the I /?* singularity 
by compensation. However, these theories depend on the Kemmer’s assumption. 
By the appearance of Tomonaga-Schwinger theory, which is based on the idea of 
renormalization of mass and charge, ambiguities which remained in the subtraction 
procedure are removed. In this way, as Case showed, it turned out to be possible 
to obtain the finite results for the nuclear magnetic moments in the scalar and 
pseudoscalar fields which were not in accord with experiences. The vector and 
pseudovector fields give divergent results owing to their tensor couplings.” 

Under these circumstances, in order to remedy the difficulties, we investigated” 
in the previous note the genetalization of the interaction Lagrangians suggested 
by Schoenberg,” and obtained nuclear potentials for the spin 0 field, though 
no explanation for the introductioa of those interactions was given. Further, 
magnetic moments of nucleons by the spin 0 field were calculated according to 


Feynman-Dyson’s formalisum (unpublished ). 
It was shown that the contributions to the magnetic moments by generalizing 


* 1. Enatsu, Prog. Theor. Phys. 5 (1950), 102. This paper will be referred to as I. 
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interactions are of small! influence on the results. That is, the usual terms have 


the form 
(const.) (¢*p,c HY), 


on the other hand, the contribution from the current which comes from the 
process that a nucleon. emits a meson by the scalar coupling and reabsorbs it by 
the pseudoscalar coupling, is given as follows: 


(const.) (Y*p.0 H¢). 


It is clear that the latter is not large compared with the former due to its operator 
f2. Such an attempt carried out hitherto, has led to the conclusion that we 
could not expect the quantitative agreement for the spin O field. 

At the present stage of our knowledge it is of interest to investigate further 
the spin 1 field theory. Contrary to the spin O field, it has the analogy with 
the electromagnetic field. It is the aim of this paper to draw attention to a classical 
aspect of the spin 1 field, and to show more insight into the meaning of the 
generalization of interaction. It will be shown that our basic ideas are expressed 
as follows, 


Maxwell equations — Extended Maxwell equations by Dirac, 
t 


Vector 


+ * 
Peudobactor| He'd equations — Spin 1 meson field equations. 


In §2 we shall briefly survey the extended Maxwell equations by Dirac.” 
He put forwaid the question with respect to the dissymmetry of the Maxwell 
_ €quations between the electric charge and magnetic pole, and showed that, although 
there is the symmetry between them, the great difference between the numerical 
values of them explains why electric charges are easily produced and not magnetic 
poles. We shall examine the formal extension of this theory for the classical 
vector meson field, and find the same relation as before. In §3, taking account 
of the dual charactor of fields, Kemmer’s vector and pseudovector fields are re- 
examined, and it is found that there is a hope that a new point of view will 
lead to an idea for modification of Kemmer’s theory. Nuclear forces are derived 
for the spin’ 1 field according to a canonical transformation method in §4. Then 
concluding remarks on the results are given in the final section. 


§2. Theory of Maxwell and its extension by Dirac 


Now, briefly, we shall sketch Dirac’s theory. His relativistic notations are 
used throughout this section. (See reference (5)). 


The ordinary Maxwell equations are 
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OF, 
eave A 
et = — Ar (1) 
O(F"). 
¥ =(0) ; 2 
ae, (2) 
where fuv= —Fig 1 a= and (F'),,=—F,. 
Dirac’s extensions are 
| AF, ; f dz 
kv 4. =—4 ue ate 
se 7 jy=—4n dye |( : ) 3,(4—2)ds, (3) 
Fi 2 
BAD ws! ang —An ae |(4 #) a(x—2)ds. (4) 
so Bo 9 ds 


The last equation asserts that the modified theory allows magnetic poles. 
It is, however, remarkable that equation (4) is reduced to the following form 


Ofre tear ae faK 2 : 
+ MY = —Athyruy, B) 
Ox” Ax* a Ox* ea ©) 


namely, in the usual notation (not the dual one, /,) equation (5) shows that. 
the current of magnetic poles is a pseudovector. 
The usual introduction of electromagnetic potentials, 


Fy = 2 — 24n 
Ore Ox’ 


(6) 


is no longer possible when there are magnetic poles, therefore, equation (6) must 
be replaced by an equation of the form 


ad, aA, 


we t 

Fu HAE, (ew, @ 
where 
Oy a, Oy OVy i 

G =s{{( » Wy _ On Be ayy aetde 8 
py—& Lema e aon Or. ) (a—y) dey de, (8) 
OGuv (4 SPL GZ (9 
Ax, Be ds ) aC ) Sy ) 


and y, is a general point on the two-dimensional sheet in space-time which is 
generated by Dirac’s stiing. 

We recall here the expression of the tensor coupling in the vector meson 
theory. The last term of equation (7) represents the source functions which 
correspond to the electric and magnetic dipoles in the meson theory. 

The retarded potentials due to a charged particle and a pole are found to be 
given by the expressions 


At Gre | 72-2) (=) ds, (10) 
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and 


ay’ ay’ A(4-Y) ye dz (11) 
Ay (4) =£ &v -\| ect 


where 
J(4)=20(4yx") for % >0 (12) 


=0 for Xo a 0. 


and €y,po is the antisymmetiic tensor of fourth rank with €p;=1 as defined in 
(I). It will be seen that potentials and dual fields satisfy the following relations, 


(F "8 = eves OAS —4rG%3 ‘ (13) 
ox* 
, ¢ 
Dak (2) =40 De [(E) (ea) dot ded BAe tes (14) 
e€ ds 9 fe) v 


Finally, Dirac’s conclusion with respect to the coupling constants is as 
follows: All charges should be integral multiples of a unit charge ¢ connected 
with the pole strength g by the relation 


ethos (15) 


As is well known, the fine structure constant is 


c 1 
=, 16 
he pe ASE (19) 
then, we obtain the value of ¢, 
o 137 
co ee ee 17 
tie 4 - Ga 


ia The great difference between ¢ and g explains why electric charges are produced 
and not magnetic poles,” he stated. This is the outline of the extended Maxwell 
theory by Dirac. 

We shall now return to the meson theory, and attempt to test the Dirac’s 


We may summarize here the main points which are necessary to the later 
discussions. 


idea. 


(1). The formal extension of Dirac’s theory for the neutral vector meson 
ficld (Proca’s field), is possible, and between a mesic charge f,; and a mesic pole 
fy there exists the similar relation that holds in his theory. In our cases a term 
A, (x) must be added to the left-hand side- of equations. (1), (8) and (14), and 


e and g are replaced by f, and f, respectively. Further, retarded potentials and 


retarded fields are modified in order to include the meson mass effect. However, 


their explicit expressions are not necessary to find the relation of f, and fy in 
our rough discussio 1s, therefore for simplicity we shall discard them entirely. Thus 
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we have for a small meson mass (x~0), 


Pag ~ Ste, (18) 


or 


ee ee) mii 1 
Ww hc 4° ee 
The numerical values are appreciated, i.e 


Z-)~05, 1 AOR 5 2 oe 
Cc 
fr) ~ 05, O25 0.125 phy, 0,083p oF nnn 
nC 


It is of interest to see that, as numerical differences between /, and /, are 
not so large, contrary to the electromagnetic ones, we could not exclude either 
of them from the beginning. But, we should not forget the following circumstances. 
That is to say, Dirac’s theory comes out from the simple model : at first, a charge 
generates an electromagnetic field around it, and there is a pole beside it, then 
the pole interacts with the charge through the electromagnetic field, or it a1ises 
in the reverse order. Hence, when an elementary particle have both a charge 
and a pole (or a mesic charge and a mesic pole), the validity of equation (15) 
(or (18)) is uncertain for the difficuities of the interaction of a particle with its 
own field, as Disac stated. ‘Therefore, we may regard the equation (18) as the 
provisional formula in the following discussion. 

(2). Fields and potentials which are dual to /,, and A, are defined by 
the formulas 


I Api 
(F) n= 5 enol! po (21) 
and 
Ca er eer. ie (22) 


We can easily examine this statement in the equations (2) and (13). These 
formulas are very suggestive if a generalization to the meson field is intended. 

In fact, we shall show later how these formulas can be applied to introduce 
the generalized interaction Lagrangians. 

(3). If we modify the Maxwell equations to allow the magnetic pole, 
source functions must be added to the wave equations for the usual potentials. 
It can further be seen in the equation (11) or (14) that the added terms are 
vectors which are considered as the product of €,,,. and a trivector. 
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§3. Field equations of the spin 1 meson 


We shall now attempt to verify the idea suggested above for the spin 1 
meson field. In accord with the previous note (1), we shall use hereafter Kemme’s 
notation, namely Greck and Latin indices run from 0 to 3 and I to 3 respectively, 
and 4)=cl, = 4%, HI). =F; - Sam —ls £u=Sam=fLsn=l. 

According to Kemmer, since there is no specification of the transformation 
of wave functions for reflexions in the spinor notation, two dual field quantities 


correspond to a spinor. For the spin 1 mseon field, they are 


(Vector) 


0g, a9 
t= e3» 23 
Ot, O%s sf @) 
| ae” ag 
=x9 9 94 
and OX. vi: 
(Pseudovector) 
005+ <s Ober bas = . 
Or, OX 3 Or, hn a. 
oy gt 
aR ot. (26) 
As in the formulas (21) and (22), we can postulate the following relations, 
as : l «3s 7 
¢é caked cnr, , — (FE, H), (27) 
grees, (A, $), 28) 
where $=¢, (6,$%3)=—A, (Lor Xeer Yor) =H, (Yew Xu Ae =—H, 
and € 0123 i. 


Then it will be seen that e ions (23) ; oe a 4 . 
‘ at equations (23) and (24) coincide with 
(26) and (25). (24) coincide with equations 


To show analogy with the electromagnetic field all quantities describing the 
meson field are denoted by the same letters, then they are 


‘ 1 dF 

url A —— A= ty 35 = 

cur ar +xA=0, div E+x$=0, (29) 
hs i} 0A 

E= a we sheds : H=curl A. (30) 


——— Patino the assumptions (27) and (28) we may say that the 
a es ait in spin 1 field into two kinds of them is not necessary and the 
ea Ae sates field equations are different representations of equations 
gis ages 4 — other words, merely the reflexion charactor is represented by 

ymmietrical tensor of fourth rank €%"°. This corresponds to the well- 


— 
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known fact that Maxwell equations are only expressed by the usual field tensors 
Fy, or by the dual ones Fy. 

Now we shall study the Dirac’s idea for the spin 1 -meson field by a simple 
example. From the equation (23) the next one is derived : 


OYas Ora OXsy 
Sb ae wee 31 
Ox, Os O04. Gh 


or in vector notations, 


bs ie +curl H=0, and div H=0. (32) 
F 
Corresponding to the currents ot poles in the electromagnetic field, we are 
able to introduce source functions into the right-hand side of (31) or (82). 
For example, if we assume the from (¢=1) 
0H 


aya +curl H=¢y*ap— ae [( grad $*) Bb—$*8( gradf)], 


(33) 
: ry | Ogh* Oy 
div H=—g $2 | 9 py_yrp |, 
ov 5 aye Mids 
then, using the Dirac wave equation for the free nucleon, Proca equations have 


the expressions 


B=— 9A _ gradg+-& g*py, 
Ot 2Qx 
: (34) 
=curl A+ b*po¢, 
2x 
and , 
curl H— OE + A = a [2 (b*p,0¢) + curl (sp.54) | 5 
at ig Wa J 
o (35) 
div E+xf{=— 5 div ($*p.8f) » 
or eliminating field quantities, 
pA—2A=— 8 | 9 (pr poy) + curl (Y*p.6 | 
2x L Or (36) 


O¢—#6= ae div ($*p,.G¥) . 
Xx 


Equations (33), (85) and (36) correspond to Dirac’s ones (4), (3) and (14) 
respectively. Hence, we may conclude that if the second unusual terms of the 
right-hand sides of (33) are disregarded, the statement (3) of $2 holds also in the 
meson field. Namely, as simply verified by reverting to the tensor notation, the 
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(34) and (35) are a trivector (vector multiplied by 


source functions of (33), 
e3T ys 


e*7®) a pseudotensor and a vector (trivector multip ied by € 
Apart from above" discussions, we shall now set up the following simple 
assumption as to the interactions between nucleons and spin 1 mesons.” 
“In the usual vector (or pseudovector) meson field, besides the usual vector 
(or pseudovector) source function in the equation (24) (or (25)), another one 
may be added. The added one is a vector (or pseudovector) source function 
which has the form 1/6e%?4,5 (or 1/6€*"°vs). The so-called tensor (or pseudo- 


tensor) source function is to be excluded.” 
As mentioned in §1 tensor and pseudotensor couplings give rise to divergent 


results, therefore we neglect them here. 


$4, Nuclear forces for the spin 1 meson field. 


As detailed calculations have been shown in the previous note (I), we shall 
refer to the results briefly. 


Lagrangians : 


1 h h* : 
Ly= -(s>) [y* (a-grad) — (grad $*-a)P+¢* * = ~ vb +20 | " 


(37) 


Iy=—(=) [ (curd A curd A) -(*4 +grad$)-( 24 +erad$)+#(4-A-$) | ’ 


(38) 
a it a 

Li=g,0e0" + & Evbon €™ Os (39) 
=—/¢*[ ¢,(a-A—9) +¢,'(0:A—p,d)] 9), (40) 

where co % ) a (-— p*¢, g* ay) 
and (ues ta = (Pop, —P*r').  (#=c=1). (41) 

Hamiltonians : 
H=HHyt+ Hy +H! + Air (42) 
ay ead 

Hy=(5-) | [¥*(a-grad$))—(grad $*-0)-+2Mprph] dV, (43) 


Hu=(5)\ | cat A curt A) +8 (A-A) + (jo) (die y+ Uh) ar, (44) 
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: (*\¢,(a- — div Il, )+¢/(¢-A—— p, div 1) | gal’, (45) 
x = 
Sey 1 o. = ° ie i } : o 
F143, = OQ (& (*h) ~—+ ( er ) ((h* 0,0) 2 


Ski P ' 
+( ee YW eh) Ph) + iH) Caron i (46) 


where 


OL 
(a =(°4 + crad é) : (47) 


The nuclear potential part in the Hamiltonian which is obtained from (42) 
by a canonical transformation” takes the following form: 


Fis, es (48) 
Ti 9 in M.a%)—] Af 
ae eer eay: -[(@ Fag (49) 
ore 4 > 4 2 
sea (Bl OE Rameau onl at ic aa 2 pa" 
" mL Lai) 
2iM Pein Bite Alles i 
+2 (0, 9,9 + 0.2 ,)( En, — Eon) |, (50) 
Pes 
Haste gt itis |i D) (a? Yea (Hd AP) 


Reese eae ® 40,2) (Ex, — Ex) | ; (51) 


It is to be noted that the direct (or contact) interaction terms vanish entirely 
in this result.* This means that in the Tomonaga-Schwinger formalisum the 
normal dependent terms of Hamiltonians cancel out in the second order nuclear 


potentials. 
Then, the expression (48) is transformed into the form which contains only 


© operator by the same method as adopted in §4 of ()®. Thus we find the 


nuclear potential 


W= -[-st+ ai" (6.6%) —o" | (6.6%) +0(r) oa 
oJ 


* In (1) the-signs of last terms of (51) and (46) were misprinted, 


The correct forms are-/7p, », a and-//,,, dir « 
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1 ! ),g@)— our. o°.P,)} 
+ (gh )sratt(e?-0?) 1G? -P,) + ( ) 


+( 2) ee—0){2 2-0) +00) Sof (0-0) 44 (Ga) 


2M r 
(52) 
where 
bed 1 -3e-nie —(o.0%)). 
Ce. > ui xr oe aot tes r c ) 
(53) 
or in the limit of WZ—> 
° 9 5 2 1 > cr 
W=-7[ -e2+g:"(@-0%)—2/"{ 5 (6-0) +07) Sa} |(G—) 
(54) 


where 7 is added to consider the symmetrical and charged theories, 
des (Ge ze) [2 or 2b OF Fare) 72: (55) 


This potential is nothing but the ordinary one which is obtained by mixing 
vector m2son with the vector coupling and the pseudovector meson with the 
pseudovector coupling. It is, however, remarked that the apperance of 1//* 
singularity in (54) is due to the inadequate approximation and by considering recoil 
effects of nucleons this singularity is removable, while the 1/7* singularity of the 
potential by the pseudovector meson by the old perturbation method came out 
from the direct interaction term, therefore it is impossible to avoid it by taking 
account of the recoil effect. Detailed account will be given in a future paper. 


§ 5. Diseussions 


One of the interests in this note is to see whether the spin 1 field 
is preferable to explain the nuclear phenomena in place of the conventional 
pseudoscalar meson theory or not. From the above calculation it is impossible 
to derive the decisive conclusion at the present stage, but we shall point out 
a few important results. 


(1) As regards the Coupling constants, it is the first problem to infer which 


is larger, (£5) ot 1 (£ i 
Arhc Arke Considering the anomalous magnetic moments of 


nucleons which have been evaluated by many authors, the estimated order of 
magnitude of them are ~ 4 in the vector theory, on the other hand, the pseudo- 


vector theory gives the reasonable ones. Hence, it seems to be. natural that we 
assume 


lal > ley’. 
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When the numerical difference betwecn them is large, we may give up the 
vector coupling in the spin 1 field and retain on!y the pseudovector one for the 
same reaSon as magnetic poles are excluded in the Maxwell-Dirac theory. As 
to magnetic moments themselves, it is necessary further investigations to obtain 
the correct results because of the ambiguity of meson currents.” 

(2) It is to be noted that, if we select the models which give rise to the 
correct signs of the quadrupole moment of the deuteron (i.c. positive sign) out 
of Kemmet’s four types of mesons with the exception of tensor couplings and 
the neutral theory, we get the next ones 


Pseudoscalar 


Pe nar festa fields. 


Charged and Symmetrical 


In fact in the latter case, it is easily confirmed in the expression (54). 

Finally, we are able to say that the spin 1 field with or without vcctor 
coupling is as remarkable as the pseudoscalar meson field. 

The author wishes to express his thanks to Prof. H. Yukawa and Assist. 
Prof. T. Inoue for their kind interests to this woik. 
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Errata (Vol. V, No. 1) 


The arguments of § 4 and §5 of (I) are incorrect. That is, the cross terms are dropped there. 
Terms of the form 


1 ; 
Ge (OS! AiO, Py + 62: Pe) —2i fd filoy-Py + a2 Ps)d} 
in two brackets of Eq. (75), and 
| Fo Si fi(G-Py+02-Pe)o 


in the brackets of Eqs. (76) and (78) should be added respectively. 
Ilowever, in the limit of A7-> c these vanish. 
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Letters to the Editor 


Notes on Fluctuation in Nuclear 
Evaporation Process 


K. Takayanagi and Y. Yamaguchi 


Department of Physies, Saitama University 
and University of City Osaka 


July 12, 1950 


In this previous paper, we have given a 
very brief discussion of fluctuation problem 
connected with nuclear evaporation processes. 
Thus in this note we shall give some com- 
plementary considerations about this problem. 

For brevity, we assume that each particle 
(j) evaporated from excited nucleus (with 
temperature 7’ and excitation energy X, 7’ 
=(X)) takes away a definite energy 


EM=E,4+V+f(X), f(X)~27, 


where E,=binding energy 
~ V =Coulomb barrier heigh 


for this 
particle. 
We will not treat here the fluctuation of the 
number of evaporated particles arising from 
the energy fluctuation of these particles (cf. I). 
Furthermore, we do not take account of the 
change of binding energies, and of mass or 
atomic number of nucleus. However, the 
possibility that the nucleus can emit various 
kinds of particles results in the fluctuation 
of the number of ejected particles of each 
kind. Let us consider this fluctuation, using 
the method given by Horning and Baumhof®. 

For simplicity, we treat the simplest case 
where the excited nucleus evaporates only 
protons and neutrons. It is easy to extend 
our treatment to more general cases, where 
the other particles may be emitted. 

Let p(z,v; X)dX be the probability that 
just ~ protons and y neutrons are emitted 
from the nucleus, when nuclear excitation 
falls into a interval (XY, ¥+dX) from. its 


initial value XY, in the course of evaporation. 
‘Then, if we put 


1 PrE+TpE 


1 khoc A Babee 2 
E™=E,+(X), 
E®=E,+V+4+f(X), 


and 


S(X)=f"(X)dX, 8, =SX), 


we find the following differential equation 


: In 
Ope, mie Pl, ¥3 S)— pple »-1:9) 
poe ple—l, v3 J @ 
Be AT ar ko oe 


In these equations I, and I, mean the 
neutron- and proton-width, respectively. If 
the ratio I’,/I’, does not depend on 
(otherwise one may take suitable average 
value of 7’,/I,), the solution of (2) under 
the initial condition 


P(r, vw; S)=dacd so 


is, as is well known, given by the double 
Poisson distribution : 


S,—-S Son Ww. (So—S) 
ar) Pie 
2’) 


If we do not interest in the number of 
ejected neutrons, we have only to use the 
margin distribution 


XS ple, v; S)=p(z, 8). 
From (2) we obtain 
a9h™ M=Px, Q-r(e-1, @,  @) 


where we introduce the new variable : 


1 
P(r, v3 SJ\=——-( 
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p(nV) 


04 en “N 


4 ° 
(/X=40 Mews 
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— x —— solution (Sa) 
ae Ope (ance with 
initial condition (6) 


_ Number of protons, % 


Fig. 1 


Comparison of the solution (5a) with the solution obtained by 


Ilorning-Baumhofi’s method, for the initial excitation: Xo=100 Mev. and 
the mass number 4=100. (We assume that #()=14+27, LM=8+27, 


[p=1 and [,,=1.5 in Mev.) 


If the ratio 7',/I, does not depend on X, 
one finds 
e=(x)=f oxex, 4 = £04 Em. 
0 q(X) I pt 
Similarly, for the margin distribution 


pv, S)= Spt, y, 8) 


we get 
Sav, R)=ply, R) =—plu-nls fk), 
In 
R= ToATs (4) 


Evidently, the Poisson distributions : 


pln, Q)=erw-o GO" 9,= Q(X), 


(5a) 
ply, Q)= eno Lo BI R,=R(X), 

(56) 
satisfy the equations (2) and (4). Thus, as 


for margin distributions, the expectation in I 
is confirmed. Furthermore, as is readily seen 
from (2), the distribution for m+v is also 
Poisson type. This distribution shows large 


fluctuation, though the fluctuation caused by 
the energy spectrum of individual particle is 
more important (cf. I). 

Now, p(z, Q) corresponds to px(Q) used 
by Horning and Baumhoff.” The solutions 
(5), however, belongs to different initial condi- 
tions than the ones adopted by them, i.e. 


p(x, Q)=0 for O>0,=Q(A0—aT,), 


T,=E,+V. (6) 
The mean value z(Q) or v(R) in our case is 
larger than the one obtained by Horning and 
Baumhof. This situation is explained by 
Fig. 1. 

Finally, we discuss the relations to the 
“mean cooling law”. If we approximate 
the equation (3) by the differential equation 

‘tp _ Op | 3/ 
0Q Ox’ Go 
with the initial condition 
Px ,Qo)=4(z), 
we get 
p(n, Q)=0(Q+7—-Q,) 
or 


Tine. (@) = Qo- Q. 
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This is just equal to the mean value obtained 
from (5): 


x(Q)= Simp(z, Q)=2.-@. 


Entirely the same result is obtained for 
neutrons. In general, similar results hold 
foreany kinds of evaporated particles. Thus, 
we believe that the contents of mean cooling 
law are considerably clarified. Namely, we 
may conclude as follows: The number of 
evaporated particles from the mean cooling 
law of differential type is the mean value. 
Actual number obeys Poisson distribution. 
Of course, these are valid only for our sim- 
plified case described above. It may, however, 
hoped that the main features of fluctuation 
would be roughly preserved in actual cases. 


1) Y. Fujimoto and Y. Yamaguchi, Prog. Theor. 
Phys. 4 (1949), 468; cited as I. 

2) W. Horning and L. Baumhoff, Phys. Rev. 
75 (1949), 370. 


On the Stars of Light Nuclei 


Y. Yamaguchi 
Department of Physics, 
University of City Osaka 

July 30, 1950 


Recently, Fermi’ has discussed the multiple 
meson production by a_ purely _ statistical 
method. He derived the distribution of 
number of produced mesons, assuming that 
possible final states occur with frequencies 
proportional to their statistical weights. Such 
a method could be used .to examine the 
disintegration modes of excited nuclei. In 
fact, as for heavy nuclei the statistical theory 
was successfully applied to various reactions”). 
However, the “stars” of light nuclei were 
not treated so far along this line. 
purpose of this note to discuss this problem. 

Firstly, let us consider the negative pion 
capture process in a C’?-nucleus. (The pion 
stars of heavy nuclei were discussed for- 


It is the. 


merly”).) If a pion is captured, its rest 
energy is dumped in the volume 2, which 
can be supposed to be of the order of nuclear 
volume (47/3)R?, R=nuclear radius, and 
contains seven neutrons and five protons in 
our case. More precisely, we may regard 2 
as an adjustable parameter. As was stated 
before, the probability of transition into the 
final state is assumed to be proportional to 
the volume S of phase space incorporated 
with the total energy 7’=ye*—B, where 
pe?=140 Mev is rest energy of a pion and 
B means the decrease of available energy 
due to binding energies. (We use Fermi’s 
notations, cf. ref. 1).) 7’ is roughly determined 
by the number vy, of alpha-particles existing 
in the final state : 

7T~60, 80 and 100 Mev 

for va=0, 1 and 2, respectively. 

For given vs, S is larger for increasing the 
number vy of particles in a final state. While 
for the smaller », the larger becomes &, 
since the increase of 7’ overcomes the decrease 
of y. Therefore, in our case the final state 
2 alpha-particles+3 neutrons+1 proton, i.e., 
the three prong-star, would have largest 
probability and prcbabilities of other states 
are smaller than this. 

Similar results are obtained for pion stars 
due to N or O nuclei. Thus we may conclude 
that the size distribution of pion stars of light 
nuclei shows a considerably sharp maximum 
at 3~4-prongs. This result is in agreement 
with the fact known from the’ recent photo- 
graphic experiments”. 

Next we shall examine the: reactions of 
light nuclei bombarded with high energy 
particles (< 150 Mev). As can be seen from 
the case of pion stars, even when 7’~150 
Mev the most probable final state would 
contain rather few prongs. If 7’ is small, 
the transitions to final states with larger » are 
completely forbidden or allowed with only 
very small probabilities. And the most 
probable number of prongs will not change 
drastically with 7. Furthermore, on account 
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of the nuclear transparency the available energy 
T of nuclear reaction shows very broad dis- 
tribution”. Thus we may expect that the 
prong distribution would be almost insensitive 
to the incident energy of bombarding particles. 

his fact was confirmed by Berkeley’s ex- 
periments”. 

These considerations ‘may be regarded as 
an evidence for what would be expected 
qualitatively from the so-called alpha-particle 


- model for a nucleus. 


1) &. Fermi, Prog. Theor. Phys. 5 (1950), 570. 
2) V. Weisskopf. Phys. Rev. 52 (1937), 295. 
3) Y. Fujimoto, S. Hayakawa and Y. Yamaguchi, 

Prog. Theor. Phys. 4 (1949), 575 and 576. 
4) D. H. Perkins, Phil. Mag. 40 (1949), 601. 
5) R. Serber, Phys. Rev. 72 (1947), 1114. 

M. L. Goldberger, Phys. Rev. 74 (1948), 1269. 
6) E. Gardner and V. Peterson, Phys. Rev. 75 

(1949), 364. 

W. Horning and L. Baumhoff, Phys. Rev. 

75 (1949), 370. 


On the Effect of 7-Ray to the 
Cosmic-Ray Underground 


S. Ogawa and S. Miura 


Institute of Theoretical Physics, 
Nagoya University 


August 6, 1920 


Recently, the production of m-meson by 
high energy y-rays (~300 Mev) is precisely 
analysed.” We have studied the effect of 
the high energy y-rays in the cosmic-ray 
phenomena. Our purpose is to obtain 
phenomenologically the character of the 7-7 
process in high energy region (~10%ev or 
more). Another aim is to make clear up the 
problem of the model in the cosmic-ray. 

In many experiments,” it has been pointed 
out that y-rays in cosmic-ray’ produce the 
penetrating particle. Tabin® has obtained 
the nuclear geometrical cross section as the 
production cross section of the penetrating 
particle by 7-rays. Considering that the cross 
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section of the ordinary pair creation is com- 
parable with the geometrical cross section of 
air nucleus, the value obtained by Tabin 
seems not so decisive—that is, if it were so, 
the ratio of the hard and soft-component in 
air would change appreciably. As- another 
effect of y-rays to the cosmic-ray phenomena, 
one of the authors has pointed out the pos- 
sibility of its effect to the positive excess of 
the hard component.” It is due to the 
asymmetry of the production of negative and 
positive z-meson by 7-rays. 

Now, s-meson and the soft component in 
eguilibrium with s-meson are the principal 
parts of the cosmic-ray underground. Nuclear 
component is of extremely small fraction. 
Evans and George® have obtained the follow- 
ing conclusion from the frequency of the 
star production in their experiment at 60m 
H.O depth underground. That is, ye-meson’s 
cross section for star production must be of 
the magnitude of 0.4 x 11-* cm? per nucleon. 
The result is inconsistent with two-meson 
theory. Another results against two-meson 
mode] is that obtained by Forrdé et al.® at 
the depth of 1000 mH,O underground. Ac- 
cording to them the intensity of the cosmic- 
ray measured by three coincidence decreases 
when 10cm thickness of lead is set between 
two counters, while it gradually - increases 
with the thickness of lead, indicating the 
maximum near 50cm. Rapid fall occurs only 
when 10cm Pb is put between upper two 
counters. 

In order to explain these results by the 
contribution of m-meson produced by 7-7 
process, we estimate the fraction Of m-meson 
relative to meson intensity at each depth 
underground. At first, we must inform about 
the interaction of z-meson with matter, espe- 
cially the collision mean free path. The 
information, however, is not so precise, though 
many experiments are attempted.” In spite 
of the lack of the knowledge, the interac- 
tion of z-meson with matter may be thought 
to be quite strong, as is seen from the con- 
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sideration correlating to the nucleon-nuclcon 
(or nucleus) interaction. Accordingly assum- 
ing the nuclear cascade process of = meson, 
we have calculated the transition curve of 
m-meson in matter based on the following 
assumption: i) a-meson with the kinetic 
energy more than its own rest energy acts 
cascade -likely.®) ii) z-meson loses its energy 
by a half, passing through one collision mean 
free path (7) and creates another z-meson. 
We remain / undetermined until the end of 
the numerical calculation. Thus, the obtained 
fraction of the intensity of z-meson whose 
kinetic energy is larger than the rest energy 
(x) in equilibrium with s-meson is presented 
in Table 1. 


Table 1. Fraction of z-meson relative to 


pe-meson 


Tn(>2e)/Ta= (bu N]-t-2)x 34 6.6 16x 10-2 
depth (in mlf,0) 60 200 1000 


¢o is the mean cross section of y-z% process per 
nucleus. VV and 4 are the number. of nucleon 
per gram and the atomic weight of the matter 
underground, respectively. We take 4=16. 


a) If the star ayent at 60mH,O depth is 
thus produced z-meson, we can estimate the 
value of %, independently on é. 


$)=4.8 x 10-*° cm?/nucleus. () 


We note that the mean energy of y-rays at 
this depth is about 10 Bev. Of course, the 
m-meson production by y-ray will be ac 
companied with the nucleon recoil, which 
may cause nuclear evaporation. The value 
(1) is also required when the star agent is 
y-ray itself. The further experiment will 
determine whether the agent is ionizing or 
nonionizing. We may expect comparable 
contribution due to these two effects. 

b) Taking the value (1), we obtain 
Tn} Ly= 0.21% 1077 at the depth of 1000 
mH,O. The z-meson in lead multiplicates 
and indicates the maximum in intensity at 
the .distance 7 or more apart. If we take 


the geometrical cross section of the nucleus 
as the z-nucleus interaction, 7~100g/cm*. 
This value of Z gives I,/I,~2% which is 
slightly small to explain the Forrd’s results. 
We remark, however, that the z-meson’s 
contribution may increase, if the multiple 
production of z-meson in z-nucleon interac- 
tion and z°-meson in the above process are 
taken into account. Considering that z-meson 
at 1000mH-,O is accompanied with shower 
(another z-meson or due to 2°-meson), we 
may be able to explain the decrease of the 
intensity caused by the insertion of 10cm 
Pb between upper counters. Because smaller 
energy part may be absorbed. We also 
remark that the value (1) is quite reasonable 
to explain the fact that the soft component 
in air shower is accompanied with the pene- 
trating particle by a few percent. De‘ailed 
paper will be reported near future. 

In conclusion, the authors should like to 
express their cordial thanks to Prof. S. Sakata 
for his kind interest in their work. 
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The Lateral and Angular Distribution 
of Cascade Showers 


J. Nishimura and K. Kamata 
Scientific Research Institute, Tokyo 
August 20, 1950 


The lateral and angular, spread of cascade 
showers are of great importance for the 
interpretation of air showers, and_ hitherto 
have been treated by various authors. The 
structure “function obtained by Moliére’ is 
the most available one at this time, but it 
was pointed out” that this function does not 
represent the structure of showers accurately. 

In order to obtain the exact structure 
function at any depth, we computed it ana- 
lytically under the assumption “ Approxima- 
tion B’’, since the inaccuracy of the Moliére 
function is mainly due to his numerical 
calculation based on the Arley’s approxima- 
tion. 

We started from the Landau’s*) equation 
and solved it analytically using the Mellin, 
Laplace and Fourier transformations. The 
course of our calculation, however, is some- 
what complicated and we shall give only the 
final resultst in this paper. 

The differential lateral structure function 
n,(Ey, FE, 7, t)2nrdr, which is the average 


_ number of electrons having the energy E 


atl 


and the deviation from the axis between r 
and r+dr at the thickness t, is given by 


m= — grief) dedpdg (tp) (-9) 
s I \ 2 7 Ep? \-P- 
B) Cx) Cer) (2), 
al 


where s, p, q are the complex variables, and 
the integration paths ( are running parallel 
to the imaginary axis. (p,q, ¢) is defined 
by, in the notation of Rossi-Greisen”, 


x My (p, qd t) 


} For the check of our calculation, we computed 
the total number of electrons from these distribu- 
tion functions, which gave the same resulls obtained 


by Snyder,) as it must be. 
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0 ee 
at +(A(s+2p+q)oo 


~Bo+2p-+)Cs+2p-+q)} wily 0 


=(£+){rG-0'9 


(oat (A(s+2p +4) + a) 


i(p—1, q, t) 


+(s+2p+q)q™i(p, q-I, t)} (2) 
with 
Mm, (0,0, t) =H, (se + Wye? ™, (3) 


Since this equation is similar to that ob- 
tained by Snyder and Bhabha-Chakrabarty”, 
it can be solved in the same way.* 

The total Jateral structure function I,(/o, 
0, r, t) is given by 

I,(Ey, 0, 7, )=) mak. (4) 
F 0 

Integration with respect to / and q yields 

an ES ; 
IT,(Bs, 0, 7, )= —qaa\) dtp r'(+y) 
x 1'(s+2p)mi(p, —2p—s, t) 


«(2) (ie) Ces) 


Similarly the differential and total angular 
structure functions are given by 


— \ \ J dsdpagl '(1+p) 


(©) 


Teo( Ko, a. 6, t)= a 
(p,q t) 


Ke) Ce), 


x gue 
ee '(1+p)l'(2p+s) 


<() ele) 
02) GO) Ge)” © 


IT, (-E,, 9, VA ore 

where M.(p, q, t) is defined by 
{3 O° 
ot’ 


— Bls+ p+ Newt 2p +a bmp, oO 


(6) 


XM2(p, —2p— 


“+ (Als-+2p-+q) +00) a+ (A(s-+2p-+q)oo 


ay A ate 
=(5 + a ){ pM? din q) t) 


+ (s+ 2pt+q)que(y, g-1, ts (8) 


* MLC y, 4) is given by pulting & ¢ (variable 
of the solution of eq. (2))—>90. 
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with 
m.(0, 0, =m, (0, 9, ¢). (9) 
For the comparison with the results ob- 
tained by others, we first computed the mean 
sjuare deviation and angle of deflection 
defined by 


fh rn ,2ardr 
ee ae ela (10) 
‘ n,2rrdr 
and 
| 0°n,2n0d0 
<0 sarge saasiena Ties (11) 
\ m.270d0 


Making use of our distribution functions 
(eqs. (1), (6)) and limiting ourselves to the 
shower maximum and E-+oo, we can obtain 
the Oe results : 


( °, ; 2. Ky’ 
Soa IN Taw et ). ( f oR ir) = 0-723(45), 
(10)’ 
x,  _ Mell, 0, NK)? _ kK 
<> a =F(0, 0B) = 9-569), 
(11)’ 


The results of va-ious authors are listed in 
Table I. 


Table I. Comparison of shower spread at high 
energies obtained by various authors. 


rar(E NG (radi. unit) eoar(2 a) (radian) 


Wuer-Wergeland) 0.074 0.33 
Moliere 0.835 0.6 
Belenky’) 0.94 = 
Janossy*") 0.724 0.570 
Roberg-Nordheim'™® 0.642 0.545 
Ours 0.723 0.569 


* These values are derived using the method 
of Roberg-Nordheim,' 


It must be noted that our values are very 
close to that obtained by Janossy which are 
the most accurate ones so far. 

The most remarkable results derived from 
equation (5) is the behaviour of the total 


Structure function near the core. Using the 


pole of I’(1+p)I"(2p+s), we can get the 
following results. 


I~? 


IT,~const. 


for s<2, (12) 
for s>2. (12)’ 

These results are similar to that obtained 
by Pomeranchuk'’” and Migdal’? but the 


definition of s is different from theirs, pro- 
vided E,4#00. In our case, s is defined by 


Ay’ (. e+ in( 2) +in(& 7 )= =a) Fee 
while their s is defined by 


ay’ (t+ in(~2) =0. (14) 
As 4/()~ —4» eqs. (18) (14) become 

smu fa(S)om()} ay 

*=t/in(~). (4y’ 


As is easily seen from these equations, s 
is slightly larger than theirs, and their results 
must be interpreted 4s the accurate ones only 
for E,=0o. Our structure function near the 
core is less steep compared with theirs. The 


.physical significance of this-result is that the 


core is composed of. the high energy particles 
having the spectrum deviated from the single 
power law when £ is finite. 

Detailed accounts and some applications 
to cosmic ray phenomena will be given in 
the subsequent paper. 

The authors wish :o express their apprecia- 
tion to Dr. Blatt, Messrs. Hayakawa and 
Fujimoto for their valuable advice throughout 
this work. 
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On the Equivalence of Pseudoscalar 
and Pseudovector Coupling Constants 
in Pseudoscalar Meson Theory 


M. Sugawara and Y. Ono. 


Physics Department, Hokkaido University, 
Sapporo 


September 3, 1950 


The problem of the equivalence of pseudo- 
scalar and pseudovector couplings in pseudo- 
scalar meson theory has been examined by 
several authors and especially Case’? proved 
that to the second order in the coupling 
constants the two couplings are equivalent 
for the calculation of nuclear forces and 
for many other simpler problems using the 
Schwinger-Tomonaga’s many-time formalism 
by carrying through the successive two 
canonical transformations upon the Schréd- 
inger equation in order to eliminate the 
divergence terms. In this short note we 
shall give a very simple proof that the two 
couplings are equivalent for the second order 
approximation. This proof consists in elimi- 
nating odd operators”? from Dirac equation 
for a nucleon in order to obtain Pauli type 
equation. This can be done by a canonical 
transformation on the Dirac Hamiltonian for 
a nucleon interacting with external meson 
field. The Schrédinger equation in the case 
of charged pseudoscalar meson field is given 
by 

ihOV /Ot= HY, (1) 
H=cp,o-p+03Me? +f (t +02U* +7-62U) 


+2(r,0-gradU* + r_o-gradU 
1 9U* 1 aU 


Tera eco 7 
with 
T4=(t,+717,)/2, r_=(t,—tr,)/2, 
x= ye /h 


where U is the meson field quantity, Jf the 
nucleon mass, 4 the meson mass, f the 
pseudoscalar coupling constant, g the pseudo- 
vector one and other notations are usual ones. 
If S is an Hermitian operator, then the 
transformation, 


y’ =e, (3) 
H’ =e He-*S—ihe(de-*5/0t), (4) 
leaves Eq. (1) in the Hamiltonian form: 
thov /Ot=H’'V’. (5) 
If we write the Hamiltonian in the form 
H=BMe?+ E+ O, (6) 


where E is an even operator and O is an 
add operator, the canonical transformation 
generated by the Hermitian operator’ 


S= —(i/2Mc?)BO (7) 
gives the new Hamiltonian” 
VES 2 B 2 
H’=BMe?+ E+ Mer? 


—er lO. [O, E]+ ine 

ik 00 B 

+ omer? oe + ont #1 

al ak 
~ MMe! 
which contains no odd terms of order (1/ Me’)’. 
Thus we can successively remove odd terms 
from the Hamiltonian to any desired order 
in (1/Me?). When the kinetic energy of a 
nucleon and the interaction energy with 
meson field are small compared with Mec’, 
above expansion may be hoped to converge 
sufficiently rapidly. It is very easy to apply 
this general procedure to the Hamiltonian 
(2) and if we remove odd terms of (1/Me’*)’, 
(1/ Me*)' and (1/Mc*)? from the Hamiltonian, 


the results are 


OB 4 veers ee 7 (8) 
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Hy =BMe?+ 37" 


+g forgtad(rs U*+r_U) 


ty big file,Ut+t-U) (9) 
and 
HH, =6Me?+ ppv 


+2 -ograd(r, U*+r_U) 


am”? g 
Biz? -plep,o-p, 7a 
grad (r,U*+r_U)]], (10) 


-for pure pseudoscalar and pseudovector coupl- 

ings respectively, where we employed static 
approximation (putting time derivatives to 
zero). The last term in (10) is small in non- 
relativistic case, the last term in (9) the 
higher order one. Considering that U itself 
contains coupling constants implicitly, we can 
see readily the two couplings are equivalent 
to f? or g® in non-relativistic case. If the 
following relation is valid ; 


(u/2M)f=g, a) 
the contributions to static non-relativistic 
nuclear forces are entirely the same for two 
couplings and the nuclear potential is for 
both couplings given by 


27 Oe Ope, Me,@ 
=i (La an Toes : 
(2 rag Pgs (a +grad™) 


a he a ot 
(Onde) — a "2 


» (2) 
which can be obtained most easily by the 
method of Green’s function. Eq. (11) was 
derived by Araki” by calculating directly the 
non-relativistic matrix element of + grad U, 
If we leave the time derivatives of field 
quantities, the equivalence condition is 
changed to the form g=(#/4M)f. In other 
types of meson fields the above method does 
not give any definite result. 


|a, 2 | 


1) K. M. Case, Phys. Rev. 76 (1949), 14. 
2) An odd operator in Dirac theory is a Iirac 


matrix which has only matrix elements con- 
necting upper and lower components of the 
wave function, while an even operator is 
one having no such matrix elements. 

3) L. L. Foldy and S. A. Wouthuysen, Phys. 
Rev. 78 (1950), 29. 

4) The Eq. (32) in Foldy and Wouthuysen’s 
paper was unfortunately misprinted. Our 
equation (8) is correct. 

5) G. Araki, Prog. Theor. Phys. 4 (1948), 193. 


On the D-Function of Non-local Fields 


Y. Ono and M. Sugawara 
Physics Department, Hokkaido University, 
Sapporo 
September 5, 1950 


Recently Yukawa” developed a non-local 
field theory with the aim of eliminating 
difficulties of ordinary local fields. Although 
the effective differences between local and 
non-local theories will appear only when the 
interaction between fields is introduced”, it 
may be interesting to study its physical 
meaning within the limits of free field. As 
is well known, the ambiguities of ordinary 
theories arise from the singularity of D-func- 
tion on the light cone, we tried therefore 
first to examine the form and behavior of 
D-function of non-local field. 

Adopting the notations used in Yukawa’s 
paper, the field operators for scalar theory 
are expressed as follows : 


U(Xus ru) = foo § Ak)" A) lew by) x 
x exp(ikyX") 18(ra+ey), 


U*( Xus ry) = fo § (dk) (QU) U* (ew’ sha’) X 

exp(iky’ X") 110(ru—0,’), qd) 

e 
where 
Ulku, Lm) = U(ky, 0 y)0 (kyk™ + x°*) x 
xX O(kyl*)d(7,0"— 2°), 
U* (key, Uy’) = U* (ky’, Uy’) (key +x") Xx 
x O(ky!l"’)8(7,/0"” — 2°). (2) 


(From “Eq. 
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These latter operators satisfy the following 
commutation relation ; 


LU(ky, by), U*(ey’, ty’) 


=~ (buh +2) -8(bylt) x 


x Blut — 2) 9 (en — ky!) 8b — ly’). (3) 


(1), (2) and (3), the following 
relation is derived, 


MEULY tp), U*(Xy’, Ty ‘I= 
— HB ry try! )-d(rur*— 2°) x 


x J (dk) EA O(kuk* +x°)d(kur") x 


exp{th.(A° —X"’)} 
=110(rypt+ry’)d(ryr” —A*)- D(X" —X"”, 7). 
(4) 
It is noticed, firstly, the new D-function 
depends only on the mean position difference 
X"—X’"’ and reduces to the ordinary one 


19,0. Gate 


X"’)= 


—S(aey i oltnk® +2°) x 


x exp{zk, (X*—X"’)}, (4’) 
when 4 tends to zero. Secondly, it is ze:o 
unless the infernal c»ordinates satisfy the 
relation ry=—ty’- 

Integrating on the k, variable, D’(X", r,) 
reduces to 
D'(X*, r,)= 5 dle 8 er — Kr) x 
x {exp(i(& X— KX”")) —exp(—i(k X— K-X"))} 
wheze 
Xp=(X, MN), =r, 75), by 
K=VeLe. 


=(k, ks), 


As this function is relativistically invariant, 
we can adopt the Lorentz coordinate system, 
in which r*=7,=0, |e|=r=24. In addition 
we take the direction of » as z-axis. Then 
we obtain 

@estars BAGH: ie 


—exp(i(aAx—KX"))} 


(&.){exp(i(k_X— AX") 


903 


1 1 f 
= 93) the dk, pl exp(t(eX— KX") 


—exp(—7(AX— KX"))} 


= + ( (iade 


sin (LR sin 0 cosg— KX") 


() 


G0 


(r 1 k) 


Qni 
--> te Io(&R sin 6) sin KX! 


2A std 0 
A Rsin@ 0(Rsin 6) 


Rsin 6x 


x | (A(R sin 0) sin KX4/K) dk 
0 : 


where 
R=|X|, r-X=rRceos), k-X=kRsinOcos¢. 


This corresponds to the ordinary expression 


D(X) =" (ba sin kR sin KX" 


0 


1 0 fos kRsin KX? 
"R OR IK 
0 


=4ni dk. (5’) 
Though the integral in (5) cannot be ex- 
pressed by an analytic function of R and _X%, 
it is remarkable that it depends on the angle 
between » and X. Especially when @=0, 
it no more depends on F#. This situation 
appears more clearly,* when we make the 
function corresponding to the d-function of 
local field. Namely, in the coordinate system 
in which r,=0, we obtain 
aD (Xn tw) La0= 9 J Ndkadk, x 

x {exp(ik_X) +exp(—ihk_X)} 

= ((27)*/ta)d(X)0(¥), (6) 
corresponding to the ordinary one 


ai DUN a) seo (2m)*/)0N)O(YIAZ), (67 


where 
Nee (elie 
Fy. (5) or (6) becomes infinite when we 


make 4 zero. This apparent difficulty can 
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be removed readily when combining .D-func- 
tion with d(ryr*—A°) (see Eq. (4)) and noting 
8 (ryr” — 22) [A= 0 (0 —A)/22° tends to 2d(7) as 
A>0. Eq. (6) does not coincide with Eq. 
(6’) as 4-+0, although the not-yet-integrated 
form (4) does, showing pathological natures 
of D-functions. We summarize as follows ; 
Field operators are commutable unless the 
internal coordinates are mutually antiparallel ; 
if they are antiparallel, the commutators are 
independent of mean position difference, 
when the internal coordinates and mean 
position difference have the same direction ; 
when they are perpendicular to each other, 
commutators are 6é-function type as usual. 
Above results are hardly understandable from 
the standpoint that the internal coordinate 
has its physical meaning as the extension of 
a particle. But all these explorations are 
still pending to the future clarification. 


1) H. Yukawa, Phys. Rev. 77 (1950), 219. 
2) Mz. Fierz, Phys. Rev. 78 (1950), 184. 


Effects of Nuclear Motion on the Fine 
and the Hyperfine Structure of 
Hydrogen 


T. Ishidzu 
Department of Dynamics, 1st Faculty of 
Engineering, Tokyo University 
September 6, 1950 


The 16-component equation for a proton 
and an electron with Breit’s interaction has 
been solved now by expansion in power 
series in a* successively (a=e*/c&), not in 
A(=m/M) as in my previous work reported”, 
and for general values of the total angular 
momentum J. The treatment is essentially 
the same as that of Bechert and Meixner”, 
but the process of ‘approximation included 
there has become clearer in this work and 


some refinement in their results has been 
obtained. 


The results have shown: The correction 
for the Rydberg energy and the fine structure 
of the hydrogen is given by the Bohr's 
simple reduced mass correction factor (1+)? 
as well as the Darwin's additional term which 
can be expressed, in units of me*, by —(Ba'/ 
8n') (1+)7%, in accord with every previous 
work? ), The result reported in my last 
letter for the energy of order f*a* in the 
same units was different from this, but it 
has been proved that some misunderstanding 
was involved in its reduction. 

The correction factor for the hyperfine 
structure is given by (1+f)~° for the two 
extreme hyperfine levels for a given 2: j=l 
+1/2, F=l+1 and j=/—1/2, F=/—1 as 
Breit and Meyerott” and others” have shown, 
but for the other two levels j=/+1/2, F=l, 
a different factor (1+ 8)-*? follows. The 
correction for the hyperfine separation of the 
ground state or more generally, of energy S 
state of the hydrogen atoms is not aZected 
by this asymmetric factor and given by the 
factor (1+8)-* or 1—38 as usual, since the 
energy expression for the S state lacks a 
term which leads to this asymmetric factor. 


) Tv.Ishidzu, Prog. Theon Phys. 5 (1950), 307. | 

2) K. Bechert u. J. Meixner, Ann. d. Phys. 22 
<1985), 525. 

3) C.G. Darwin, Phil. Mag. VI, 39 (1920), 537. 

4) G. Breit and G. E. Brown, Phys. Rev. 74 
(1948), 1278. 

5) G. Breit and R. E. Meyerott, Phys. Rev. 72 
(1947), 1023; ibid. 75 (1949), 1447. 

6) G. Breit, G. E. Brown, and G. B. Arfken,” 
Phys. Rev. 76 (1949), 1299. 


On the Nuclear Stars 


Y. Yamaguchi 
MS! 

Department of Physics, Osaka City University 
October 1, 1950 


Recently, the nuclear evaporations with 
higher excitation energy have been discussed 
by many authors'-®. Thus it seems useful 
to compare with them. 
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First of all, HLP, in which various theories 
are checked in comparison with huge collec- 
tion of experimental results, is the most 
important material in this field. In this paper 
it is concluded that as for the energy-tem- 
perature relation the Fermi-gas model is more 
preferable than to the liquid drop model. 
This conclusion is consistent with the em- 
pirical information on nuclear level densities”. 
Furthermore, if we assume that the smaller 
nuclear radius proposed by Serber et al.®, 
the so-called Bardeen’s remark”:’ seems 
unnecessary (while HLP used the large nuclear 
radius of Bethe)’. 

HLP, HB, FY, CP and LC are all based 
on Weisskopf’s evaporation theory”. On 
‘the contrary Goldberger® has adopted the 
Monte Carlo method throughout the reaction 
process induced by high energy neutrons and 
has tried to check the statistical evaporation 
theory. I think, however, that his check is 
somewhat questionable because he used the 
Monte Carlo method until regions of too 
low energy to accept its validity. 

Next, CP, based on the liquid drop model, 
may be rejected by reason of its wrong energy 
distribution of star prongs, as discussed in 
HLP. Also the level density formula used 
in CP has little experimental support”. They 
have taken into account the fluctuation of 
numbers of evaporated particles, and have 
tried to interpret the 90Mev neutron induced 
‘yeaction and negative meson stars. Never- 
theless, if we include the effect of nuclear 
transparency” ®):).9)."), which gives very 
broad distribution of initial excitation energy 
of the nucleus, we can hardly obtain such 
results as given in CP. The differences 
between light and heavy nucleus-stars?' 
can not be explained in CP. More refined 
treatments of meson stars are given by us’? 
(as for light nucleus-stars, sce ref.))), 

HB and CP overlooked the change of 
binding energies in the course of evaporation 
process. This effect was included in FY and 
LC. FY has a defect that it is used rather 


crude approximations for various quantities. 
The treatment in LC is the most accurate 
and complete, though in this paper the 
fluctuation of numbers of emitted particles 
is not considered. The elementary considera- 
tions of this fluctuation are given in FY and 
ref. 16). 

In highly excited nuclei, there exist some 
special circumstances (e.g., the thermal ex- 
pansion’, depression of Coulomb _ barrier 
height’®, change of nuclear surface tension’, 
etc.), and these make the treatment of very 
large stars more complicated. 

In short, the present stage of nuclear 
evaporation theory may he compared with 
that of the cascade theory in cosmic-ray 
physics, and it is most desirable that more 
accurate calculations which take account of 
all possible effect and use empirically well 
determined quantities (binding energies, etc.), 
should be performed. 

In the star-problem the theory of light 
nucleus-stars is particularly poor. Heckrotte 
and Wolff? analyzed the carbon reactions 
on the basis of Weisskopf’s theory. But the 
nuclei such as C, N or O are too light for 
the statistical theory to be valid, although 
the tendency of prong-number distribution’? 
is explained from the statistical cansidera- 
tions’. 

Recently the striking differences between 
n-p and p-p scattering at high energies have 
been found», Thus the calculations on 
nuclear transparency of high energy nucleons 
(refs. 2), 6), 8) and 12) ; these are based on 
old ‘experimental results) should be revised. 
Since the n-n force is remained unknown at 
present, it is deeply desired. to analyze ex- 
perimentally and theoretically, ¢.g., n-d scat- 
tering at higher energy. 

The author wish to thank Mr. Fujimoto 
for his helpful discussions. 
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On the Positron Theory of Vacuum 


Y. Katayama 
Department of Physics, Kyoto University 
May 24, 1950 


The problems of photon self-energy and 
gauge invariancy are the appropriate examples 
of the formal difficulties in * the present 
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quantum field theory and we have no satis- 
factory answer notwithstanding the various 
methods of remedy proposed by many authors, 


such as the mixed theory, the regulator 
method and other technical methods of cal- 
culations. 


In these circumstances, we have a litile 
opportunity to re-examine the positron theory 
of vacuum analysed by Dirac’? and Heisen- 
berg? et al. in some years ago and some 
reflections on this problem have been pro- 
posed by a few authors in our country.» 

We also treat this theory from the tho 
roughly covariant view point developed by 
Tomonaga and Schwinger and then clarify 
the situation of these problems. 

At first, 
under the influences of the external electro- 


magnetic field satisfying the following Dirac’s- 


equation 
[ru (Ou—te Ay) +m] g(z)=0, qd) 


and then construct the vacuum expectation 
value 


<[d 2’), $@")] >.=—-SM@’, 2"), @) 
by means of S®(z’,2’") given by some 
integral equations. If there would exist no 
real electron field, this function is also given 
by the power series of e, that is 


De’ —2 rp A ple) “Ses 2043) 
2ozal 


TpAp (xa) San—2""). (3) 


Following Dirac’s program, one can get 


we consider the electron field — 


the physically significant density matrix R@ 


by means of subtraction of any singular term, 
RO’, "= SOG’, 2!) —-UM(2’, 2”). 
(4) 


And we must determine U™ by the re-_ 


quirements of the equation of motion and 
the conservation law. Although it will be 
more difficult to find this subtration terms 


f 
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rigorously than in Heisenberg’s first treatment, 
we can easily get the first several leading 
terms as follows : 


he s—5)= ula+ 9° na) SE) 


bse \ dal (dp) arms 


ce r)* 


Trtv iheDs —m) Fy, @— S + za) 


+(5 — 8a) ryip, P(e 5 +20) 
OF waa +20). 
Oxy | 


0’ (p? +m? *)u(: 5 te 2 ae (5) 


+2a(a—1)rp 


a 


Zz zZ 
where a =a2to, a’ =2—- and u(a’, xv") 
Pd 


are given by the integral equation : 


uta’, a )=1tiel Ay (x Vag (a’’’, a’’)dx,/"" 

(6) 
Expanding u(2’, #’”) by coupling constant-e, 
we have 


u(2+ 5 2-5) 
see 4 (o- 
eee zd, cae MT + Dy, Ab nZa** Zu 


2 
(6") 


and for the case of the product of this func- 
tion by the field independent some functions, 
this is written approximately 


spr fe n tdb,- 


n= 24 


u(a+ oe oy 
=exp| icf Ane — $+ 2a)daz,|. (7) 
That is, the first term of (5) becomes 


exp| ie \ “Andaz, |8C2) 


1 , 1 
a (dpye”[ira(pa —e\ A ,du) —m | 


i| 
xa} (p—ef Adu)? +n? | (8) 


which coincides with the first term of Heisen- 
berg’s calculations and the subtraction terms 
of Umezawa and Kawabe’s consideration.” 
However, for the other subtraction terms we 
can not modify so simply and the physical 
meaning of these terms are not so clear than 
the first. 

According to this subtraction procedure 
(5), we have the gauge invariant results and 
no charge renormalization in the usual co- 
variant theory. And the non-gauge terms, 
such as the photon self-energy, those of the 
y-decay of neutral scalar, pseudo-scalar and 
-vector meson and the photon-photon scatter- 
ing are given correctly in the subtraction 
terms as S(a’, 2”), that is, 


Mas [U™ (a,x) rp] Ap 


= x) (dp) [6(p?+m*) + ape’ (p?-+m?) | Ap* 


Ee eee 
8x (log m 1)J,A;, 


i [UM (a, 2)]% ey 


Se 1) 
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and 


PSU Ce, %) pels U,= 7oneg An A pU,. (9) 


Detailed account will appear in this journal. 
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Dr. YOSHIO NISHINA, His Sixtieth Birthday 


Dr. Yoshio Nishina’s sixtieth birthday will come round on December 6th, 1950. On this 
happy occasion, I, as one of his old pupils, should like, by way of expressing my hearty 
congratulation, to celebrate not only what he has done in physics, but also the great and 
important part he has played and will still play in the development of general culture and 


_science of this country. 


When he graduated from Tokyo Imperial University in 1918, he intended to devote him- 
self to the study of electric engineering. But after he entered the Institute of Physical and 
Chemical Research here, he came to be interested in pure physics, and so during his long 
stay in Europe from 1921 to 1928, too, he gave himself wholly up to that study. 

Does it not show that he has a mind both very brilliant and flexible that Dr. Nishina 
who had been educated to be an electric engineer proved to be an excellent scholar in pure 
physics? His research works cover a very wide scope, extending over not only theoretical 
physics and experimental physics, but even biology, agriculture, and medical science. Even 
within the limit of physics, his research activities cover almost all its branches such as 
quantum theory, theory of elementary particles, X-rays, atomic nucleus, and cosmic rays,— 
the result of his extraordinary enthusiasm and energy to carry out research works. 

It is needless as well as impossible to enumerate here one by one all the results of these 
many-sided researches of his, for he has done a very important and significant work which 
by itself will be enough to immortalize his name. 

And that is his study of the scattering of hard X-rays or y-rays by a free electron. 
During his stay in Europe, Dr. Nishina made a theoretical study of this problem with Dr. 
Klein, a Swedish theorist, and evolved from it the famous Klein-Nishina ‘formula. How 
widely this formula is being used, how important a part it is now playing, in the field of 
nuclear and cosmic ray physics, is a well-known fact to anyone who has any knowledge of 
the present-day physics. 

Apart from these excellent research achievements, his existence has a special significance 
to us physicists of Japan, firstly because it was no one but he that brought the seeds of new 
physics into this country and secondly because deep impressions and strong stimulations were 
given to young students and researchers by him through his character and his enthusiasm at 
once persistent and energetic. If it were not for ‘him, the younger generation in every 
branch of science in Japan would have been quite at a loss. 

Dr. Nishina, who had spent his days in Europe with Lord Rutherford at Cavendish 


_ Laboratory, and with Dr. Bohr at Copenhagen, had grown with the growth of new physics, 


absorbing a fresh and active atmosphere of its centres, and personally going through new 
methods of its study. So, when he came home later, he tried to introduce and transplant 
these methods and atmosphere into this country. How pleasant and exciting it was to us 
younger people to listen to his stories of his own experiences at those centres of new science 
in Europe! How deeply impressed and influenced we were by his zealous academic interest 
through his stories and deeds! 

He did not only encourage younger people morally, but, as soon as Nishina Laboratory 
was made in 1931 in the Institute of Physical and Chemical Research after his return from 
Europe, he set about his great work of making this laboratory a research centre of sblospic 
nucleus and cosmic rays. We can well. imagine how hard it must have been to realize that 


ii 4 
which was economically poor as well as lacking in 


general sympathy and understanding towards pure studies. A laboratory ae gs an — 
in new physics, whether it be nuclear physics or cosmic ray physics, is a la ee ory a 
very large scale which is furnished with a monster of a crane above our head and a giganti 
. a few hundred tons, and which has an electric 


high-tension apparatus weighing as much as ; : ; 
source and a machine shop attached to it besides—in short, it looks like a factory, so to 


speak. How laborious a task it must have been, and how painstakingly he a have work 

to enlighten the people and make them understand that a research in pure p ea “aa 
such a big “factory.” But at last his earnest campaign of enlightenment succeed , a = 
nuclear research laboratory was made in the Institute of Physical and Chemical Resear 

two cyclotrons were set there in 1935. The result of his campaign was not limited to that 
only, but some universities came to have their own “factories” of nuclear physics furnished 
with such majestic apparatuses as cyclotrons and so forth, for people at large were pursuaded 
by him and began at last to perceive the importance of the study of atomic nucleus. 

Here in his Nuclear Research Laboratory, Dr. Nishina wanted to carry out his great plan 
of research, not only in pure nuclear physics, but also in its wide application to biology, to 
agriculture, and even to medical science. Here, however, the Pacific War broke out, and it 
became difficult for him to continue his research, much more so to carry out his new plans 
and new ideas. Thus his grand conception collapsed. : 

Since the end of the War he has been and is, as the president of the Scientific Research 
Institute that succeeded the former Physical and Chemical Research Institute, making efforts 
to carry out its aim by good management though under extraordinarily difficult conditions. 
Thus he has turned now his elastic and versatile ability and untiring energy from his 
academic life to a new life of an administrator of a research institute. And he is, I hope, 
_ competent enough to accomplish the difficult task which is, in our post-war Japan, even more 
important and difficult than research itself. 

He is now also endeavouring, as vice-president of Japan Science Council, to contribute 
to science administration and to promote international co-operation in science. We now only 
wish and hope that the great scientist of world-wide fame who is also the father of new 
physics of Japan will be as healthy and as active as ever, because we expect much from him 
in these hard times of Japan and the world. ' 


great plan in this country of those days, 


Sin-itiro ToMONAGA 


Noy. 20th, 1950 


Obituary Notice : 


In the cource of the proof of this number, on Wednesday eyening, January 10, 1951, 


Dr. Y. Nishina was stricken by a cancer of the liver. Death came a few months after going 


to hospital in Tokyo. Here we offer our cordial condolence. 


Jan. 11, 1951 M. K. 


